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SUMMARY 


A computer program, the Propeller/Nacelle Aerodynamic Performance 
Prediction Analysis (PANPER), was developed for the prediction and analysis of 
the performance and airflow of propeller-nacelle configurations operating over 
a forward speed range inclusive of high speed flight typical of recent prop- 
fan designs. A propeller lifting line, wake program was combined with a com- 
pressible, viscous center body interaction program, originally developed for 
diffusers, to compute the propeller-nacelle flow field, blade loading distri- 
bution, propeller performance, and the nacelle forebody pressure and viscous 
drag distributions. The computer analysis is applicable to single and coaxial 
counter-rotating propellers. The blade geometries can include spanwise varia- 
tions in sweep, droop, taper, thickness, and airfoil section type. In the 
coaxial mode of operation the analysis can treat both equal and unequal blade 
number and rotational speeds on the propeller disks. The nacelle portion of 
the analysis can treat both free air and tunnel wall configurations including 
wall bleed. 

The analysis was applied to many different sets of flight conditions 
using selected aerodynamic modeling options. The influence of different pro- 
peller-nacelle-tunnel wall configurations was studied. Comparisons with 
available test data for both single and coaxial propeller configurations are 
presented along with a discussion of the results. 
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INTRODUCTION 


Brief History of the Problem 

The recognition by NASA and Hamilton Standard in 1974 of the potential 
for improved fuel efficiency of a high speed propeller relative to the conven- 
tional gas-turbine engine, for aircraft which cruise at high subsonic speeds, 
has established a renewed interest in this propulsion concept. This interest 
has resulted in a major participation in the technical evaluation of this 
concept by the Hamilton Standard Division (HSD) of the United Technologies 
Corporation (reference 1). A photograph of a conceptual model of the high 
speed propeller (Prop-Fan) is presented in Figure 1. In initial model tests 
conducted in a wind tunnel at the United Technologies Research Center (UTRC), 
under a NASA contract to Hamilton Standard, sufficiently high values of 
propeller efficiency were measured to encourage further interest in the con- 
cept (reference 2). The requirement for a computer analysis to assist in the 
design of the high speed propeller and the spinner-nacelle was recognized, and 
the task for the development of such an analysis was awarded to UTRC. Since 
these initial model tests, several different model designs have now been 
successfully tested by Hamilton Standard at UTRC and at NASA Lewis Research 
Center. The models tested have varied from blades of basically straight 
design to blade designs with severe blade sweep. Thus the requirement that 
the analysis be able to handle a high range of blade geometry designs is well 
established. 

Due to the aerodynamic complexity associated with the high speed propel- 
ler, previously existing analyses had clearly identifiable limitations. 
Development of a propeller-nacelle system for operation at high subsonic or 
transonic speeds resulted in several problems due to the transonic flow regime 
and the contemplated propeller design. Of particular concern were the 
problems associated with calculating the correct local inflow conditions at 
the high speed propeller blades due to the interference effect of the nacelle 
and the strong compressibility effects on the airfoil characteristics and 
induced inflow calculations. The severe blade geometry designs also have a 
significant impact on the induced velocity calculations for the propeller 
(reference 3). Figure 2 illustrates the large variation of blade geometry in 
a fixed coordinate system with changes in blade pitch angle which must be 
handled correctly. It is recognized that the high speed operating conditions 
and geometric configurations of the Prop-Fan compromise many of the assump- 
tions associated with the methodology used for conventional, low speed propeller 
designs. Analyses based on higher level modeling assumptions (full potential, 
Euler, and Navier-Stokes) are desirable and have been developed to some degree. 
However, these analyses are computationally expensive to run. As a result, 
the propeller analysis to be described was developed to model the propeller 
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induced flow field using a computationally efficient lifting-line theory with 
ad-hoc aerodynamic modeling features to approximate the influence of the 
transonic operating conditions and three-dimensionality of the Prop-Fan 
design. 

In addition to the difficult propeller aerodynamic problem, the 
propeller-nacelle flow field is expected to be influenced by viscous effects. 
Since friction losses due to the nacelle boundary layer will be considerably 
higher at transonic speeds than at subsonic speeds, a reliable propeller- 
nacelle flow field prediction method should have the ability to account for 
viscous phenomena. 

A procedure, which includes viscous effects, has been developed at UTRC 
for predicting circumferentially averaged flow in compressor stages 
(references 4 and 5). The procedure solves a set of viscous flow equations in 
the flow region of interest and thereby simultaneously accounts for blade 
effects and endwall effects. (In a propeller-nacelle field, the nacelle is 
analogous to the compressor endwall.) 

For application to the high speed propeller problem, the above analysis 
for the blade-nacelle interaction was expanded to include blade and wake 
induced aerodynamic effects through the coupling with a lifting line propeller 
model and an appropriate propeller wake model. This analysis was also modi- 
fied to treat the propeller nacelle problem under free flight conditions or 
wind tunnel conditions including wall bleed by using appropriate boundary 
conditions. 

The United Technologies Research Center, along with the Hamilton Standard 
and Sikorsky Aircraft Divisions of UTC, has been actively engaged in the 
development and evaluation of computer analyses for predicting the induced 
effects of propellers and single and dual counter-rotating helicopter rotors. 
In particular, a propeller performance analysis was developed at Hamilton 
Standard based on an earlier UTRC analysis, which has recently been adapted to 
high speed propellers. With this background, existing computer methods for 
lifting line and wake modeling techniques were combined into a new propeller 
analysis to handle the special requirements of the high speed propeller, and 
the propeller analysis was combined with the blade-nacelle analysis to provide 
a self-consistent propeller-nacelle performance prediction analysis. 
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Technical Background 


The usual method of approaching the propeller design problem is based 
upon lifting line-vortex theory or momentum theory approaches to determine the 
blade loading distribution (references 6 through 19), described briefly in the 
following section. Although these approaches, as commonly used, can be 
applied successfully for certain flow situations, they do suffer two serious 
drawbacks. First, either these approaches do not account for the interaction 
between the nacelle and propeller flow field at all, or they account for the 
interaction via empirical correction factors. Secondly, usual propeller 
analyses ignore cascade effects along the inner portion of the propeller disk 
which may be important for propellers having a large number of blades. 
Considering the interaction problem first, it should be noted that the blade 
nacelle interaction flow represents a highly complex nonlinear process in 
which the nacelle and blades mutually interact to a significant extent. Both 
the physical contour of the nacelle and the viscous boundary layer developing 
on the nacelle surface can affect the propeller-nacelle performance. For 
example, the nacelle shape and viscous effects can modify the propeller wake 
development thereby modifying the induced velocity at the propeller flow. In 
addition, the nacelle can modify the pressure distribution and the resulting 
velocities through the propeller blade row, particularly in the inner region 
of the blade disk thereby again affecting propeller loading and performance. 

Similarly the flow about the nacelle is affected by the presence of the 
propeller blades. The nacelle drag (viscous drag due to the nacelle boundary 
layer and drag due to the pressure field) is expected to be influenced by the 
presence of the blades. Although the propeller-nacelle interaction problem 
has not been addressed in detail, in propeller calculations to date, similar 
interactions have been addressed in compressor design analyses. In partic- 
ular, a computer program developed at UTRC (the UTRC ADD code), which calcu- 
lates the circumferentially averaged flow through a compressor stage and which 
includes the effects of endwalls, has been used successfully for the com- 
pressor through flow problem (references 4 and 5). The code solves the entire 
flow field at once and does not require an iteration between the viscous wall 
flow and the nominally inviscid flow in the center of the stage. This same 
procedure also takes into account the cascade effects which occur in a 
compressor stage and which as previously mentioned, are expected to occur in 
the inner portion of a propeller disk. Therefore, based upon previous exper- 
ience with this compressor stage code, the propeller-nacelle interaction 
problem was considered to be a reasonable application for this code. 
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Propeller Analysis 


The field of low speed propeller and rotor aerodynamics is abound with 
literature describing the need for and use of variable inflow performance 
analyses and the use of vortex wake models (e.g., references 5 through 33). 
UTRC has recognized this requirement and has helped develop these. analyses to 
a high level of sophistication. Today most all rotor and propeller manufac- 
turers use some type of a variable inflow model based on wake modeling for 
accurate performance predictions. 

Goldstein (reference 6), Theodorsen (reference 7), and other early inves- 
tigators realized the need for accurate descriptions of the propeller inflow 
distributions for performance predictions based on consideration of the 
propeller wake. Early lifting line performance analyses (references 6 through 
8) used helical, uncontracted Goldstein type wake models (classical undis- 
torted wake models) to obtain variable inflow distributions for performance 
predictions. These analyses gave better predictions than simple momentum and 
blade element analyses, but for many applications they still left a large gap 
between the analytical and actual test results. Distorted wake geometries 
which more accurately represent the actual physical wake geometry were 
required. 

For static thrust conditions, investigators attempted to compute wake 
distortions numerically using free wake interaction analyses where the wake 
was free to move under its own induced influence (references 14 and 15). 

These analyses, although more accurate than undistorted wake analyses, are too 
time consuming in terms of computer costs to be used in standard design pro- 
cedures. Therefore, several investigators studied the results of these calcu- 
lations and flow visualization studies (references 21 through 27) in an 
attempt to develop accurate empirical wake description formulae. In parti- 
cular, pioneering efforts by Gray (reference 25), Landgrebe (reference 26), 
and Ladden (reference 17) experimentally defined the wake geometry character- 
istics of hovering rotors and statically thrusting propellers. Their results 
were generalized to give simple analytical expressions for the distorted wake 
geometries based on the operating conditions and rotor parameters. The gen- 
eralized wake equations of Landgrebe and Ladden have now been used for several 
years by both the helicopter and propeller manufacturers in lifting line 
performance analyses for static thrust conditions with some modifications for 
each application. A sample of the distorted wake and classical wake geo- 
metries is represented in figure 3. 

At low or moderate propeller forward flight speeds, the need for induced 
inflow distributions and vortex wake modeling is also of importance. At these 
flight speeds a much larger portion of the propeller wake transport velocity 
comprises the forward flight speed. Because of this, axial wake distortions 
do not play as important a role in the wake geometry's influence on the inflow 
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distribution at the propeller blades as in the static thrust case. However, 
propeller self-induced radial contraction effects on the wake geometry are 
still of importance for accurate descriptions of the inflow distributions. 
Although normally neglected in previous operational propeller analyses for 
moderate speed flight , it is recognized that the nacelle does displace the 
propeller wake and alter the flowfield at the propeller disk. However, this 
effect is stronger in the region of the blade root and not near the blade tip 
where it would be most influential on the blade loading. 

At high speeds, accurate performance calculations also require the pre- 
diction of detailed inflow distributions which are dependent on vortex wake 
modeling. The propeller self-induced wake distortions of lower flight speeds 
are small at high flight speeds and the wake transport velocity is dominated 
by the propeller forward flight speed. At these speeds other factors may 
affect the wake geometry more than the induced effects of the wake on itself. 
Such influences as the effects of nacelle body blockage (or accelerations) can 
distort the wake geometry, particularly near the inboard regions of the 
propeller blades. Also, the effect of compressibility can alter the wake's 
influence at the propeller disk. The influence of the propeller wake and the 
resulting inflow distribution at the propeller blades is thus important for 
accurate performance predictions. 


Nacelle Analysis 


In previous compressor analyses for the circumferential average flow, two 
broad assumptions are generally made. The first assumption, or quasi-steady 
assumption, states that the flow is steady (constant in time) when examined in 
a coordinate system rotating with the angular velocity of the propeller. 

Thus, at each radial position, temporal and circumferential averages of the 
flow are the same. This assumption is equally applicable to a compressor or 
propeller for axial freestream flow. The second assumption states that the 
flow is primarily inviscid and implies viscous forces are small compared to 
blade induced forces on the airflow. This assumption is also generally 
applicable to either the compressor or propeller problem (except in the 
immediate vicinity of the compressor endwall or propeller nacelle). Based on 
these assumptions, the circumferential averaged equations of motion may be 
derived. In the past, two methods have been developed to solve these equa- 
tions: the matrix through flow analysis (references 34 and 35), and the 

streamline curvature method (references 36 and 37). Both procedures solve for 
a "core" flow, which is inviscid, and account for the presence of blade forces 
by arbitrarily specifying work and loss distributions at the blade trailing 
edge. Once the core flow is calculated the annulus wall boundary layer 
(nacelle boundary layer) may be calculated from the known pressure distribu- 
tion. If the growth of the annulus boundary layer significantly affects the 
pressure distribution, a solution of the strong interaction problem exists and 
requires an iteration between two flow fields: the core flow, and the 

boundary layer flow. 

Recently, this strong-interaction problem was solved at UTRC in a more 
direct manner without iteration between flow fields (references 37 and 38). 
This viscous method, which forms the basis for the UTRC ADD code, can be 
described in the following manner. First an orthogonal coordinate system is 
constructed from a plane potential flow solution such that the stream function 
forms the coordinate normal to the wall and the velocity potential forms the 
coordinate tangent to the wall. The equations of motion may be greatly 
simplified by assuming that the velocity normal to the potential flow 
streamlines is small compared to the streamwise velocity. Such an assumption 
is quite reasonable as the potential flow streamlines are expected to be a 
good approximation to the actual streamlines. This procedure reduces the 
governing viscous flow equations to a parabolic system of partial differential 
equations which can be solved by a forward marching numerical integration 
procedure. Furthermore, it can be shown that the resulting solution has the 
same order of accuracy for viscous flows as the streamline curvature method 
has for inviscid flows. Thus the viscous-flow strong interaction with the 
wall (nacelle) boundary layer is solved without the need for iterations 
between different flow fields; i.e., the inviscid core flow away from the 
nacelle and the nacelle boundary layer flow. 
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Combined Analysis 


The approach taken to create the combined propeller-nacelle analysis was 
to use a modified version of the original UTRC ADD code which includes the 
required refinements necessary for the nacelle-propeller problem. The calcu- 
lation of the blade forces was handled by creating as a portion of the 
combined analysis an expanded propeller lifting line solution procedure based 
on existing UTC propeller and helicopter performance analyses. This combined 
analysis consists of two portions (propeller and nacelle) which compute the 
combined propeller-nacelle performance by interfacing the required transfer of 
internal data (flow field properties, wake geometry displacements, and blade 
forces) which will influence the respective solution procedures from the two 
different portions of the analysis. This procedure allows the combined 
analysis to have the ability to also calculate isolated propeller or isolated 
nacelle performance for comparisons with the combined propeller-nacelle con- 
figuration performance predictions. A simplified diagram of the program task 
structure of the combined analysis is presented in figure 3. The details of 
the technical approach of each portion of the analysis and the coupling pro- 
cedure of the combined analysis are presented in the following sections. 
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TECHNICAL APPROACH - PROPELLER 


Overview 

In the following sections, a method of analysis for the predictions of 
the integrated and local spanwise blade air loading for the high speed 
propeller is presented. The analysis presented for the high speed propeller 
portion of the combined propeller-nacelle code is a logical extension of the 
development of the Prescribed Wake Method of Landgrebe (reference 26) for 
hovering rotors which was modified and applied to statically thrusting 
propellers by Landgrebe and then Ladden (reference 17). Although this anal- 
ysis is directed towards the high speed propeller flight regime, no restric- 
tions have been made in the assumptions which would limit it to only high 
speed flight. Briefly the Prescribed Wake method is derived utilizing lifting 
line theory, and incorporates a wake model consisting of a finite number of 
trailing vortex filaments. The trajectories and positioning of these fila- 
ments must be prescribed. Figure 4 is an illustration of two types of 
prescribed wake models used for hovering rotors (similar representations are 
used for propellers). Once the position of the wake is fixed, a matrix solu- 
tion is solved utilizing the Kutta-Joukowski and Biot-Savart relations, and 
the airfoil lift properties. This solution yields the blade circulation dis- 
tribution and the corresponding induced velocities. The blade element 
velocity diagram is then constructed, and with the use of two-dimensional 
airfoil data corrected for three-dimensional tip effects, the airloading 
distribution (lift and drag) is obtained. The total thrust and power are then 
established through spanwise integration. 

The use of the propeller analysis in conjunction with the nacelle code is 
handled in a coupled manner, whereby the nacelle's influence on the propeller 
is incorporated through modifications to the noninduced inflow at the 
propeller disk and to the wake model. The theory of the analysis as outlined 
in the following sections does not make direct reference to the nacelle's 
influence on the velocites. It is thus understood that the noninduced 
velocity terms include the nacelle's influence in the following propeller 
related sections. 
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Coordinate Systems 


Within the lifting line propeller analysis there are three coordinate 
systems used. Two of these coordinate systems are fixed to the same origin. 
These are a cartesian and a cylindrical coordinate system whose origin is 
defined by the intersection of the propeller's axis of rotation and the blade 
pitch axis of rotation (which are always assumed to intersect). The axial 
coordinate directions for both of these coordinate systems are identical and 
coincident with the axis of rotation of the propeller, positive in the direc- 
tion of flight. The blade pitch axis of rotation is coincident with the x 
axis in the cartesian coordinate system. The relationships between the x, y 
plane and the r, <t> plane are the conventional relationships for right handed 
coordinate systems. Figure 5 illustrates these two coordinate systems. The 
cartesian coordinate system is used to define the input blade geometry, while 
the cylindrical coordinate system is used for computations involving the wake 
and blade geometry because of the symmetry of the problem in this coordinate 
system. 

The third coordinate system is called the blade element coordinate 
system. The local origins for this system are defined by a reference line 
running along the blade in the spanwise direction. Generally the quarter 
chord line, which is not necessarily straight, is used for the reference line. 
This is a locally orthogonal coordinate system whose unit directions are 
defined in the local spanwise, chordwise and normalwise directions. The local 
spanwise direction is defined as positive in the direction from the root to 
the tip along the local tangent to the blade reference line noted above. The 
local chordwise direction is defined to be perpendicular to the local spanwise 
direction and in the plane of the zero angle of attack of the local airfoil 
section, positive in the direction of the trailing edge to the leading edge. 
The local positive normalwise direction is defined as the direction of the 
local vector cross product of the spanwise to the chordwise direction 
vectors . 

Figure 6 is an illustration of this coordinate system. The blade element 
coordinate system is the computational coordinate system used in the lifting 
line analysis solution. The reasons for using this coordinate system will be 
explained in the following sections. The transformation relationships for the 
various quantities used in the analysis between the various coordinate systems 
are shown in Appendix A. 



Propeller Lifting Line Theory 


The concept of a prescribed wake lifting line theory applied to 
propellers consists essentially of assuming that the blade of each propeller 
can be represented by a segmented bound vortex lifting line located along the 
propeller blade quarter chord line with a spanwise varying concentrated circu- 
lation strength which at each segment is proportional to the local blade lift 
(Kutta-Joukowski Law). The wake is assumed to be modeled by a system of 
discrete trailing vortices whose circulation strength is a function of radial 
variations in the blade loading (lift) distribution. A finite length wake is 
used which is of sufficient length to approximately model an infinite length 
wake. Figure 7 is an illustration of this modeling procedure. There are N 
segments modeling the blade which can be arbitrary in length. The circulation 
strength of the bound vortex over each segment is assumed to be constant, and 
the values of the strength are determined by the aerodynamics of the problem 
which are discussed in the next section. Because the flow field sensed by the 
single propeller under steady axial flight conditions is time independent, 
there is no azimuthal variation in blade loading for the single propeller 
configuration. The wake shed by the blade can be modeled by a series of 
trailing segmented vortex filaments shed from the junction points of the blade 
bound vortex segments whose circulation strength is constant over the complete 
length of the filament and equal to the difference in the circulation 
strengths of the adjacent bound vortex segments. Figure 7 also illustrates 
this feature. The trailing filament segmentation is defined by a specified 
azimuthal step size (Ai|> ), each segment modeled by a straight line vortex seg- 
ment. The location of the trailing wake segment end points are prescribed by 
using various types of wake models. 

The influence of the bound and trailing vortex elements at any field 
point is computed by using the Biot-Savart Law for finite length, straight 
line segments of constant strength (reference 8) as shown in Appendix B. In 
this analysis the calculation is done in the cylindrical coordinate system and 
then transformed to the blade element coordinate system. The induced velocity 
vector, Vfg, at a local blade element in the local blade element 
coordinate system due to any trailing segment can be written as 


7 tr s v tr s + V ™c 'c * w T«n 'n 


( 1 ) 


where the subscripts s, c and n denote the spanwise, chordwise and normalwise 
directions respectively. Any one of these induced velocity components can be 
expressed as the sum of the products of the trailing segment circulation 
strength and the appropriate component of the geometric influence coefficient 
computed from the Biot-Savart Law as noted in Appendix B. For example the 
normalwise induced velocity from the 1th segment of the kth filament of the 
mth blade is 
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( 2 ) 


v TRn sr T R (K,m)GC TRn (i,k,m) 


Because the trailing filament circulation for a single propeller configuration 
is the same for all blades at a fixed radial position and constant along the 
trailing filament, the induced velocity component at a field point can be 
expressed as a summation over all of the blades and filament segments of the 
kth filament. 


, t«„ = r T.C>>-2Xec T , n (k,i, m ) (3) 

Letting 

6 C TR n M ^ | GC TRn (k , £,m) (4 j 


the normalwise induced velocity component at a blade element point due to the 
kth filament of all blades then can be written as 

V TR n * r TRM 6 C TR n M (5) 


For the discrete modeling process, the trailing circulation is the difference 
between adjacent bound circulations. Making this substitution in the above 
equation, summing the influence for all filaments and regrouping the influence 
coefficients for like bound circulations into new influence coefficients 
yields the total normalwise induced velocity component for the complete 
trailing wake geometry at a blade element point in terms of the bound circula- 
tion from N segments and the respective wake influence coefficients 

v ™ n *£ r ti) 6 ' < W il (6) 

Computing the bound vortex influences at a blade element point from all 
of the blades yields a similar expression for the normalwise induced velocity 
component in terms of the bound influence coefficients, GCg n and bound 
circulations. 



r (i)Gc Bn (i) 


(7) 


12 



Combining the above two relationships to define the total induced velocity 
component at a blade element point yields for the normalwise component. 


V 


n 


v a + 


'TR n 



r 


(i)GC n (i) 


( 8 ) 


Similar expressions exist for the chordwise and spanwise components of 
induced velocity. Of course, if the blade element point is near a filament 
segment (bound or trailing wake) the influence for the particular segment 
must be handled correctly due to the strong singularity. Generally a core 
radius is assigned to the segment and if the blade element point is within 
that radius the influence is neglected or modified in a prescribed manner as 
described in a later section (Vortex Core Modeling). 

Thus, for a prescribed wake geometry and lifting line geometry, the 
induced influence can be defined in terms of the unknown bound circulation 
distribution. Blade element theory is then used to define the relationships 
necessary to solve for the bound circulation. 
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Blade Element Aerodynamics 


The modeling of the propeller blade by the lifting line approach defines 
the inflow and the effective angle of attack at each blade segment in the 
aerodynamic model. Tabulated airfoil data is used to relate the effective 
angle of attack at each blade segment to the local lift and drag character- 
istics. Also, the use of tabulated airfoil data, acquired from two-dimen- 
sional airfoil tests, inherently accounts for the chordwise vorticity distri- 
bution (related to chordwise pressure distribution) and the Kutta condition. 
This permits the blade segments to be represented by bound vortex segments 
which have a circulation strength representative of the integrated chordwise 
vorticity distribution. For a lifting line model, the unknown quantity is 
then the spanwise bound circulation distribution which is related to the lift 
distribution through the Kutta-Joukowski relationship as will be described. 

In figure 8 the relationships between the relative velocity vectors and 
the geometric quantities in the local blade element coordinate system are 
illustrated. It is assumed that the propeller is operating in steady axial 
flight. This blade element coordinate system is the computational coordinate 
system and the solution for the blade forces is done in this coordinate 
system. One reason to use this coordinate system is that the angle of attack 
is small and the chordwise velocity is the dominant velocity term in this 
coordinate system. This makes the solution numerically easier to obtain in 
the blade element coordinate system. Another reason to use this coordinate 
system is that the resolution of the blade geometry quantities and velocities 
into this coordinate system implicity handles the concept of yaw or skewed 
flow aerodynamics. The concept of skewed flow aerodynamics (in the two- 
dimensional sense) basically states that the pressure force on an airfoil 
section is independent of the spanwise flow velocity component, and thus a 
function of only the velocity components in the plane normal to the local 
chordwise direction (reference for local blade element coordinate system). In 
the relationships to follow, the aerodynamic quantities used are assumed to be 
those which correspond to the local blade element coordinate system. Each 
blade element section is treated as a two-dimensional section with the 
influence of the other sections transmitted through the induced flowfield. 

Linear Aerodynamics 

The aerodynamic relationships for the blade forces include nonlinear 
behavior with certain quantities. The nonlinear solution is obtained by using 
the linearized solution with the appropriate nonlinear iteration techniques as 
will be described in the next section. 

In the blade element coordinate system the local velocity vector diagrams 
appear similar to the equivalent velocity vector diagrams of a statically 
thrusting propeller. In this coordinate system the solution on the blades is 
independent of time. For this reason, the linearization procedure and 
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solution techniques are best performed in the local blade element coordinate 
system. The linearizations follow directly from those presented in reference 
17 for the statically thrusting propeller. 

Consider a local blade element section, at station i, of chord length 
c^ and thickness t^, whose zero angle of attack line is at an angle 
®Bi the local noninduced velocity vector in the plane defined by the 

normalwise and chordwise unit normals. The local lift force, per unit span, 
L^, can be related to the circulation strength of the bound vortex, at 
the quarter chordline by the Kutta-Joukowski relationship: 

Li*/>V ( r, (9) 

where is the magnitude of the total velocity vector normal to the blade 
element section (including induced terms). 

This relationship can be rewritten in terms of the local section lift 
coefficient Cfc£, as: 


r. S -L C V C*. 
*i 2 ' 1 


( 10 ) 


Assuming operation below stall, and linearizing the local section lift 
coefficient as the combination of linear section lift curve slope a^, and 
effective section angle of attack a^, the relationship for the circulation 
becomes: 


r i * 


iVi v i a i 


( 11 ) 


This lift curve slope, obtained from tabulated airfoil data, is of course a 
function of the airfoil type and Mach number for the local section. It is 
noted that the linear lift curve slope assumption and other linearizations 
introduced below are removed by the nonlinear techniques to be described in 
the next section. 

The magnitude of the local velocity vector* in the normalwise chordwise 
plane is just, 


In this section and the following section, upper case V's with subscripts 
denoting component direction always denote noninduced velocity, while lower 
case v's with subscripts denoting component direction always denote induced 
velocity. is the magnitude of the total velocity (noninduced and 
induced) . 
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( 12 ) 


Assuming that in the blade element plane the induced and normal velocities are 
small compared with the chordwise noninduced velocity yields the following 
approximation for the section circulation. 


- _L 
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(13) 


From figure 7 it can be seen that the blade element 
just: 



section angle of attack is 


(14) 


If the chordwise induced velocity is neglected, the section angle of attack is 
assumed to be small, and the normal velocity is also small with respect to the 
chordwise value, the blade circulation is further approximated as 


where 


r = J- 
■ 2 


C i a i V c, 


e B * 




(15) 


(16) 


The normalwise induced velocity, v n £ is a function of the combined wake 
geometry from all blades and the blade circulations, where the normalwise 
geometric influence coefficients are computed using the Biot-Savart Law (see 
section entitled: Lifting Line Theory). Then from equation (8), the blade 

circulation of section i can be written as: 


r i • 


>i c i( 
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(17) 


This relationship is valid under the above assumptions for each blade 
element section, thus a system of simultaneous linear equations in terms of 
the unknown blade circulations for each section can be written as a matrix 
equation in the form: 


MW-W 


(18) 
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and the solution for the circulations can be obtained directly using standard 
solution techniques, either iterative or direct techniques. Because the 
number of unknowns is relatively small, the analysis uses a Gauss-Jordan 
reduction technique. 

Once the blade circulations are known, the section angle of attack is 
known. Using this angle of attack, the actual sectional lift and drag 
coefficients are obtained from tabulated airfoil data (C^ and C^) for each 
section. The components for forces in the blade element coordinate system are 
then computed from this information. Transformation of these forces to the 
cylindrical coordinate system and with the appropriate integrations in the 
axial and rotational directions results in the thrust and torque respectively 
for the propeller blades. The thrust and power coefficients are computed as 


T ' pn V 

(19) 

p 

p ^>n3D3 

(20) 


Nonlinear Aerodynamics 

The assumptions used to obtain a linear solution for the blade aero- 
dynamics can be invalid for several reasons; the linearization of the airfoil 
data may be inappropriate for high Mach number or stalled conditions, the 
small angle assumptions may not be appropriate, and neglecting the induced 
velocities to determine certain quantities may introduce further error. The 
cumulative effect of all of the assumptions may introduce significant 
differences in the solution for certain conditions. To remove all of the 
assumptions used in the linearized aerodynamic model, a technique based upon 
the linearized aerodynamic solution is used. This nonlinear solution 
technique is identical to that presented in reference 17. For this iterative 
technique, the nonlinear aerodynamic relationships are expressed in terms of 
both the linearized relationships and the difference between the linear and 
nonlinear terms, where the difference is obtained from a previous solution 
iteration level. The formulation follows. 


The assumption that the local section lift coefficient is linearized is 
removed by rewriting the section lift coefficient as a linear portion and a 
correction to the linear portion: 

( 21 ) 


where 
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( 22 ) 
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The superscripts denote iteration level. Here, the local lift coefficient is 
obtained from the tabulated airfoil data as a function of the blade element 
section airfoil type, local section blade element Mach number computed 
including induced velocities, and the nonlinear angle of attack. The 
assumption that the total velocity can be approximated by the chordwise value 
is removed by expressing the total velocity as 


where 





( 23 ) 
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( 24 ) 


represents the correction term to the original linear assumption. 

The small angle assumption is removed by computing the section angle of 
attack as 

<>'* 8 , ; + (25) 


where 


and 




( 26 ) 


c/~'= tarT 1 
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( 27 ) 


is the nonlinear correction term. Combining the relationships in the original 
equation for the blade section circulation and expressing the induced velocity 
in terms of the geometric influence coefficients yields: 



is the combined nonlinear correction term. 
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If the combined nonlinear correction term, Cp^, is set to zero, the 
linearized equation for the blade circulation results. If the above equation 
is rewritten in matrix form in terms of the unknown blade circulations, a 
system of simultaneous linear equations results where the correction term is 
lagged in iteration level. 

M {r 1 } = {c +C/" 1 } (30) 


This system of equations can be solved by direct or indirect methods for 
each iteration level where the previous solution is used to obtain the latest 
nonlinear correction terms for each iteration. The initial correction terms 
are obtained from the initital linearized solution. In the current analysis 
the Gauss-Jordan reduction technique is used to obtain the solution at each 
iteration level. 

The resulting local blade forces in the chordwise and normalwise direc- 
tions are obtained in the same manner as is done for the linearized solution, 
and the resultant integrated forces are obtained. 

Skewed Flow Drag Model 

As noted in an earlier section (Blade Element Aerodynamics), the calcu- 
lation of the blade forces in the blade element coordinate system removes the 
necessity of considering skewed flow in an explicit manner in the aerodynamic 
relationships for the pressure lift and pressure drag forces. However, the 
drag coefficient data available for airfoils generally does not distinguish 
between the pressure drag and friction drag components. Since the friction 
drag should be dependent on the total velocity, a skewed flow drag model is 
formulated below using the tabulated drag coefficient data. This model 
attempts to include the additional component of skin friction drag neglected 
by the blade element aerodynamic model described in the previous section. 

Figure 9 is a vector diagram of the drag force and components of the drag 
force at a local blade element station for a condition with a velocity compo- 
nent in the spanwise direction. From this diagram it can be seen that the 
total drag force, D (per unit area), has been broken into two portions, a 
pressure drag force, Dp r , and a skin friction drag force, Df r . The pressure 
drag acts only in the direction of the flow in the plane defined by the chord- 
wise and normalwise directions, while the friction drag is assumed to act in 
the direction of the total velocity. The conventional method of neglecting 
the spanwise velocity component, and computing the drag force, with the 
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pressure and friction forces combined in the drag coefficient (C d ), yield 
the following expression for the components of drag per unit area in the blade 
element coordinate system. 


D S s ° (31) 

Dc *-jP C„V c (V c 2 +V„ Z ) I/J (32a) 

D "*ypCt|V I ,(Vc 2 -fV„ 8 )'* (33.) 


Separating the drag coefficient into pressure and friction portions and 
expanding yields 


°C s i^ v c[ C d pr (v C 2 + V n 2 ) , ic dfr (V c 2 +V n 2 )' /2 ] (32b) 

D n = X^ Vn [ C d pr (v C 2 + v n 2 ) +% fr (V c 2 + V n 2 )' /2 ] (33b) 


for the conventional formulation. 

If the friction drag and pressure drag components are known separately 
it can be seen from the diagram in figure 9 that the components of the drag 
force should be 

D S = y/)V s q lfr (Vc 2 + Vs 2 + V n 2 ) ,/2 

D c = J P * c [c dpr (V c 2 + v n 2 )' /2 + C dfr (v 2 + v s 2 + V n 2 ) ,/2 ] 

D n = y ^Vn[q jpr (Vc 2 + V n 2 ) l/2 + C dfr (V c 2 + V s 2 + V n 2 ) ,/2 ] 


(34) 

(35) 

(36) 


In this section, V c , V n and V g denote velocity components which 
include noninduced and induced velocities. 
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From these relationships it can be seen that for large spanwise 
velocities, the drag force components can be significantly different as 
compared with the conventional formulation. The above formulation has been 
incorporated into the analysis as an optional feature by assuming that the 
skin friction drag coefficient is approximately represented by the drag 
coefficient of the section at zero angle of attack, C<j 0 , and thus the 
pressure drag coefficient is just the difference between the total drag 
coefficient and C £ j c . 


c dpr s *t) C do 


(37) 


Tip Relief Models 


The high speed propeller tip experiences high subsonic and transonic flow 
conditions at normal design operating conditions. These flow conditions, 
coupled with the high degree of three-dimensionality of the problem and the 
use of two-dimensional airfoil characteristics, make the application of the 
lifting line model inappropriate without some type of tip relief scaling 
procedures applied to the airfoil characteristics in the tip region. Two such 
models have been incorporated into the propeller analysis. The first is a 
model based on unswept fixed wing theory. The second is a model based on 
conical flow theory applied to swept wings. 

Ivvard^T^ReUef Model 

For flow conditions where the tip Mach number is greater than one and the 
forward speed Mach number less than one, it is necessary to correct the pres- 
sure forces in the propeller tip Mach cone region. Figure 10 illustrates the 
region on the blade under consideration. In reference 41, charts for a 
scaling factor which is a function of the ratio of blade chord to blade 
radius, radial location on the blade, and tip Mach number are presented which 
can be used to correct the spanwise variation of lift and drag due to lift. 
These charts where derived from the results for fixed wings presented by 
Eward in reference 42. In the current analysis this scaling factor is 
included as an option to the aerodynamic model. It should be noted that this 
model was developed from fixed wing results for application to conventional 
propeller blades without sweep, thus the application to swept blades is ques- 
tionable. 

The application of the scaling function to the lift and drag due to lift 
is done by determining the spanwise location where the tip Mach cone inter- 
sects the blade at a specified fraction of the chord from the leading edge 
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(generally the trailing edge is used). All blade stations from this location 
to the blade tip are modified by the appropriate scaling value if the local 
free stream Mach number is greater than or equal to one. These scaling values 
are tabulated and are an integral part of the analysis. Both the spatial 
location used to define the point for the definition of the tip Mach cone 
location and the fraction of the chord for the intersection location can be 
varied in the analysis. An analytical description for the tip relief scaling 
values is also available in the analysis and it gives slightly different 
values than the original tables. 

£o_nica_l £l£wJTheor£ Ti£ R_e_lie^_M£d_el 

Generally, high speed propeller designs require significant blade sweep, 
therefore the application of the tip relief model based on Evvard's model may 
not be appropriate for these applications. Conical flow theory for thin swept 
wings with sweep can be used to provide a tip relief model for application to 
the swept propeller blade. The model as developed for application to swept 
propeller blades in based on a tapered, swept, fixed wing configuration which 
assumes that the wing leading edge is linear and subsonic and that the wing 
trailing edge is linear and supersonic for a constant supersonic freestream 
Mach number distribution. Furthermore, the span of the wing is assumed to be 
of sufficient length to avoid overlap of the tip Mach cones on the wing 
surface and the wing tip is squared off. These assumptions result in the 
existence of two regions on the surface of the wing for which a scaling func- 
tion based on the ratio of the three-dimensional section lift to the two- 
dimensional section lift can be obtained. The regions on the surface of the 
wing are depicted in figure 11 which also defines the wing planform shape and 
the coordinate system used in the formulation. On the surface of the wing in 
region I along any section chordline of length c, the chordwise pressure 
distribution based on conical flow theory (reference 43) for a subsonic 
leading edge in the above configuration can be integrated along the chordline 
to obtain the appropriate section lift in the region, 
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where 8 ■ (M 2 -!)^ 2 , m ■ 8 cot A^g, k ■ (1-m 2 )*/ 2 


and E* (k) 


is the complete Elliptic Integral of the second kind of modulus k. 


22 



On the surface of the wing in region II the chordwise integration of the 
pressure distribution, approximately a constant deficit along the chordline 
from the value at the edge of this region, is performed to obtain the 
appropriate section lift in the region. 
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These two section lifts are combined to obtain the total section lift. The 
two-dimensional section lift for the same freestream Mach number is known. 



(40) 


The scaling function is then just the ratio of the three-dimensional to the 
two-dimensional section lift. 
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This scaling function can then be applied to the supersonic region of the 
propeller tip by assuming that the length of this region is equivalent to the 
semi-span of a supersonic wing with a freestream Mach number equal to the 
propeller tip Mach number. This scaling function is included as an option in 
the analysis. Since the model as derived is based on a constant Mach number 
and the propeller blade senses a radially varying distribution, an approxima- 
tion for this effect can be incorporated by scaling the relationship of 
equation (40) by the ratio defined below based on Prandtl-Glauert compressible 
scaling rule. 
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(42) 


This feature is also included in the analysis as an option available to the 
user. 
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Airfoil Data 


The use of tabulated two-dimensional airfoil data in the lifting line 
analysis allows for the real effects of compressibility and airfoil type to be 
modeled. The two-dimensional effects for transonic flow operation are thus 
implicity incorporated in the analysis through the use of transonic airfoil 
data. This, of course, means that to obtain reasonable results within the 
limitations of the assumptions of the analysis, the airfoil data must be as 
accurate as possible. The current analysis has available two different sets 
of isolated airfoil data for the NACA 16 series airfoil. This type of airfoil 
is currently used for high speed propeller designs. The details of these 
isolated airfoil data sets as used in the analysis are explained in the 
following sections. 

For the inboard regions of the high speed propeller designs, the use of 
cascade data is more realistic than the use of isolated data. Because the 
inboard sections are generally thicker airfoil types, the NACA 65 series 
airfoil data for cascades has been included in the analysis. Again, two sets 
of data are available. This data is also explained in more detail in a 
following section. 

Hamilton Standard NACA Series 16 Airfoil Data (Manoni) 


Lift, drag and pitching moment data were obtained for a Series of NACA 
16-series airfoils having a thickness ratio range from 0.02 to 0.10 and a 
design lift coefficient from 0 to 0.6 over a Mach number range of 0.4 to 1.25 
and an angle of attack range of -4° to +8°. This data was obtained in the 
UTRC main wind tunnel for the Hamilton Standard Division (HSD) of United 
Technologies Corporation as reported in reference 46. It is generally called 
the "Manoni" airfoil data after the author of this reference. A generaliza- 
tion of lift and drag coefficients for the NACA Series 16 airfoil family span- 
ning a complete Mach number range, based on the data of the reference supple- 
mented by other existing data, was made by HSD personnel for application in 
propeller performance analyses. The transonic similarity parameter approach 
was used because it enables one to generalize data for a limited number of 
sections. Then, the generalization can be used to define aerodynamic charac- 
teristics for a large number of airfoil sections over a complete Mach number 
range. In the transonic similarity form, the aerodynamic equations of motion 
are treated as in dimensional analyses to obtain parameters controlling the 
forces acting on a family of airfoils. This results in a group of equations 
of which the following are of primary interest: 


?L a = C Ua (t/c) l/3 = f,(C,h/t) 

(43) 

c 0p =c 0p (,/c) ■ 

(44) 


24 



^Dc l 2 s C Dc l 2 ^ 1/3 = tyCiM) 


(45) 


C = (M 2 -|)(t/c )" 2/3 (46) 

The above symbols are defined in the List of Symbols. 

In order to make the generalization applicable to higher thickness 
ratios, additional adjustments had to be included. This resulting generaliza- 
tion spans the following ranges: 

Mach Number 0.3 to 1.25 

Angle of Attack -4® to +8® 

Thickness Ratio 0.02 to 0.40 

Design Lift Coefficient 0 to 0.7 

and has been implemented into this propeller-nacelle performance analysis. 

Published NACA Series 16 Airfoil Data (NACA) 


In reference 47 two-dimensional airfoil data for NACA Series 16 airfoil 
sections are presented. The data was developed from a compilation of all the 
data available up to the early 1950's and the effort was somewhat handicapped 
by the lack of data throughout the transonic Mach number region. The data 
cover a range of thickness ratios of 4 to 21 percent and design lift coeffi- 
cient of 0 to 0.7. The lift data are presented in the form of curves of lift 
coefficient, C^, as a function of freestream Mach number and airfoil design 
lift coefficient. Data from the reference was partially computerized by HSD 
personnel spanning the following ranges: 


Mach Number 
Angle of Attack 
Thickness Ratio 
Design Lift Coefficient 


0.3 to 1.15 
-6® to +12® 
0.03 to 0.09 
0 to 0.2 


The data was extrapolated to obtain 0.03 thickness ratio. The thickness ratio 
and design lift coefficient ranges can be extended from the data provided in 
the reference. This data is generally labeled the "NACA" airfoil data. 

Cascade Airfoil Data (NASA SP-36) 


Since the inboard sections of a high speed propeller design can approach 
gap to chord ratios which represent cascade airfoil conditions, a set of 
cascade airfoil data correlations has been incorporated in the analysis for 
typical high speed propeller design requirements. The cascade correlation is 
based on the NACA 65 series A = 1 meanline. This cascade correlation is valid 
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for gap to chord ratios greater than 0.8 for a Mach number range up to 1.0. 
The correlation is valid for ±10.0 degrees in angle of attack and the thick- 
ness to chord ratio available up to 0.1, with a numerical design lift coeffi- 
cient range of -3.75 to 3.75 (camber angles of ±90 degrees). The description 
of the theory used to relate the cascade correlation parameters to isolated 
airfoil parameters so that it can be used in the lifting line solution 
follows . 

T^h_e 0 £y_-_C_a s^a d^e_C o^r£ e ion 

The theory used for the cascade correlations used in this subroutine is 
given in NASA SP-36 (reference 40). It is based on the observation that the 
exit flow angle from a blade is close to the metal angle, or trailing edge 
angle, a 2 * (see figure 12). Thus, only the deviation angle 6 Q is correlated 
with cascade parameters. This basic correlation has the form 


a 2 = a s - 



where the deviation angle 6 q has the form 

^o = K Ss ' * &oo + M a <£ c /cr b 


(47) 


(48) 


The parameters appearing in the above equation were obtained by curve fitting 
to the charts in reference 40. 

The minimum loss incidence angle i ^ has the form 

mo 

•m 0 = K it' k 8$ *oo + n< £c (49) 


where again the parameters were obtained by curve fitting the charts in 
reference 40. 

The minimum loss is based on the correlation of reference 40, 

ZBm = 2C $(D)cOS* 0| /C0S 3 Ct2 (5°) 

where $(D), a correlation parameter, is obtained by curve fitting and the 
diffusion parameter D is given by 
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cos a, i cosa, 
0 = cos a 2 + T cr “ 


tan Q i 



(51) 


Equations (47) through (51) are the basic relations used in reference 40. 

These relations have been modified and expanded as follows. The first is 
a back pressure correction given by 


* i tan q g 

° 2C (pUj)2/(/oUj)| 


(52) 


The second is a calculation of loss at any incidence angle given by 


2 B 


^ [»m] 


(53) 


where i is the incidence angle for twice the loss which is determined by 
curve fitting. These empirical relations yield predictions in good agreement 
with Pratt and Whitney Aircraft cascade predictions. 

Theory - C L > C D Correlations 

The Cl, C d data correlations are based on thin airfoil theory and the 
data given in reference 48 and in particular the series 6 airfoils which were 
developed for low drag at the design lift coefficient Cj^. The 
theoretical lift curve slope is given by 


C Lqo = 27T (54) 

Since the data shows that the actual lift curve slope is less, the theoretical 
lift curve slope is modified by an empirical factor 

c La 0 = 27r(l-€) (55) 
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Thus, for thin airfoils, the lift is given by 


C L = 


C Lo +C i-ao a 


( 56 ) 


as shown in figure 13. 

For thin airfoils Cl 0 is primarily a function of camber angle 
It can be shown that the design lift coefficient is related to the equivalent 
circular arc camber by the relation 


C i d = <k(l-0.05<£ c 2 ) 


(57) 


Thus, Cl q is determined empirically by the ratio of design lift coeffi- 
cients, 
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and Cl q is determined from data for a reference airfoil. At large angles 
of attack, the Cl departs from the linear relation. Hence, experimental 
data was fitted to an equation of the form 


C L =C Lo +C Lao a +A(a-a 0 ) n 


(59) 
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Finally, for Mach numbers less than the critical Mach number M cr , the 
Prandtl Glauert rule is used. 


C Lt 

C L = / a- M<M Cr ( 62 ) 
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The drag coefficient in the drag bucket is fitted to an equation of the form 


C D ' C 0 Bo 


1 + AC 
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/ Q Q min \ 
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and outside the drag bucket to an equation of the form 
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(64) 


where the respective quantities are shown in figure 13. 


The conversion of lift and drag coefficients (Cl, Cp) to an exit 
angle (a 2 ) and loss coefficient (Zg) is based on the resolution of the 
forces whose directions are shown in figure 12. The axial and tangential 
forces T2 and T^, for a two-dimensional cascade using compressor termi- 
nology and the definition of Cl, Cp are given by 


T 2 = n b y />, U e 2 c ( C L sinae - C 0 cos a e ) 


(65) 
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where the effects of skewed flow are reflected in Cl and Cp predictions 
based on the local blade element coordinate system. Dividing these equations 
through by (1/2 pU e ^) Ngc we have the force coefficients, 


C 2 * Cl sin a«“ C D cos ae 
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(69) 


(70) 
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Rewriting equations (69) and (72) in terms of cascade parameters we have 
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C<£ = 
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where 


The static pressure ratio across the blade row is given by. 


/5l\ s fS. f2L s / l+(y-l/2)M? \ y/yH 

\ P| ) Pqi ^02 P l Po * / 


and M 2 is given 


yM ' /Z 

^ B (l + (y-l/2)M?) y/y ‘ X 


M 2 . ( C -^— I ) U1 


(76) 


(77) 


(78) 


30 



t 


where 


(79) 


((l+( r -|/2)Mf) l/Z ) 
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The axial and tangential momentum induced velocity can be related to the 
thrust and torque. 


W* = T* /{ N^g [(pUj)2 * (/DU* ijjj’ (80) 

v{N B g[(pU 2 ) 2 +(^U z ),]} (81) 


hence 



(82) 


(83) 


For many blades the momentum or gap average induced velocities approach the 
local induced velocity. Finally, we note that the resultant angle of attack 
can be approximated as 


Qi 



/ Sing, -w 2 /U| \ 

\cosa|-w^/U| ) 


(84) 


which can be used to obtain the effective Cl and Cj) from the cascade 
correlations. 

jL°JL u JLi£ n _j Algorithm 

For the cascade correlation the exit air angle a 2 and loss coefficient 
Zg are known. The lift and drag coefficients and induced velocities can be 
determined directly as follows. The total pressure ratio i® given by 
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equation (77) and the downstream Mach number M 2 by equation 78. Then the 
compressor force coefficients when corrected to account for the absence of the 
induced velocity, 



(85) 


are determined by equations (73) and (74). Then the induced velocities are 
determined by equations (82) and (83). The resultant velocity is determined 
by equation (75) and the induced angle by equation (84). Then the resultant 
angle is given by 


Qg ~ Qj CXj 


( 86 ) 


and equations (69) and (71) solved simultaneously for and with C^, 

C^, and a g known. 

NASA_6^A_=_l__S£rJ L e£ Ai£foi_ls 

The NASA equivalent circular arc camber angle for the A = 1 meanline is 
arbitrarily defined so as to produce good data correlations. The leading edge 
and trailing edge metal angles are defined relative to chord angle a g by 

a* = tarf nnn * = tan' 1 (0.4212 C ld ) 

I Vdx/xs 0.005 a ( 87 ) 


a? s tan" 1 (•—) =tan _l (-0.2343 Q.) (88) 

2 \dx/x*o.95 


and the camber angle by 


( p c = 4 ton'H 0.1103 Ci d ) 


(89) 
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Hence, we have 


a* * a s + tan '(0.4212 0^) 
aj = a s + tan" '(-0.2343 C^) 


(90) 


with aj* and c^* given, equation (90) may be solved for a g and C]^. 

Thus the camber is given by equation (89). 

£asc*id£ Ai£foi_l_ Data_(Me HorO 

An alternate source of cascade data is provided in the analysis based on 
the theoretical and empirical correlations of Mellor (references 49 to 51) and 
the data of reference 52. This cascade correlation is for the NACA 65 series 
airfoil sections and also includes several other airfoil types. In the 
current analysis the module which contains the correlations has been 
programmed to provide only the 65 series correlation, but this can be easily 
modified at some future date. The correlations as developed by Mellor are for 
incompressible flow, and as such should be used with the recognition of that 
limitation in the application to the high speed propeller. The correlations 
in the analysis are functions of camber, gap-to-chord ratio, thickness ratio, 
blade stagger angle, Reynolds number, and angle of attack. 

Cascade Correction for Isolated Airfoil Data (Flat Plate Theory) 

Included in the analysis is an analytic cascade correction which may be 
applied to isolated airfoil data for which actual cascade data does not exist. 
The correction procedure as used in the analysis comes from reference 44 and 
is based on analytical results presented in reference 45 for two dimensional 
flat plate airfoils and can be used for cambered airfoils if the local geo- 
metric blade angle 0, and the angle of attack, a are measured to the zero lift 
direction. 

This correction consists of modification of the isolated airfoil angle of 
zero lift due to the influence of camber by a difference Aa 0 , and thus 
modifies the angle of attack for the table lookup of the lift coefficient. 

The resulting lift coefficient is then scaled by the cascade correction factor 
K, which is a function of the local gap to chord ratio, g/c and local geo- 
metric blade angle 0. 
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Cj i = KC £ (a-Aa 0 ) 


(91) 


The values of K used in the analysis are presented in figure 14. 'Because of 
the approximate nature of this model, it is preferable to use actual cascade 
data if possible. 

Cascade Correction for Isolated Airfoil Data (Mellor) 


As a consequence of providing a cascade correlation using the correlation 
of Mellor (reference 49), the same correlation can be used to provide a 
cascade correction for isolated airfoils. The correction is an optional 
feature in the analysis which can be used in place of the correction based on 
flat plate theory. The correction is obtained by calculating the ratio of 
section lift coefficients of the NACA 65 series at the desired gap-to-chord 
ratio versus an infinite grfp-to-chord ratio for the isolated airfoil operating 
conditions. This ratio is then used to scale the isolated airfoil section 
lift coefficient obtained from the appropriate data base of isolated airfoil 
characteristics at the operating conditions. 

Cajc_acie and^ i»£l£ted_A^r^oi l_Dat£ I^nt_e_rp£l_at_io_n 

An interpolation procedure is available in the analysis to provide 
airfoil characteristics along the blade in the transition region for which 
cascade data does not exist and isolated airfoil data does not rigorously 
apply. The interpolation procedure in this region is based on the gap-to- 
chord ratio (t). The airfoil characteristics are obtained by scaling the 
isolated airfoil characteristics for the particular blade element section by 
the use of an interpolation function. 


C A " C l T cn f i (T> 
ISO 

C d * C d M T ) 
® a ISO a 
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This function is based on the concept that the airfoil characteristics in the 
transition region must approach the isolated airfoil characteristics from the 
cascade values in a continuous and asymptotic manner. A function which has 
this behavior is 


f(T) 


a + T n 
1 + T n 


(93) 
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where the exponent (n) on the gap-to-chord ratio is defined by the user. The 
value of the exponent controls the rate at which the function asymptotically 
approaches unity for large gap-to-chord ratios. Larger values of the exponent 
increase the rate. The value of the parameter (a) is determined in the 
analysis by the ratio of the cascade to isolated airfoil characteristics 
(R c ), where the cascade data is obtained from the existing correlations for 
a gap-to-chord ratio at the boundary of the correlation data (og) . 


a - R c + a B n (R c -1) (94) 


This modeling option can be used in place of the cascade correction for 
isolated airfoils described in the previous sections. 
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Wake Modeling 


The use of a lifting line model in this analysis requires that the wake 
geometry be prescribed for computational purposes. As noted earlier the wake 
is represented by a system of trailing vortex filaments. Within this analysis 
it is possible to prescribe internally several different wake geometries or 
input the wake geometry from an external source. The reasons for this 
versatility are twofold. First, at the conventional operating conditions for 
the high speed propeller, the true shape of the wake is not known, thus the 
need for the ability to include a versatile range of wake geometries. It is 
however clear that the wake shape at these speeds is not significantly 
distorted. Secondly, at low speed or static flight conditions the wake shape 
is highly distorted as shown in reference 26, thus the need for wake models 
which can vary significantly in shape is well documented. The use of the 
appropriate selection of the wake geometry model for the flight condition 
under investigation is important for accurate propeller induced velocity 
predictions. Preferred wake models for the various flight speed regimes will 
be indicated. However, the final selection of wake model should be based on 
the results of a wake sensitivity study which remains to be conducted in 
future applications of the analysis. The wake models which can be used in the 
analysis are described in the following sections. 

Classical Wake Model 


Early lifting line models, for either static or forward flight condi- 
tions, generally used what is called a classical wake model. This model 
consists of defining the axial wake transport velocity as the addition of the 
forward speed and the momentum induced velocity for the thrust and flight 
condition being investigated. 


Vt 


= V 0 + v. 

0 i mom 


(95) 


No radial wake contraction is used. The resulting wake shape is an uncon- 
tracted helix for which the pitch rate depends on the flight condition and 
thrust level. Figure 15 illustrates the model. This model can be used in the 
analysis for any flight condition. Although it is known that for static 
thrust conditions this type of wake model will give incorrect answers, it is 
in the analysis because the self-induced distortions of the propeller wake 
decrease with forward speed. The classical wake or the classical wake as 
modified by superposition of the nacelle influence (to be discussed) can be 
used for the high speed propeller application. 

Modified Classical Wake Model 


Variations of the classical wake model can be used in the analysis if 
desired to approximately account for both the nacelle influence and/or the 
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self-induced wake distortion. One model allows for a radially varying axial 
wake transport velocity which replaces the uniform value used in the classical 
wake model. 


v, = f„(x> ■ V 0 +v imom 


(96) 


where the function f y (x) represents a spanwise velocity scaling distribution 
which varies radially along the blade. As in the classical wake model there 
is no radial wake contraction. The function f y (x) is normally obtained from 
knowledge of the nacelle influence at the blades. If desired, it can also 
account for the difference between the propeller self-induced velocity span- 
wise distribution and the momentum values. Figure 16 illustrates this wake 
model. The use of this model is assumed to be more representative of the wake 
of the high speed propeller nacelle configurations than that generated by a 
classical wake model. However, alternate procedures will be described. 

Generalized Wake Model 


For static thrust conditions, wake distortions are significant in terms 
of propeller performance predictions. Because of this fact, a wake model 
which has been shown to accurately model the statically thrusting propeller or 
hovering helicopter rotor wakes (reference 17 and reference 27) can be used in 
the analysis. This wake model is known as the generalized wake model. It is 
based on empirical results, and has radially varying axial wake transport 
velocities and radial contraction of the wake. Figure 16 illustrates the 
general character of this wake model for a statically thrusting propeller (see 
also figure 4). The wake geometry varies in both the axial and radial direc- 
tions and is highly dependent on wake age, i|> w . The wake age of a point in 
the wake, i|> w , is the azimuth interval of blade rotation between the current 
blade azimuth and the azimuth at the time the wake point was emitted from the 
blade. A detailed explanation of the wake model is available in reference 17 
or reference 26. The equations for the tip vortex coordinates normalized by 
the blade radius are noted below: 
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in terms of the generalized wake parameters. 
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The cross sections of the vortex sheet from each blade are essentially 
lines rather than discrete points, as is the case for the tip vortex cross 
sections. The vortex sheet cross sections are assumed to be linear, the axial 
position of a vortex sheet at a given azimuth can be defined by two points. 

For simplicity, the two points are the imaginary extensions of the cross 
section to r * 0 at one end and to r 8 1.0 at the other end, as shown in 
figure 17. These two points establish the intercept at the axis of rotation 
and the slope of the vortex sheet. It was assumed that the radial position of 
such intercepts (e.g., point A of figure 17) is linearly proportional to the 
radial coordinate of the intersection of the vortex sheet with the vortex 
sheet boundary (point B). With the exception of the immediate vicinity of the 
blade (^ w ( 2ir/b), the vortex sheet boundary was assumed to be equivalent to 
the boundary formed by the locus centers of the tip vortex cross sections 
(equivalently, tip vortex streamline). For i|> w less than the blade spacing, 
2x/b, the boundary was faired from the point of maximum circulation on the 
blade (point 5 of figure 17) to the tip vortex boundary. The constant of 
proportionality was assumed to be the ratio of the radial position of the 
origination of the filament streamline at the blade (point C) to the radial 
position of the vortex sheet boundary at the blade (point D). The equations 
for the vortex sheet coordinates normalized by the blade radius are noted 
below: 
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The generalized wake coefficients (A, X,K ^ ,1^ >K i r =i>^2r=l ,^2r=o) 
define the wake shape must be determined in advance and input to the analysis. 
The values can be determined from experiment or from reference 17 for static 
propellers or reference 26 for hovering rotors. 
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Because the wake geometry must change smoothly from the highly distorted 
shape under static thrust conditions to a basically undistorted shape at high 
forward speeds with a nacelle, the generalized wake geometry noted above has 
the radially varying velocity distribution (neglecting the momentum induced 
velocity) superimposed on it to account for the flight speed and the influence 
of the nacelle. Thus if the wake shape can be determined for low* forward 
speeds from some source of information, it is expected that this model will be 
useful in modeling that flight regime’s wake characteristics with the appro- 
priate input wake coefficients. 

Nacelle Influences on the Wake Geometry 

The presence of the propeller nacelle can alter the local flow at the 
propeller blades and wake, particularly near the nacelle walls. The wake 
model described in the Modified Classical Wake section is based on the assump- 
tion that the nacelle does not influence the radial location of the wake geo- 
metry and that the axial transport velocity is constant for all axial loca- 
tions. In reality this is not the case. In the analysis, it is possible to 
apply a second approximate wake perturbation model which is a function of the 
influence of the nacelle over the complete region exterior to the nacelle 
boundary. This allows the wake to accelerate and decelerate, expand or 
contract to follow the flow around the nacelle. This model is based on the 
superposition of the nacelle's influence (e.g., as determined from the nacelle 
streamline solution) on one of the previously described wake models aa used to 
represent the isolated propeller wake. That is, once a wake model is selected 
for an isolated propeller, the wake can be corrected for the presence of the 
nacelle. The equations for the nacelle influence are presented in a later 
section (Nacelle Wake Corrections). Figure 18 illustrates this concept. It 
is recognized that this model, based on the principle of superposition is only 
an approximation to the actual wake. However, it is believed to be the most 
accurate method available, short of prescribing empirical values from measured 
data, which does not require calculating the complete coupled propeller- 
nacelle wake interaction problem numerically. 

Wake Input Models 

Since it may be necessary to prescribe wake geometries which are not 
compatible with any of the currently available wake models, the analysis will 
accept wake geometry input compatible with the cylindrical coordinate system 
mentioned in an earlier section. Thus, if it is determined that low or high 
speed flight conditions require currently unknown wake geometries for the 
propeller-nacelle configurations under consideration, the analysis will be 
able to accept these geometries. 
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Wake Rollup Modeling 


The representation of the wake of a propeller blade as noted earlier 
consists of modeling the shed sheet of vorticity with discrete vortex fila- 
ments. The actual wake, however, is known to roll up into a tip vortex and 
inboard sheet. The tip vortex has, in general, a strength equal to the 
maximum circulation strength on the blade, as shown in figure 19. The occur- 
rence of a rolled up root vortex has not been clearly demonstrated for 
hovering rotors or statically thrusting propellers, where wake distortions are 
most important. The modeling of this phenomena, if it exists at all, is 
generally not considered. For conditions where wake distortions are 
important, it is generally the distortion of the tip vortex which dominates 
the flow field characteristics. In the analysis, it is possible to prescribe 
the wake rollup for either the inboard or tip region if so desired. The age 
of the wake rollup and the concentration of circulation strength are deter- 
mined from input quantities. Generally the tip vortex rollup is modeled to 
occur within 10 to 15 degrees behind the blade. The use of inboard wake roll- 
up is not justified from a physical point of view as noted above; however, for 
conditions where a significant length of wake is required (static thrust) it 
is possible to reduce the computation time in the analysis by specifying 
inboard rollup to occur in the far field wake without sacrificing computa- 
tional accuracy. 

Vortex Core Modeling 

If in the course of the use of this analysis, the prescribed wake model 
is such that a vortex filament segment is positioned close to a field point 
(blade segment center) where the induced velocity is being calculated, the 
predictions using the Biot-Savart Law for the magnitude of the induced 
velocity may exceed those which actually occur in real life for a discrete 
line vortex such as the tip vortex or where the vortex sheet is modeled by 
discrete vortex filaments. To resolve this problem, the analysis includes 
a vortex core model which is normally used. This model assumes that the 
variation of the velocity magnitude in the core is constant with variation 
in the direction normal to the filament segment and equal to the velocity 
at a specified core radius. Figure 20 illustrates this model for a finite 
length vortex segment. Although it is known that the velocity in the core 
of a discrete vortex decays from the peak value at the core radius approxi- 
mately linearly to zero at the center, this simple model is felt to be as 
valid as any other, on the basis that if a close passage does occur, the 
modeling of the blade and wake by discrete vortex filaments is in itself 
inappropriate. 
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Compressibility Considerations for Induced Velocity 

When the performance of two-dimensional airfoils is found from the air- 
foil data, the effects of compressibility are included, as the data is tabu- 
lated as a function of Mach number. To use the airfoil data for the case 
where the propeller is operating at high forward Mach numbers, compressibility 
effects on the induced velocity must also be determined. 

The conditions for which the effects of compressibility on the induced 
velocity need to be considered are outlined in reference 41. In summary, if 
the local section Mach number (total of the rotational and forward Mach 
numbers) at the blade exceeds the local speed of sound, it is necessary to 
consider the effects of compressibility on the induced velocity. The proce- 
dure to incorporate the effects of compressibility is also outlined in the 
above reference for vortex models. To compute the effect, it is necessary to 
consider the finite speed of the pressure disturbances in comparison with the 
speed of the blade section being analyzed. This leads to the rule of 
forbidden signals when some of the sections are operating at section Mach 
numbers above one. Thus the effects of the line vortex shed from a section 
operating at a section Mach number above one are considered only at points in 
the region of its aft Mach cone, where the apex of the Mach cone is attached 
to the section from which the line vortex is shed on the propeller blade. 
Figure 21 (from reference 41) is a pictorial representation of this concept 
and illustrates the blade and Mach cone which defines the limits of the 
possible induced influence of the trailing line vortex. For the case of a 
high speed propeller operating with blade section Mach numbers greater than 
one but forward Mach numbers less than one, the initial portion of the trail- 
ing line vortices shed from these blade sections will have no influence at the 
propeller blades. However, after some finite period of time the influence of 
a portion of the line vortex will reach the propeller disk. Figure 22 (also 
from reference 41) illustrates this concept, where the position of the blade 
and the Mach cone boundary are plotted for different rotational positions of 
the blade (time). It is shown in this figure that for a forward Mach number 
less than one, the induced influence of a portion of the line vortex does 
eventually reach the blade after some finite time (between time C and D in 
figure 22) . 

To calculate the induced inflow velocity of a propeller operating at a 
tip Mach number exceeding one and a forward Mach number less than one, it is 
necessary to evaluate the effects of the trailing vortices that lie in the 
zone of influence only. Thus, unlike the incompressible flow case where 
effects of the vortices are considered from the propeller plane to infinity 
instantaneously, it becomes necessary to discriminate against the effects of 
the vortices outside the zone of influence or the Mach cone as shown in figure 
21. Thus the portion of the trailing vortex filaments which has no influence 
on a blade segment at a particular time should not be considered when the 
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calculation of the geometric influence coefficients using the Biot-Savart Law 
is made as noted in the section entitled: Lifting Line Theory, This approach 

has been included in the analysis. The concept of the induced influence of a 
trailing line vortex limited by compressibility considerations can also be 
applied to the bound vortex segments representing the blade. The validity of 
applying this concept to the bound vortex influence has not been established; 
however, it is available in the analysis if so desired. 

As noted above, the concept of wake compressibility as outlined in 
reference 41 assumes infinite propagation of the induced field of the trailing 
line vortices shed from blade sections operating at section Mach numbers less 
than unity and finite propagation (speed of sound) of the induced influence 
when the sections are operating at section Mach numbers greater than unity. A 
second feature, which is not normally used, is provided in the analysis for 
exploring the influence of finite rather than infinite propagation speed for 
blade sections with a section Mach number less than one. That is, a provision 
is included which assumes that the influences of all bound and/or trailing 
wake elements propagate at the speed of sound. 
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Coaxial Theory - Equal Blade Number 
and Rotational Speeds 

For the coaxial, counter-rotating propeller analysis, the assumption of 
an axisymmetric steady state flow problem with respect to the blade coordinate 
system must be removed. In the blade coordinate system, the circulation solu- 
tion becomes periodic with half blade spacing for equal blade number and 
rotational speeds. This feature is demonstrated in figure 23 where it can be 
seen that the geometry of the problem is periodic with half blade spacing. 

The half blade spacing is the azimuthal interval of blade rotation between the 
two propellers required for the relative geometry to repeat itself. The use 
of lifting line theory is still valid within the constraints of the other 
assumptions noted in earlier sections. In this analysis the only major 
conceptual difference between the single and coaxial solutions at the blade 
rows is the above noted periodicity. The coaxial analysis assumes quasi- 
steady aerodynamics. Thus the relationships used to determine the unknown 
blade circulations only need to be modified to include the solution period- 
icity. The details of such a deviation are straightforward and only the 
resulting relationships are represented herein. 

Coordinate Systems 

Within the coaxial lifting line analysis, definitions of the coordinate 
systems are identical with those in the single propeller analysis. The 
cartesian and cylindrical coordinate systems are referenced with respect to 
each blade row, with the definitions of positive quantities consistent with 
the respective direction of rotation of each blade row. In the analysis there 
is no distinction between which blade row is rotating in which direction, only 
that they are counter-rotating; this is because the problem is symmetric in 
nature. The definitions of velocities computed at a blade row are consistent 
with that blade row's coordinate system. Figure 24 illustrates the above 
noted coordinate systems. For calculation purposes, the analysis requires 
that the hubs of the two propellers be separated by a finite axial displace- 
ment. Propeller one (in terms of input ordering) is assumed to lie in front 
of propeller two, and a positive hub displacement positions propeller two 
behind propeller one (figure 24). The blade element coordinate system for 
each blade row is consistent with the respective cartesian and cylindrical 
coordinate systems. 

Lifting Line Theory 

As noted earlier, the addition of a counter-rotating propeller requires 
the inclusion of geometric periodicity and hence circulation periodicity in 
the relationships used to compute the solution. For the calculation of the 
induced velocity at a field point using lifting line theory, the relationships 
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derived earlier must be expanded to include the periodic nature of the circu- 
lation solution of each blade row and the wake influence. The resulting 
relationship for a component of induced velocity, in this case the normalwise 
component, at an arbitrary field point is: 
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where it is assumed that there are identical number of blade stations (N) at 
each blade row and that there are L propeller positions between half blade 
spacings for a given blade azimuth step size. The term, propeller position, 
means the discrete azimuthal location of the propeller within a specified 
rotational interval, defined by a finite azimuthal step size. In this case, 
the interval is the half blade spacing, and the discrete azimuthal locations 
are defined by the azimuthal wake segmentation interval. 

Blade Element Aerodynamics 

The effect of including a second propeller blade row on the aerodynamic 
relationships is to expand the number of unknown bound circulation strengths, 
and hence the size of the matrix of equations which must be solved using 
either the linear or nonlinear solution. The relationship for the unknown 
circulation at station i for propeller position j of propeller p for the 
nonlinear solution becomes: 
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where the correction term has the appropriate modifications for the 

coaxial propeller condition. Dropping the superscripts and the correction 
term yields the linear solution. From this relationship for the unknown bound 
circulation at each station on each propeller blade for each propeller posi- 
tion a matrix of NxLx2 simultaneous equations can be written, and the same 
solution methods can be applied as for the single propeller condition. 

The resulting blade forces can be computed and integrated for each pro- 
peller at each propeller position. The resulting values can be averaged over 
the propeller positions to obtain the steady thrust and torque. For the 
nacelle code these steady forces are the values used for the coaxial case 
since the nacelle analysis assumes steady flow conditions. . 
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Other Considerations for Coaxial Propellers 


The use of the different wake models and airfoil data for the coaxial 
propeller configuration is the same for the single propeller configuration. 
There is only one restriction in the analysis which limits the use of these 
portions of the analysis to single propeller configurations. That restriction 
is the application of the rule of forbidden signals to the bound vortex which 
is noted in the section entitled: Compressibility Considerations on the 

Induced Velocity Calculation. The analysis will not allow the application of 
this modeling feature to the coaxial configuration. The reason is that the 
formulation for this model as applied to the bound vortex has not been 
developed for the coaxial configuration because this model is in itself 
questionable. 

Because the wake of the leading propeller can pass close to the lifting 
line representation for the blades on the following propeller, care should be 
taken to inspect the resultant blade loadings to look for large gradients in 
loading due to close blade-vortex interactions. The validity of the solutions 
where these close passages occur is questionable. The use of different core 
radii for the wake vortex systems can influence the answers when these close 
passages occur, and so careful selection of the appropriate values should be 
made for the coaxial configuration. 
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Coaxial Theory - Unequal Blade Number 
and Rotational Speeds 


For the coaxial configuration, the introduction of unequal number of 
blades and unequal rotation speeds on the propeller disks creates, an addi- 
tional complexity in the analysis. The concept of periodicity of the 
propeller blade/wake geometry and periodicity of the circulation solution 
still exists, but the definition of azimuthal interval must be replaced with a 
time interval as the quantity which defines this periodicity. This is due to 
the fact that the propeller blades do not necessarily rotate at the same 
speed. In general, it can be shown that the time interval for which the 
propeller blade and wake geometry are periodic is, 
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where b is the maximum of the number of blades on each of the propeller disks. 
This relationship reduces to the conventional one for coaxial propellers with 
equal number of blades and rotational speeds on each disk; the periodic 
azimuthal interval is one half of the blade azimuth interval. 

The relationships for the formulation of the numerical solution procedure 
are no different than those for the conventional coaxial configurations. 
However, the implementation of the formulation in the computer code must be 
made in terms of time intervals instead of azimuth intervals. For unequal 
rotational speeds on the propellers, this results in different azimuth 
intervals used to define the propeller positions and wake geometries. 

The analysis includes this feature to treat the unequal number of blades 
and rotational speeds with certain limitations. These limitations are related 
to computer program storage requirements and are not conceptual in nature. 
Large variations between the number of blades on each propeller disk, or large 
variations in rotational speeds, result in periodic time intervals which 
require a large number of unknowns in the solution formulation. This increase 
in the number of unknowns increases the physical storage needed in the program 
and the computational time required to obtain the solution. The current 
analysis is limited to combinations of blade number and rotational speed which 
result in periodic time intervals with corresponding azimuthal intervals on 
the respective propeller disks of no less than five (5.0) degrees. This 
limitation can be relaxed by re-dimensioning the computer code if deemed 
necessary at some future time. The analysis controls the selection of the 
actual time increments used in the definition of the wake geometries based on 
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the periodic time interval relationship defined above. The analysis attempts 
to find a time increment which is at least one third of the periodic time 
interval for the particular selection of blade numbers and rotational speeds 
and which does not result in an azimuthal interval for either propeller less 
than the current limit of five degrees. The user can also limit the time 
interval to result in an azimuthal interval on either propeller which is not 
greater than a specified value. 

The analysis has been coded in such a manner that the unequal blade 
number and rotational speed feature does not result in a significant computer 
overhead penalty when applied to the more conventional or coaxial operating 
configurations . 
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TECHNICAL APPROACH - NACELLE 


Overview 

A large and important class of problems is represented by viscous 
turbulent flows over axisymmetric nacelles with the following properties: 

1) moderate to high Reynolds numbers, 2) strong inviscid-viscous interactions, 
3) large regions of inviscid flows, and 4) significant "elliptic" effects due 
to streamline curvature and streamline divergence. Accurate and reliable 
calculation of such flows is of immediate and practical interest in the 
prediction of flows over a nacelle for either external or internal flow. 

This section consists of several parts. First, a procedure is presented 
for constructing an orthogonal coordinate system from the planar potential 
flow solution using the stream function as a normal coordinate and the veloc- 
ity potential as the streamwise coordinate. Since the potential flow stream- 
lines approximate the viscous flow streamlines, the thin channel approxima- 
tion, references 4 and 5, may be applied to the Navier-Stokes equations 
written in this potential flow coordinate system to reduce them to a set of 
parabolic equations which can be solved by forward marching numerical methods. 
Second, a method using conformal mapping based on the Schwartz-Chr istof fel 
transformation is presented for obtaining the plane potential flow solution 
for the flow over an axisymmetric nacelle. This procedure has the important 
property that the inverse solution is obtained directly through the transfor- 
mation so that equipotent ial lines and streamlines can be used to construct 
the coordinates. Third, a solution algorithm based on the two point box 
scheme is presented for solving the viscous flow equations. It should be 
noted here that the term, duct, is used even when considering the external 
flow problem of the propeller-nacelle geometry. For these external flow 
problems, the outer wall of the duct becomes the outer limit for the computa- 
tional region with the appropriate modifications to the boundary conditions, 
and the inner wall becomes the nacelle profile. And finally, because of the 
particular requirements of the nacelle application, streamline curvature 
corrections to the thin channel approximation are necessary (reference 33) and 
are included in the analysis. 
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Analysis 


Streamline Coordinates 


The equations of motion are solved in a streamline orthogonal coordinate 
system (n, s, 4>) . This coordinate system is generated by first solving for 
the plane potential flow (independent of $) through a duct with the same 
cross-section that the annular duct makes with the meridonal plane. The 
normal coordinate n is the stream function and the streamwise coordinate s is 
the velocity potential. Rotation about the axis of symmetry produces an 
axisymmetric orthogonal coordinate system uniquely suited to solve the 
problem. 

The (n, s) coordinates are related to the physical coordinates (r, z) 
through Laplaces' equation. 


d 2 n d 2 n 


<J 2 s d z s 
dT 2 " + dz z 



( 105 ) 


( 106 ) 


The metric scale coefficients are the same in both directions and are equal to 
the inverse of the magnitude of the potential flow velocity, V. 


V* 




( 107 ) 


Lengths along the streamlines and potential lines are given by 


dx =ds/v 


( 108 ) 


dy = dn/V 


( 109 ) 
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I 


and curvature of streamlines and potential lines are given by 


1 _ 

av 

(110) 

R s 

dn 

1 

j3V_ 


R n " 

<Js 

(111) 


In this coordinate system (n, s) are independent variables and (r, z) 
dependent variables. While one could solve this problem numerically, an 
alternate, more efficient procedure is available. Solution of the inverse 
problem (r(n,s), z(n,s)) can be obtained by conformal mapping using the 
Schwartz-Christof fel transformation (reference 54). If a curved wall duct is 
approximately represented by straight line segments in the W plane to form a 
many sided polygon, the Schwartz-Christof fel transformation may be used to 
transform the interior of this polygon into the upper half of the C plane 
(figure 24). Then a point source at the origin of the C plane transforms into 
the potential flow through the duct in the W plane. The only approximation 
used in this approach involves the segmentation of the geometry, which would 
occur in any finite difference solutions of equations (105) and (106). 

The complex potential for a source at the origin of the C plane is given 
by 


F = in£ = s + in 


( 112 ) 


where 

(113) 

(114) 


— C0< s < ® 


0 <n < 7T 


In the W plane, using the Schwartz-Christof fel transformation, we have 

N 


^.4- n (C-d,)'V* 
dC C 1*1 


(115) 


where bj is the location of the Ith pole in the € plane corresponding to the 
Ith corner of the polygon in the W plane. In equation (115), the a^'s are 
known, but the bj's must be chosen to fit the duct. If the bj's are known, 
equation (112) and equation (115) may be combined to obtain 


51 


(116) 


dw _ d(z+ir) _ dw_ /dF 
dF ’ d(s+in) dC / dC 

Integration of equation (116) with n constant, produces a streamline, and 
integration of equation (116) with s constant, produces a potential line. The 
metric scale coefficient is equal to the inverse of the magnitude of the 
potential flow velocity and is given by 


V = 


dF 

dw 


. Jil . 
dw dw 


(117) 


The unknown b^'s are determined iteratively by comparing the calculated 
length along each wall with the actual length. In the £ plane, the inner wall 
cooresponds to n * 0 and the outer wall to n = n as shown in figure 25. 
However, since the poles also lie on the C axis, the numerical integration of 
the wall contour using equation (116) is obtained by setting n * e and it = e. 
Hence for bj assumed, the calculated length along each wall for the vth 
iteration and Ith pole is 


u 




(118) 


The actual length along the wall to the Ith pole is known and given by 


Li 



(119) 


and the error in calculating the location of the Ith pole is 





( 120 ) 


Therefore, to minimize the error by Newton's method, we have from equation 
(118) and equation (119) the recursion formula 
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S I M * , »S I M + V I W (L I -L I W ) 


( 121 ) 


where S I V+ ^ is the location of the pole in the C plane in terms of (n, s) 
coordinates. The value of is obtained from equation (112), 


bj* 1 * explsj"*' )cos(n) (122) 


where n * 0 for the inner wall and n * if for the outer wall. If initial 
values for by are assumed, equation (118) may be integrated along each wall 
to obtain Lj and Vj V . Then, equations (121) and (122) are used to obtain 
the next guess. This procedure is repeated until <_ e Q . 

Equations of Motion 

The equations of motion may be written in the streamline coordinate 
system (n, s, $) with the corresponding metric scale coefficients (1/V, 1/V, 
r). Assuming that the plane potential flow streamlines approximate the real 
flow streamlines, the thin channel approximation (reference 55) may be made by 
taking 


U n /U s :0(«) 


(123) 


and dropping terms of 0(e 2 ). 


The resulting equations of motion are given by 


d± 9 

an ‘ V 


(124) 


v_ au» v_ _au 5 _ iL + iP- 

9 an as ~ 9 as an " r as as" 


V 

<7 




(125) 


f 1 


V 
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= 0 


(126) 



3V 3r 

an ” r an 


ap 

an 


v_ _at _ x. il dUs * ? u » u * ii- 

g an as ” g as an r • ds’ 

V__a_ / gr n » \ r n ^ ar t h (127) 

g an \ V / r an + \f 


V d* 31 . V «» « , V a f«nn\ , + i»( U8) 

g an as g 3s an 8 an ' v ' V V 


MT £ <VU,) * T„ (129) 

(130) 

Cp£ ( V £H 

p 8 pRT 


(131) 


(132) 


I- i r * Cp ln(T/T r ) - Rin (p/p r ) 


(133) 


where the stream function ifr has been defined to satisfy the continuity equa 
tion and the normal velocity U n is given by 
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(134) 


^ - . 9 ^ U n 
ds “ v 


Equations (124) through (133) form a set of eight first order partial differ- 
ential equations and two algebraic equations which may be used to solve for 
ten unknowns. The boundary conditions for this problem are given by 


u s (0 t s) = 0 
u$>(0 # s) = 0 

> 

Qn(0,s) = 0 

^( 0,$)*0 


(135) 


for the inner wall and 


U s (7r,s) = 0 (Internal flow) 

= u« (External flow) 

U^(tt,$) 8 0 


(136) 


Qn( 7r » s)= 0 

^(7r,s)=^(tr) 

for the outer wall. 

The remaining unknown is the viscosity. For laminar flow, Sutherlands 
law is used to obtain the molecular viscosity p(T). For turbulent flow 
empirical relations are used to obtain the effective turbulent viscosity Uj.. 
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Turbulence Model 


The turbulence model used for these calculations is an equilibrium 
turbulence model in which the turbulent stress is proportional to the local 
rate of strain of the mean flow based on simple mixing length concepts. Two 
turbulent regimes are used according to the observed characteristics of tur- 
bulence described in reference 56. In regions far removed from the wall 
(freestream), the turbulence is assumed "wake like" with an eddy viscosity 
which is a function of s only. In regions near the wall, the turbulence is 
strongly affected by the wall, and a wall turbulence model is used. 

In the freestream two scale factors are required to describe turbulence; 
a length scale and a velocity scale. The length scale, according to mixing 
length theory, is proportional to the largest turbulent eddy size. This 
length scale is taken to be the duct height. The velocity scale is a measure 
of the kinetic energy of the turbulence and according to mixing length theory, 
reference 57, is proportional to a velocity difference. For duct flows with 
arbitrary velocity distributions, this velocity difference is taken to be the 
difference between the maximum and mean velocity. Hence the freestream 
turbulent viscosity is given by 


b I** 


P-U 

2 



(137) 


where Poo and U*, are the maximum values at a given s. Since the mean 
velocity U is defined by 



(138) 


equation (137) reduces to a wake mixing length model given by: 


Mi ' x - £ F L — tu — °) 


(139) 


where x is an empirical constant determined by experiment. 
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Near the wall, the turbulence model derived by Van Driest (reference 58) 
is used. 


M w 


= i + k 2 Y* 2 [i-exp(-Y7A*)j 


(140) 


where U + and Y + are defined by 


u/'^r U4D 

Y *= P„vVT w /p /fi w (142) 


Two empirical constants are used < and A*** which are taken to be 0.40 and 26.0, 
respectively. 

The empirical constant x remains to be determined. It is noted that two 
boundary layers exist in the duct. Hence for uniform freestream flow, equa- 
tion (137) reduces to 


Mt sX A. u - 5 * 


(143) 


where 6* is the boundary layer displacement thickness. The value for x has 
been obtained by Clauser (reference 59) and is taken to be 0.016. 

Initial Conditions 


The initial conditions for this problem require specification of all 
dependent variables i|>, U g , U^, P, p, T, I, T ng , t., q n as a function of n. 

In the freestream, it is assumed that T ng , T n ^, and q n , flow swirl angle are 
zero. Then the remaining variables are uniquely defined by P Q (n), 8(n) , 

T Q (n), iKir) using the isentropic flow relations and equations (124) and 
(126) to insure that the continuity equation and radial momentum equation are 
satisfied. 

In the nacelle boundary layer, T ns (n), T n <|,( n )> q n (n) mu8t specified in 
such a manner as to be consistent with the turbulence model. This is 
accomplished by using Coles' law (reference 60) for the mean velocity profile 
assuming colateral flow (6 s constant). Thus we have 
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u '* u c +2n c sin 2 (Z- |) 


(144) 


with Coles' shape parameter II c given by 


n c 



(145) 


The function U c + (y + ) is the contribution from the law of the wall. This is 
obtained by integrating 


— ■ — *i 
w dY * 


(146) 


using equation (140) for the turbulent viscosity p^.. The friction velocity 
U is obtained from Coles' friction law 


U. 




+ 2.2 + 2 


(■ 



(147) 


If the displacement thickness 6 and shape factor H 12 are input parameters 
describing the boundary layer we have 


!_ H,2 + I 

8* S Hjg — I 0*8) 


* 4 ^ 

and Coles' law (equation (147)) can be evaluated for U . 

The temperature distribution is obtained from the velocity distribution 
using the relation 
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(149) 


_T_ 

T- 


♦ 


p •/* Izl 

' r 




With the mean boundary layer profiles calculated, x ng , x^ and q n are obtained 
using the relations for turbulent viscosity Up and stress (equations (129) 
to (131) and (140)). 

Inviscid Flow Calculation 


The solution for the inviscid flow around the nacelle is used to calcu- 
late the propeller wake corrections due to the influence of the nacelle on 
the propeller wake and is also used to calculate the influence of the nacelle 
on the flow at the lifting line. 

For inviscid flows with no work or heat addition, the total temperature 
T q , total pressure P Q , and angular momentum (rU^) are conserved stream- 
lines. These quantities, however, may vary from streamline to streamline. 
Therefore along a streamline the isentropic flow relations may be used. 
Furthermore, since the coordinate n is the stream function for the potential 
flow around the nacelle, the curve n(r,z) * n Q approximates the inviscid 
flow streamlines. Hence P Q , T Q , and (rU^) are functions of n only. 

Thus the following equations can be used to calculate the solution for the 
inviscid flow. 


T 0 (n) 
T(n,s) * 


I + 


r- 1 
~r 


M‘ 


( 150 ) 




B 

2 


M 



( 151 ) 


rU^ * r,U^ ( (n) 


( 152 ) 
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In addition the inviscid flow must satisfy the continuity equation, equa- 
tion (124), and the normal momentum equation, equation (126). 


d± . s 

dn “ v 


(153) 


dp . / ?U S av u <ft 2 dr 
dn " V p r dn 


which may be expressed as 



(154) 


(155) 


dP 

. 2 r 

cos 2 a dv sin 2 a dr 

D J 

1 


dn 

" r- 1 

V an r an 

r ( 

l U) 1 

j (156) 


The solution of equations (153) through (156) produces the inviscid 
compressible flow field. This solution is obtained by the following 
transformation of variables. 

Let 


X = 



(157) 


and 


a = 2 


Y’ COS 


dv 

an 


a tt + T 


sin 2 a 


dr 

dn 


(158) 
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Equation (109) becomes 


3* + { 0 + ¥sr| ln(Po, l} x '° 

which is an ordinary first order linear equation. 

The solution is given by 


<*>(n)=expj o /ady + ^ *"(^h)| 


n 

4>(n)= / a<£(y)dy 

0 


X=(x H + 4>(n))/tf>(n) 


The boundary conditions for the solution of the problem are 

*( 0 ) = o 

p (0= P© * 

For the case of the nacelle in a wind tunnel, the boundary 


^(o) * 0 » 

V'(n): 'J'n * 


since mass flow is conserved in the tunnel. 


(159) 


(160) 


(161) 


(162) 


(163) 


conditions are: 


(164) 
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Perforated Wall Bleed Model 


For some tunnels, the use of a zero normal flow boundary condition is 
inappropriate when the tunnel is run using bleed. This type of operation can 
be modeled as described below. 

If one treats a single hole in a perforated wall as an orifice, then the 
mass flow can be derived in terras of the plenum stagnation conditions and the 
local static pressure inside the tunnel, Holman (reference 61). Then an 
expression for the mass flow added to the tunnel flow is given by 




A h 

C 

A_ 




✓YRTa 




(165) 


where P-j and T<j< are the plenum conditions, P is the local tunnel static 
pressure, A^/A g is the ratio of the hole area to surface area, and C the 
effective discharge coefficient which is a property of the perforated wall. 

If the tunnel static pressure is greater than the plenum total pressure, the 
mass flow bleed is out of the tunnel. Under these conditions, Pj and T<j< 
are taken from the wind tunnel conditions, and P is the plenum pressure which 
is assumed known. 

The mass flow bleed is related to the stream function by 


i* .£<£H A 

3s V P r U r 


(166) 


Equations (165) and (166) provide the boundary conditions for a perforated 
wall relating two dependent variables and P in terms of the characteristics 
of the perforated wall and the plenum conditions. 

Since perforated walls have very high rates of bleed which suck off or 
blow off the boundary layer, the boundary conditions were taken to be 


* 


J+l 

KL 



V p r 0 r 


As 


(167) 


ns 


KL 


* 0 


(168) 


Tn *KL 


0 


(169) 
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Nacelle Wake Corrections 


The presence of the nacelle can significantly alter the velocity distri- 
bution at the propeller blade. The flow, and hence the streamline, can 
accelerate or decelerate to follow the shape of the nacelle. Thus if a wake 
model is known for an isolated propeller, then the isolated wake should be 
corrected for the presence of the nacelle by addition of the nacelle's 
influence on the streamline due to its presence in the flow (figure 18). 

The nacelle wake corrections are determined by comparing the calculated 
streamlines for the flow over the nacelle to the classical wake streamlines. 
Since the vorticity shed from the lifting line is convected by the stream- 
lines, the location of a vortex element is given by 


z = Z 0 + u co t 


(170) 


r 


= r„ 


(171) 


* 0 + At 


(172) 


for the classical wake where U,,, is the freestream velocity and 0 the propeller 
rotational speed. The location of the vortex filament in the flow field about 
the nacelle is given by: 


z = z 0 + /u 2 dt 

O 


(173) 


r 


r 0 + /u r dt 
o 


(174) 


\i> z 4/ + 
T o 



(175) 
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The wake corrections (displacement from classical wake) for each streamline 
passing through the lifting line at the point r Q is given by 


Az(r 0 ,t) =/(U z -u 00 )dt (176) 

■'o 

Ar(r 0 ,t) =/u r dt (177) 

•'o 

A = f Q (-^ -n)dt (178) 


Since the lifting line theory uses i|» as in independent variable referenced 
to the blade at t = o (i|> Q = o) . 



(179) 


where the dimensional time is 


dt = 


ds 

vu s 


(180) 


and V is the metric scale coefficient. 
Blade Force 


The action of the propeller blades on the circumferential average flow is 
treated by using body forces (force/volume) in the streamwise and tangential 
momentum equations, equations (125) and (127) and by using a dissipation func- 
tion in the energy equation, equation (128). A rigorous derivation of these 
terms is given in reference 5. 
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The forces produced by the propeller are calculated using lifting line 
theory as derived in the propeller portion of the report. These forces 
(force/span) act on the lifting line and are resolved in the axial, tan- 
gential, and radial direction so that thrust and torque (power) can be calcu- 
lated. Resolution of these forces in the (n, s, <J>) directions is given by 


f s = f* COS 6 +f r sin 9 (181) 


f n = f z sin0 - f r cos0 ( 182 ) 


f *“" fp (183) 

where 6 is the angle that the coordinate streamline makes with the z axis 
(figure 25) and where it is noted that (f g , f n , f^) are forces applied by the 
propeller to the flow field and (f r , f 2 , f p ) are the forces applied to the 
propeller by the flow. An assumption that the dissipation of energy is small 
compared to the work done is made, and thus the rate of dissipation of 

energy of the blade boundary layer is neglected. Furthermore, it is assumed 
that the normal force is small. 

The body forces (force/volume F g , F^) are calculated from the lifting 
line forces (f g , f^) by assuming uniform distribution over the blade chord 
passing through the lifting line and projected onto the (r,z) plane. Thus 


Ps 



u 

c cos a 


(184) 


iL- f * 

ds " c cos a 


(185) 


where a is the angle between the chord line and z direction. 
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Nacelle Drag 


The nacelle drag is composed of two parts; pressure drag and friction 
drag. The pressure drag is the resultant force in the axial direction 
produced by the static pressure acting on the body surface area. This force 
Dpr is given by 


Dpr 



(186) 


where ? is the static pressure acting normal to the body surface and d£ z 
is the differential area in the axial direction. For axisymmetric bodies, as 
shown in figure 27, equation (186) becomes, 


Dpr(x) = 7 T r , 2 P| 


+ 


/ 


x 

P sin0 27rrdx-7rr 


2 


P B 


(187) 


The first term is the contribution due to the pressure acting on the area 
within the centerbody. The second term is the contribution of the static 
pressure acting on the nacelle surface and the last term is the contribution 
of the base pressure acting on the nacelle afterbody. Analytical treatment of 
the separated flow on nacelle afterbody is not sufficiently developed to 
predict the base pressure drag. Therefore it is assumed that the base pres- 
sure P R equals the local static pressure P. 

Following the usual convention, the drag coefficient is defined using the 
maximum projected area in the axial direction. 



D pr 

^Uflo(27rrmox^ 


(188) 


66 



The nacelle friction drag is produced by the nacelle boundary layer. 

This is determined by integrating the axial component of wall stress acting on 
the wetted area on the nacelle as shown in figure 27. 


r* 

Dfr (x) = / T ns cos0 27rrdx 


(189) 


The friction drag coefficient is then given by 


C 


Dfr = 


Dfr 

J ? U <X)( 27rr rox) 


(190) 


where the projected area rather than wetted area is used to keep it consistent 
with the definition of pressue drag coefficient. 

Numerical Solution 


With the relationship between and the mean flow specified, equations 
(124) through (133) can be solved by an implicit forward marching numerical 
integration scheme. These equations are first linearized by expanding all 
dependent variables in a Taylor series expansion in the marching direction 
(s), and terms of 0 (4s 2 ) are dropped. Finite difference equations are then 
obtained using the two point centered difference scheme of Keller (references 
62 and 63). The resulting matrix equations are (10 x 10) block tridiagonal 
and are solved by block factorization using the method of Varah (reference 
64). 


ftiis numerical solution is second order accurate in the n direction, 
first order in the s direction, and linearly stable. The 4s step size is 
limited not by linear stability conditions, but by the required accuracy in 
the Taylor series expansion in s. 
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Streamline Curvature Analysis 


The derivation of governing code equations which are solved in the ADD 
code is given by Anderson (reference 4). This derivation consists of 
constructing an orthogonal coordinate system from the potential fl*ow solution 
and parabolizing the Navier Stokes equations by assuming that the velocity 
component normal to the inviscid flow streamline (streamwise coordinate) is 
small compared to the streamwise component of velocity. As shown in reference 
4, one important effect of this assumption is to replace the actual viscous 
flow streamline curvature with the approximate potential flow streamline 
curvature. Furthermore, it is noted that the potential flow streamline curva- 
ture was calculated from two dimensional incompressible potential flow 
obtained by conformal mapping. These approximations are valid for a large 
variety of flows (reference 65). However, when the streamline curvature 
departs significantly from that of the incompressible plane potential flow, 
then these approximations lead to errors in the prediction of the flow field 
(reference 53). In addition to a viscous flow analysis, an approximate 
inviscid flow analysis was derived which was based on this streamline curva- 
ture approximation. This inviscid flow analyses is capable of treating 
axisymmetric compressible swirling flow with variations in total pressure and 
total temperature from streamline to streamline (i.e., rotational flow). Both 
the viscous and inviscid analysis are dependent on this streamline curvature 
approximation and show significant errors under the same circumstances. The 
purpose of this section is to derive streamline curvature corrections to 
reduce the error in calculating the solution of both the viscous and inviscid 
flow equations. 

Ba_sic_ E^at_ion£ 

The basic parabolized Navier Stokes equations used by the ADD code 
(reference 4) and including all convection terms are given below. These 
relationships neglect the body forces for the development of the curvature 
corrections and are presented in a normalized form. 


Continuity Equation 


w-° 

(191) 

6PU. , 

V dn 

(192) 

GPl*, i* 

V 0S 

(193) 


S Momentum Equation 




GPU. dU, 
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(194) 
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<t> Momentum Equation 


GPU S *U* GPU, 
v as + V 


GPU n dU* , GPU$U* dR A { GPUnU* dR j . A. 

V dn + VR ds j VR dn j ' R dn V V / 


(195) 


n Momentum Equation 


Energy Equation 


(GPUs dU n 6PU n dUn GPUnU* dV I 
j V ds * V dn V 2 ds ( 

GPU S 2 dV GPU* 2 dR I 6 dn 
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(196) 
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(200) 

Relations 
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Curvature 


Velocity Components 


Vorticity 


1 dV 

K n * t n 

(207) 

K * R s dn 

(208) 

k *'T 

(209) 

U n = U sin Q 

(210) 

U 8 = (ucos 0)cosa 

(211) 

U* s (Ucos0)sina 

(212) 

* r- £ (Ru *' 

(213) 

"» * - t ■& (RU *' 

(214) 

-.■•■liW-iWl 

(215) 


The terms in the brackets in equations (195) and (196), the tangential 
and normal momentum equations, are dropped in the ADD code equations but are 
retained for the purpose of deriving streamline curvature corrections. The 
curvatures are given by equations (207) through (209). K n is the curvature 
of the coordinate potential line; K g is the curvature of the coordinate 
streamline; and is the curvature of the swirl flow. 
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The velocity components are given by equations (210) through (212). It 
is noted that the angles are defined relative to the (n, s, <j>) ADD code 
coordinates. With no swirl, = 0 and 0 is the angle between the velocity 
vector and the streamwise coordinate S. 

The vorticity components are given by equations (213) through (215). If 
there exists variations in total pressure and total temperature from stream- 
line to streamline, then by Crocco's theorem (reference 66), these vorticity 
components are not zero. 

The ADD code equations used in reference 4 are given by equation (192) 
and equations (194) through (202) with the terms in brackets in equations 
(195) and (196) dropped. These equations are based on the assumption that the 
actual streamlines lie close to the assumed streamlines which in this case are 
the coordinate streamlines obtained by a conformal mapping solution of the two 
dimensional incompressible potential flow equations. Under these conditions 
U n is small compared to U s and the terms in the brackets in equations 

(195) and (196) may be dropped. These equations are parabolic differential 
equations (see reference 4) and can be solved by forward marching numerical 
integration. 

If the terms in the brackets of the normal momentum equation, equation 

(196) , are retained, then the set of equations is elliptic. This can easily 
be demonstrated by substituting the equations for the stream function, 
equations (192) and (193) into equation (196). The resulting equation will 
have a 3^/3S 2 term which will make the equations elliptic. If, however, the 
terms in the bracket are treated as known, the equations remain parabolic. 

Streamline Curvature Correction 


When the terms in the brackets of equation (196) are neglected, the 
actual streamline curvature is taken to be the streamline curvature K s of 
the coordinates. Three sources of error are then introduced. First, the 
effects of axisymmetric flow on curvature are neglected and replaced by two 
dimensional effects. For many gas turbine ducts the ratio of inside radius to 
outside radius approaches one and this effect is small. However, for the high 
speed propeller this effect is significant. Second, the effects of compress- 
ibility are neglected. Again, however, in curved ducts as choked flow is 
approached, the effects of compressibility on streamline curvature may be 
significant. Third, viscous effects are neglected. These viscous effects are 
principally due to displacement of the streamline produced by the growth of 
the boundary layer. Near separation, the boundary layer grows very rapidly 
and significant changes in streamline curvature may occur. 
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The normal momentum equation, equation (196) can be written in the 
following convenient form; 


where 


/ K s + AK S \ 2 FU» 2 dR 

“ V V / ru s ' R dn 


s0 

y Mf dn 


(216) 


AK S = -V 
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The term AK s is the curvature of the actual streamline relative to the 
coordinates. It can be seen in equation (216) that the form of the ADD code 
normal momentum equation is preserved. Furthermore, it should be noted that 
if AK g is known a-priori, the ADD code equations remain parabolic. 

The set of equations, equation (192) and equations (194) through (202) 
with equation (216) replacing equation (196) and with AK g included can be 
solved as an elliptic set of equations using a multi pass method by lagging 
AK S . This procedure, however, may incur severe stability problems. An 
alternate procedure, adopted here, is to neglect viscous effects but include 
axisymmetric and compressibility effects by calculating the streamline curva- 
ture of compressible axisymmetric potential flow and using it in place of 
(K g + AK 8 ) in equation (216). Since this streamline curvature is known a- 
priori, the set of equations remains parabolic and can be solved in a single 
pass. 
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If there are no gradients in total temperature and total pressure and no 
swirl, then by Crocco's theorem the flow is irrotational and the vorticity 
components, equations (213) through (215), vanish everywhere in the flow 
field. Under these conditions equations (192), (193) and (215) may be 
combined to yield a single equation for the stream function. 



a / 1 d±\ j_ /j a*\ 

an \PG dn f as Ipg ds/ ° 


(218) 


where the density p 


is given by 
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U 8 and U n are obtained from the stream function using equations (192) and 
(193). For incompressible flow p * poo is constant. Equations (218) and 
(219) replace the two momentum equations, equations (194) and (196). Equation 
(218) is an elliptic equation and requires that ^ be prescribed on the 
boundary. 
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Equations (222) and (223) and all following equations contain a coordinate 
stretching function for the normal coordinate which is given by 
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(224) 


where B is an arbitrary parameter. 

Equations (218) and (219) can be solved iteratively. If p(n,s) is 
assumed known, then equation (218) can be solved by line relaxation. If we 
choose the line of unknowns to be along the normal coordinate (S ■ constant) 
and J is an index for the streamline and K is an index for the normal line, 
then the finite difference equation for equation (218) is given by 
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where 
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(226) 


Treating ip ^ as the unknown and collecting coefficients we have 
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and equation (225) reduces to 
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with the boundary conditions 
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The matrix equation (231) is tridiagonal and can be solved by Gaussian 
elimination using the following algorithm. 
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Then the new guess for (i|» j£) is given by 
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where fl is a relaxation parameter (0 < fl < 2). When an iterative guess for 
the whole field is known ^(n,x), a new iterative guess for the density is 
obtained from equations (192), (193) and (219). This iteration procedure is 
repeated until 
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Once the solution for ♦ is obtained from equation (218), the flow angle 
is determined from equations (210) and (211). 
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The normal pressure gradient is determined by differentiating equation (204) 
where % is constant for the flow field. Then 


an 
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Equating equation (239) to equation (216) where 


0, results in 
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where K c is the streamline curvature for axisymmetric compressible flow 
which is to be used in equation (216) in place of (K s + AK S ). 
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The ADD code contains an approximate inviscid flow analysis which is used 
to obtain blade forces using blade element strip theory (see references 5 and 
65). This inviscid flow theory is derived using the same streamline curvature 
approximation that is used in the viscous flow analysis. Since gas turbine 
flows are swirling compressible flows with variations of total temperature and 
total pressure from streamline to streamline, they are rotational flows. The 
approximate inviscid analysis in the ADD code was developed to treat these 
flow fields. The purpose of this section is to derive the streamline curva- 
ture corrections for this analysis. Later it shall be shown that if exact 
streamline curvatures are used in the analysis, the analysis is also exact. 

For inviscid flow, Q n , E , E^ are zero. Then equations (195) and (197) 
with the use of equations (192) and (193) reduce to 


d* dl d* dl 
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The form of the solution to equations (241) and (242) can be written as 
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Then using equation (205), 


n T *n T (*) (245) 

e T = © T (*) (246) 


Equations (243) through (246) demonstrate the well known theorem that in 
inviscid flow, I, (RU^), Up, Op are constant along streamlines but may vary 
from streamline to streamline. The normal momentum equation (equation (196)) 
is written in the following form 

an 2 9 r Ke 9 i aR . •% "1 2 
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The approximate inviscid solution assumes that the potential flow streamline 
used to construct the coordinates approximates the actual streamline. Hence 


y = n 
6 = 0 
AK$ 3 0 


(248) 


and the equations listed above are easily integrated according to the current 
method. Clearly, if the flow angle 0 and the curvature K c from equation 
(240) replacing (K s + AK S ) can be specified then a more accurate solution 
for the inviscid flow can be obtained. 
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COUPLING PROCEDURE 


Assumptions Affecting the Coupling 

As noted in the Introduction, the coupling of the two separate procedures 
for the nacelle and propeller portions of the analysis requires the transfer 
of pertinent data affecting the respective solution procedures of the 
analysis. The ability to use this solution technique to obtain the combined 
propeller-nacelle performance prediction is based upon the following major 
assumptions related to the coupling process which are noted below. 

As was noted in the technical approach for the propeller solution pro- 
cedure, the flow field is steady with respect to the propeller blades; how- 
ever, this is not true when the point of reference is the nacelle wall 
boundary. Because the nacelle analysis does assume time-averaged properties, 
an assumption concerning the propeller's influence must be made in order to 
use the nacelle analysis and to avoid the necessity to calculate the local 
time varying induced flow field about the nacelle due to the rotating 
propeller and wake. The assumption is made that the inclusion of the effects 
of the propeller on the nacelle and surrounding flow field can be handled in a 
time-averaged manner by including only the addition of the local blade forces 
at the blade row locations and that this model will adequately approximate on 
a time-averaged basis the effects of the combined blade and trailing wake 
system in the flow, to predict the flow field about the combined propeller- 
nacelle configuration. 

To calculate the forces on the propeller blades it is necessary to assume 
that the nacelle flow field solution is representative of the actual flow 
field. Since this is a basic assumption of the nacelle portion of the 
combined analysis, no further explanations are necessary concerning the 
nacelle portion of the analysis in terms of the coupling procedure. The 
assumptions concerning the influence of the nacelle on the wake geometry were 
described in the section entitled: Nacelle Influence on the Wake Geometry, 

and are not repeated here. 
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Description of the Combined Analysis Solution Procedure 


To obtain the combined propeller-nacelle performance solution, the inter- 
facing of the propeller and nacelle solution procedures requires the transfer 
of data between the two portions of the analysis. In the combined- analysis , 
the propeller portion requires that the local inflow conditions (velocity, 
density and speed of sound) at the blade row(s) be specified and that the wake 
geometry be corrected for the presence of the nacelle, if so desired. The 
nacelle portion requires that the lifting line geometry and local blade forces 
(radial, tangential, axial) per unit length be specified for the propeller 
blade row(s). Figure 28 is a simplified flow diagram of the combined solution 
procedure with the required transfer of data noted, as developed for this 
application. The major tasks performed in the analysis associated with the 
pertinent features of the combined analysis are also shown in this flow 
diagram. 

The structure of the solution procedure is shown in the diagram. Within 
the analysis there are two ways in which the combined analysis can be used. 
First, the program can be run with only the inviscid flow field calculated to 
define the inflow and wake perturbations for the propeller portion of the 
analysis to obtain the combined propeller and nacelle performance; however, 
the blade forces are not included in the inviscid nacelle solution. This may 
be desirable for cost considerations if it is shown that the viscous solution 
does not significantly change the propeller performance predictions as 
compared with the inviscid propeller performance predictions, and, if the 
inclusion of the blade force’s influence on the nacelle performance is not 
desired. Second, if the viscous solution is desired, the analysis can also be 
made to cycle through the viscous portion a selected number of times. This 
mode will include the blade force’s influence on the nacelle performance. The 
selection of the number of cycles required to obtain compatible solutions must 
be determined from the user’s experience. Based on the experience gained in 
the development of the computer analysis, it is believed that generally only 
one cycle will be required. However, this assumption remains to be proven in 
actual practice. It should also be noted that for the viscous solution's 
influence at the blade, it is not currently possible to separate, from the 
viscous velocity field solution, the portion of the inflow due to the addition 
of the blade forces. To resolve this problem, the calculation of the flow 
field due to the nacelle at the blade row is performed at the blade leading 
edge in the viscous solution procedure. This approach removes the blade 
force's influence on the flow field at the propeller because of the parabolic 
nature of the problem. That is, the inclusion of the blade forces, distri- 
buted over the blade surface, is not felt at the leading edge of the propeller 
blades, only downstream of this location. 
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INITIAL APPLICATION OF THE ANALYSIS 


The initial application of the analysis was performed under Contract 
NAS3-20961 and was limited to a single propeller-nacelle configuration and 
test condition. The results obtained from this application indicated that the 
inflow velocities due to the nacelle body were incorrectly predicted at the 
blade row. This finding was verified by the use of an alternate nacelle analy- 
sis to predict the velocity profile at the blade row location. In addition to 
this problem, the use of the computer program requires the selection of 
various propeller aerodynamic modeling features which can significantly affect 
the performance predictions. As a result, Contract NAS3-22142 was established 
to refine and modify the analysis to remove known problems and to expand the 
scope of the application of the analysis to include different propeller 
nacelle configurations, test conditions, and modeling features. Following 
this effort additional applications of the analysis have been made under 
Contract NAS3-22257 to configurations and operating conditions for which test 
data is available for comparison. 

In general, the results of the initial application of the analysis will 
not be reported herein because of the limited nature of that study. It is 
sufficient to note that this single application led to the more expanded 
application and refinement effort for the contracts noted above. The results 
of this expanded application are reported below. This study includes a 
comprehensive investigation of the effects of the propeller aerodynamic 
modeling features, the sensitivity of results to the airfoil characteristics 
defined for the same airfoil section from different sources, the application 
of the analysis to different operating conditions, and the application of the 
analysis to different propeller-nacelle configurations which include single 
and coaxial geometries. 
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NACELLE FLOW FIELD PREDICTIONS 


The combined analysis requires that both the flow field around the 
nacelle body and the blade loading solution be predicted with reasonable 
accuracy. Both solutions are predicted based on the super-position of the 
individual solutions in a cyclic manner. The refinements to the original 
axisymmetric duct flow field solver (ADD Code, references 4 and 5) to treat 
axisymmetric nacelle configurations with or without wind tunnel wall influence 
(including tunnel wall bleed) have resulted in an analysis (PANPER) capable of 
treating a wide range of high speed propeller configurations when combined 
with the propeller blade loading analysis. The verification of the ability of 
the analysis to predict the flow field perturbations around the nacelle for 
these high speed conditions was performed by comparing predictions for an 
isolated nacelle at Mach Numbers of .2, .7, and .8 with the predictions 
obtained using an established transonic, inviscid flow field solver for 
axisymmetric bodies (RAXBOD, references 67 and 68). The results of this 
comparison are presented for the .8 Mach Number condition. The results at .2 
and .7 Mach number will not be presented, but are similar in the degree of 
correlation. 

In figure 29, the predicted nacelle wall pressure distributions for the 
SR-3 Nacelle design at .8 Mach Number are shown for both the PANPER and the 
RAXBOD codes. An example of this propeller blade planform is shown in figure 
30 for illustration purposes. Superimposed on figure 29 is the corresponding 
nacelle profile and the approximate location of the blade lifting line. The 
prediction for the pressure distribution using the PANPER analysis differs 
from the RAXBOD prediction near the nose of the nacelle. This is due to the 
use of an optional smoothing process in the PANPER analysis applied to the 
nacelle profile coordinates resulting in slightly different local profile 
coordinates. Further downstream, past the blade row in the region of the high 
curvature of the nacelle profile, the predicted solutions also exhibit some 
differences. The RAXBOD analysis predicts a lower pressure distribution 
corresponding to a region of transonic flow. The PANPER analysis, which is 
not a transonic flow analysis, clearly does not predict the same pressure 
distribution. It should be noted that the PANPER analysis yielded nearly 
identical results for both the inviscid and viscous pressure distributions for 
the other lower Mach number conditions. These comparisons indicate that the 
general character and level of the nacelle wall pressure distributions are 
reasonably predicted, although regions of supersonic flow are not correctly 
modeled. However, this discrepancy is probably not significant for propeller 
performance predictions in the region of the blade lifting line location as 
will be shown in the next section. 
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During the course of this investigation test measurement data for the 
flow field upstream of a non-lifting propeller, obtained in the NASA Lewis 
Research Center test facilities, was made available to UTRC and has been 
presented in reference 69. The results of applying the nacelle code to these 
test results are included herein as an additional example of the code's 
predictive capability. This study also illustrates the need to include the 
effect of blade thickness in the nacelle flow field solution procedure when 
using the analysis to predict the flow field upstream of the blade row. The 
condition investigated was at .8 Mach Number for the non-lifting SR-3 
propeller/nacelle configuration. The flow measurements were obtained using 
laser velocimetry techniques at an axial position .09 propeller radii upstream 
of the propeller blade row. Three different nacelle flow field solutions were 
obtained to demonstrate the effect of different streamline curvature models 
and the presence of the blades in the solution procedure. The results in 
figure 31 present the measured and predicted Mach Number distributions along 
the blade lifting line location for this condition. The results of the 
different streamline curvature models, without blade thickness, illustrate the 
improvement in the predictive capability of the nacelle code when the two- 
dimensional streamline curvature model of the original analysis was modified 
to an axisymmetric streamline curvature model for propeller application. The 
inclusion of the blade thickness in the axisymmetric streamline curvature 
model further improves the flow field predictions near the high solidity 
region of the propeller blades where the presence of the blade should be most 
strongly felt. The discrepancy in the overall level of the flow field (Mach 
Number) away from the nacelle wall may indicate that the actual condition 
tested is not a true 0.8 Mach number condition. If the analysis were run to 
match the implied freestream Mach Number of the test results, the correlation 
would be excellent over most of the range of radial positions measured. Near 
the nacelle wall, the results would differ slightly. This difference could be 
due to either the numerical model or a slight difference in the location of 
the measured data compared with the position used in the analysis. Another 
possible reason for the difference is that the SR-3 propeller is never truly 
non-lifting in a local sense. The blade's integrated lift may be zero, but 
the local blade loading can be non-zero. Thus, there could be the effect of a 
local lift induced flow field in the test results. The analysis was run 
without the blade forces, and thus there are no lift induced velocities 
included in the predictions. 
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The use of the nacelle induced flow field at the propeller blade as 
predicted by the nacelle portion of the PANPER analysis requires that the 
locations along the blade at which the velocity components are calculated 
be specified. Traditional Prandtl lifting line theory under uniform free- 
stream conditions requires that the angle of attack at the wing quarter-chord- 
line be defined using the freestream velocity and the induced velocity. The 
propeller lifting line solution is an extension of this concept to the 
propeller environment which includes the rotational and nacelle induced 
velocity fields. Near the nacelle, the nacelle induced flow field varies in 
the radial and axial directions (axisymmetric assumption). Thus, the nacelle 
induced flow field used in the propeller portion of the analysis is affected 
by the selection of the location (blade quarter-chordline, leading edge, etc.) 
at which this information is calculated. The selection of these locations has 
not yet been clearly defined. In figure 32, the axial component of the 
nacelle flow field, as predicted using the nacelle portion of the PANPER 
analysis, is presented along the blade leading edge and the blade quarter- 
chordline with and without the inclusion of the blade thickness in the stream- 
line curvature model. The results indicate that the axial velocity is rela- 
tively insensitive to the selection of the prediction location for the solu- 
tion obtained neglecting the blade thickness influence in the streamline 
curvature model. However, when the blade thickness is included in the 
solution procedure, the axial velocity component shows a strong sensitivity to 
the selection of the prediction location. It should be noted that the radial 
velocity component does not exhibit this strong sensitivity as shown in figure 
33. 


The sensitivity of the axial velocity to the selected prediction location 
introduces an additional uncertainty in the use of the analysis. Since the 
propeller portion of the analysis already recognizes the high solidity regions 
of the blade near the hub walls through the use of the cascade characteris- 
tics, the use of the blade thickness influence in the prediction of the 
nacelle induced flow field solution is actually "double accounting" for the 
blade thickness. Thus, for propeller performance predictions, the most appro- 
priate way of using the analysis may be to either; include the blade thickness 
in the nacelle induced flow field solution without using cascade charac- 
teristics in the lifting-line propeller solution, or conversely, use the 
cascade characteristics in the propeller solution without the blade thickness 
influence in the nacelle flow field solution. However, for the first method 
the selection of the prediction location (leading edge, quarter chordline, 
trailing edge) for the nacelle induced velocities as used in the propeller 
analysis has not been adequately investigated. As a result, the analysis is 
generally run for propeller performance predictions using the second method of 
neglecting the blade thickness influence in the nacelle induced flow field. 

For nacelle induced flow field predictions away from the propeller disk, the 
nacelle analysis should be run including the effects of blade thickness as was 
done for the results presented in figure 31. This is necessary because the 
lifting-line models do not account for thickness in the propeller induced 
field. 


85 


SINGLE PROPELLER PERFORMANCE PREDICTIONS 


Aerodynamic Modeling Features 

The development of the propeller portion of this analysis is* based on the 
application of methodologies which have been extended to the limits of their 
operational applicability to attempt to provide a useful, numerically 
efficient performance analysis. To provide this type of analysis tool using 
the incompressible lifting line model required the inclusion of many different 
aerodynamic modeling features. These features attempt to include in a 
rational manner the effects of compressibility on the propeller tip loading, 
the effects of compressibility on the propeller wake induced velocity field, 
the influence of the nacelle flow field, the effects of the high solidity 
inboard blade sections using cascade data or corrections to isolated airfoil 
data for the cascade influence, the effects of highly swept blades, and the 
selection of different sources of isolated airfoil and cascade character- 
istics. This large number of modeling features has resulted in an even larger 
number of combinations of these features which can be used when predicting the 
propeller performance. In addition, these features do not address the 
question concerning the inclusion of the blade thickness influence in the 
nacelle flow field solution procedure noted earlier. Obviously many combina- 
tions are unrealistic, however there are still a large number of realistic 
combinations which can be used. This investigation has attempted to define a 
reasonable combination of these features by applying the analysis to different 
propeller-nacelle configurations and operating conditions to study the 
results, examine the sensitivities to modeling options and compare the predic- 
tions with test data where available. These comparisons are presented to 
provide some useful insight into the best selection of these modeling 
features . 

Through out this report the terms net and apparent efficiency are used. 
For clarity they are defined herein. The net efficiency (n) is that measured 
in the tunnel by the force balance for the total system installation. In a 
simplistic sense this is the total of the propeller forces and nacelle body 
forces. The apparent efficiency (n^) i8 the component due only to the 
propeller forces generated in the presence of the nacelle. Since the effect 
of the actual body is to reduce the efficiency, the apparent efficiency is 
usually higher than the net by some small correction (Ari). 


n A - r\ + An 
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The analysis can only predict the apparent efficiency since the nacelle body 
force predictions are valid only up to the termination of the body model or 
the separation point. Since separation usually occurs for the afterbodies 
used for these studies, the model was truncated reasonably far downstream from 
the blade row, but did not have closure on the geometry. Estimates for the 
back pressure can be made, but they are not of sufficient accuracy for this 
study. For comparisons without test results, the apparent efficiencies are 
used. In this report, the choice of the use of apparent or net efficiencies 
for comparison with test data was based on the effort required to adjust 
either the data or the predictions. Where the amount of test data was small 
compared to the number or prediction cases run, the data was correct to 
apparent efficiency. For the cases where the amount of data was larger, the 
predicted efficiencies were corrected to net values. These corrections were 
supplied by NASA. 

The results presented herein start with the application of the analysis 
to selected propeller and nacelle configurations to compare the different 
aerodynamic modeling features. The influence of the nacelle on the solution 
is neglected for these preliminary applications. The intention is to define 
and illustrate the sensitivity of the propeller performance predictions to the 
modeling features. The analysis is then applied to conditions for which test 
data are available. For this application the effects of the different 
modeling features are again studied with the number of features limited in 
scope based on the first series of applications. However, because test data 
are available, the nacelle influence is included for this study. And finally,' 
the analysis is applied to configurations in selected wind tunnel facilities 
to identify any potential problems in the use of the different facilities with 
the model scaling selected. 


Single Propeller-Nacelle Configurations 

Fundamentally the single propeller-nacelle configurations studied 
represent variations in blade design and nacelle shape. The variation in 
blade design include blade sweep, taper, twist, solidity, and thickness 
changes for the different blade geometries. The nacelle geometries studied 
were either the area ruled or conical shaped profiles. These nacelle profiles 
are shown in figure 34. The propeller geometries were the NASA designated SR- 
2, SR-1, SR-1M, and SR-3 designs; representing respectively, straight to 
moderately swept blades (0, 30, and 45 degrees). The difference between the 
SR-1 and the SR-1M designs is the blade tip twist and section camber. The 
airfoil sections were the NACA Series 16, transitioning to series 65 near the 
blade root. The blade design characteristics for each propeller configuration 
are presented in graphical form in figures 35 to 37. The plan forms were 
shown earlier in figure 30 for these designs. 
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Single Propeller-Nacelle Operating Conditions 

The range of operating conditions studied was based on the range of early 
unpublished test data provided by NASA. The range in free-stream Mach number 
was from 0.6 to 0.8. For the initial studies, which investigated, sensitiv- 
ities to the modeling options, the variations in rpm at a fixed power loading 
(constant power) provided advance ratio variations in the range of 2.6 to 3.2. 
For the performance map predictions shown later, the advance ratio variations 
were in the range of 2.5 to 4.5 with blade angle variations in the range of 58 
to 63 degrees depending on the propeller configurations and free-stream Mach 
number conditions (0.6 and 0.8). The range in power coefficients studied has 
from near windmilling conditions to as high as 3.0 depending on the combina- 
tion of conditions. The design condition, when referenced herein, corresponds 
to an advance ratio fo 3.07 with a power coefficient of 1.69. 

To help in determining when operating conditions represent supersonic tip 
speeds, the following relationship can be used to find the advance ratio which 
corresponds to a tip Mach number of 1.0. 


“Tip 




1 + 




1/2 


- 1/2 
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For the 0.6 free-stream Mach number condition there was no advance ratio 
condition studied which had supersonic tips. At the 0.7 Mach number 
condition, supersonic tip conditions existed for advance ratios less than 
3.079. However, at the design point advance ratio (3.07) for this Mach 
number, the fraction of blade radius which experiences supersonic flow is very 
small (0.003R). For the 0.8 Mach number condition all advance ratio 
conditions less than 4.189 had supersonic tip Mach numbers. At the design 
point the outer 25 percent of the blade is operating in supersonic flow. 


Effect Of Compressible Tip Loss Models 

The analysis currently has four compressible flow tip loss models 
available along with the non-tip loss model. These models were described in 
an earlier technical section, but for convenience they will be briefly 
described here. The first model is the Evvard model (tabulated solutions) 
which is based on fixed wing, rectangular tipped geometries applied to the 
non-rectangular geometries of the high speed propeller designs. The second 
model is an exact analytical integration of the relations used in conical flow 
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theory for the rectangular wing tip. In theory this model is identical to the 
Evvard model which is available in tabular form, however there are small 
differences in practice. The third model is based on conical flow theory for 
swept wings assumed to be operating at a free-stream Mach number equal to the 
propeller tip Mach number. The fourth model is a variation of this model 
where the tip loss correction is further scaled using the ratio of the Prandtl 
Glauert compressibility corrections based on the propeller's local section and 
tip. Only the tabulated Evvard and the swept conical flow models were studied 
in depth and Mach numbers compared with the non-tip loss model results. The 
analysis was applied using the SR-1, SR-2 , and SR-3 propellers with the area 
ruled spinner and nacelle geometry at; (1) the design advance ratio for a 
variation in free-stream Mach number at the design power coefficient, (2) the 
design advance ratio with a variation in the design power coefficient for a 
fixed Mach number, and (3) the design power coefficient for a variation of the 
advance ratio at a fixed Mach number. For some of the conditions investigated 
using the swept conical flow tip loss models, solutions were not obtained 
because the combination of blade geometry and operating condition violated 
assumptions used in the model. For this study, the analysis was applied using 
the cascade data correlations from reference 40 and the cascade correction 
based on flat plate theory was applied to the isolated airfoil data of Manoni 
over the blade transition region. The compressible wake model was used for 
all of the above cases. In addition, the influence of the nacelle was 
neglected for this study and thus the predicted results should not be 
considered indicative of the actual propeller performance capabilities. The 
results are presented separately in figures 38 to 40 for each propeller 
configuration. In these figures, the propeller efficiency is plotted on a 
greatly expanded scale as a function of the free-stream Mach number to 
illustrate the differences in the predictions using the various tip loss 
models . 

The results of the application of the analysis to the unswept propeller 
blade design (SR-2) are presented in figure 38. At the 0.7 Mach number (tip 
Mach number 1.005), all tip loss models appear to function in the analysis, 
but the results are identical. This is because the supersonic region is too 
small (0.3 percent of radian) relative to the size of the lifting line 
segments along the blade and does not activate any of the models. The swept 
conical flow tip loss models are not used at the 0.8 Mach number for this 
propeller geometry because the assumptions used to develop the model are 
violated; the trailing edge of the blade must be supersonic. The effect of 
the Evvard tip loss model for all of the conditions at this Mach number is to 
reduce the propeller efficiency for a specified power level. When the 
analysis was run at a specified blade pitch the same effect was noted, 
although these results are not presented. 
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For the SR-1 propeller, a 30 degree swept blade design, the results of 
the application of the different tip loss models are shown in figure 39. 

Again, the lower Mach number condition (0.7) does not activate any of the 
models and the swept conical flow tip loss models can not be applied for this 
design at the 0.8 Mach number conditions. The Evvard model again results in 
the reduced efficiency compared with the non-tip loss model. 

For the SR-3 propeller, a 45 degree swept blade design, the results of 
the application of the different tip loss models are presented in figure 40. 

At the 0.7 Mach number the swept conical flow models are again not applicable 
because the supersonic trailing edge requirement of the model is not met. At 
this Mach number the difference between the non-tip loss and the Evvard tip 
loss model results are small. The results for the application at 0.8 Mach 
number vary depending on the loading condition and advance ratio. The reduced 
loading condition (Cp*1.5) at the 3.07 advance ratio exhibits no significant 
difference between the swept conical flow tip loss models and the non-tip loss 
model while the use of the Evvard model results in a reduced efficiency 
prediction. At the design loading condition (Cp*1.7) for this advance ratio 
the swept conical flow tip loss models do result in a difference in the 
predicted efficiency at the 0.8 Mach number. The conical flow model, with the 
variable Mach number distribution scaling, results in relatively little effect 
on the performance. For the lower advance ratio case at this Mach number the 
use of the conical flow models results in reduced efficiency, and both models 
yield about the same level of reduction. For all cases at this Mach number 
the Evvard model results in reduced efficiencies. 

As a result of this study the decision was made to only use the Evvard 
tip loss model in the applications to follow. This decision was made because 
the conical flow model is sensitive to propeller trailing edge geometry and 
Mach number as evidenced in the above results. This makes it unacceptable for 
general application. This selection of the Evvard model is not meant to imply 
that it is better than the swept conical flow tip loss models, but that as 
implemented, the conical flow tip loss model is not reliable as an investiga- 
tive tool for all of the geometries studied in this investigation. 


Effect of Cascade Models 

The modeling of the propeller blade lift and drag characteristics is a 
critical aspect of the propeller portion of the analysis. The performance 
predictions can be no better than the lift and drag characteristic data used 
in the analysis. Within the analysis the inboard blade regions (NACA series 
65 airfoil sections) operate at large chord-to-gap ratios for the high speed 
propeller designs under investigation by NASA. Isolated airfoil character- 
istics do not represent these cascade type blade sections. Two different 
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sources of cascade data have been provided in the analysis to model these 
blade section lift and drag characteristics. Both models are correlations, 
the first based on the NASA method of reference 40, the second on the method 
by Mellor (references 49 to 51). In addition the analysis has corrections 
which can be applied to isolated airfoil data based on either flat plate 
cascades or the Mellor correlation. 

To study the influence of the cascade data on the performance predic- 
tions, the analysis was applied to the SR-1 propeller geometry using both of 
the cascade correlation models at the 0.8 Mach number condition for a 
variation in the advance ratio with a fixed blade pitch angle (60.2 deg.). 

For both cascade correlations used in the inboard regions, the flat plate 
cascade correction was used to correct the isolated airfoil characteristics 
along the blade transition region to attempt to model the airfoil character- 
istics in this region. In addition, a case was run with the isolated airfoil 
characteristics in this region corrected based on the Mellor correlation. 

(The inboard and transition regions are shown in figures 35 through 37.) This 
case used the Mellor correlation model, resulting in three different cascade 
characteristic models for comparative purposes. The variation in rpm 
represents the range of advance ratio from near 3.0 to windmill at a fixed 
Mach number. Again, no nacelle influence was included for this study. The 
predicted propeller efficiency, power coefficient, and thrust coefficient are 
presented respectively in figures 41 to 43 for these conditions. The plot of 
efficiency versus propeller rpm (figure 41) indicates a maximum difference 
between the predictions for the three cascade models of about 2.5 points in 
efficiency except at the lowest rpm studied. At this windmill condition the 
difference in the predicted efficiencies is quite significant because the 
propeller is essentially developing zero power as shown in figure 42. The 
outboard regions of the blades are so lightly loaded in relation to the 
inboard regions that differences in the inboard cascade data result in very 
significant differences in efficiency. For the more highly loaded conditions 
near the design advance ratio the outboard sections carry most of the 
airloading and thus the differences in the cascade data on the inboard regions 
do not have as significant of an effect. At these higher rpm conditions 
(lower advance ratio) the results obtained using the different cascade models 
do not show a consistent behavior in terms of the differences in predicted 
efficiencies for the variation in rpm; except that the efficiencies obtained 
using the NASA correlation model are always the highest or equal to the 
others. Over this rpm region the average variation in predicted efficiency is 
just under 1.0 point. The predicted power coefficients using these different 
models are shown in figure 42 to illustrate the windmill effect as the 
propeller rpm is reduced. The variation of the power coefficient with 
propeller rpm is more consistent than the efficiency was. For completeness, 
figure 43 shows the variation of the thrust coefficient with rpm for these 
different cascade models. Note that at the windmill condition the propeller 
is creating a drag force. 



Effect Of Transition Region Interpolation Model 


The transition region on the propeller blade represents the change from 
the NASA 65 Series airfoil sections on the inboard region of the blades to the 
NASA 16 Series airfoil sections on the outboard region. The analysis can 
currently model the airfoil characteristics in this transition region by two 
methods. The first method assumes that the airfoil characteristics in the 
transition region can be modeled by using the isolated NASA 16 Series airfoil 
characteristics scaled by an analytic function based on either flat plate 
cascade theory or the Mellor correlation to incorporate the cascade influence 
as a function of the local gap-to-chord ratio. Thus with this model the 
cascade influence diminishes concurrently with the transition to the outboard 
section characteristics. However, this model ignores the fact that the 
airfoil section in this region is neither a 16 or 65 NACA series, thus any 
cascade correction to the isolated 6 series data is not a true representation. 
The second method attempts to rectify this deficiency. It is based on an 
interpolation function which uses the inboard cascade data and the outboard 
isolated airfoil data as end point conditions to provide a functional 
relationship between spanwise location (gap-to-chord ratio) and section 
airfoil characteristics. In the analysis the behavior of the interpolation 
function can be varied via input of the exponent on the gap-to-chord ratio as 
noted in an earlier section of this report. The sensitivity of the propeller 
performance predictions to the modeling of the transition region character- 
istics is also shown in figures 41 to 43. In figure 41 from the previous 
section, the predicted propeller efficiency is presented as a function of 
advance ratio for the SR-1 propeller design using different cascade correla- 
tions along with a variation in the transition region model. The variation in 
the transition region model is shown between the flat plate and Mellor model 
correction to the isolated airfoil characteristics using the Mellor cascade 
correlation for the inboard cascade data. The two different correction models 
result in relatively small differences in the predicted efficiencies, except 
at the windmill condition. Near the design advance ratio the difference is 
about two points in efficiency. Figure 44 shows the effect of applying the 
interpolation function to the SR-1, SR-2 , and SR-3 propeller designs at design 
point. For this study, two different exponent values were used for both the 
NASA and Mellor cascade correlations. For comparative purposes, the flat 
plate cascade correction to the isolated airfoil characteristics was also 
included. As was the case for the previous studies, the nacelle influence was 
neglected. Comparing the results obtained using the various transition region 
models holding all other variables fixed indicates that the differences 
between the models are small, a maximum of about 2 points in efficiency. 
Further comparison of only the interpolation model results for the two 
different exponents in the same manner shows an even smaller difference on the 
predicted efficiencies. The effect of replacing the cascade correction to 
isolated airfoil data with the interpolation function in the transition region 
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varies for the different propeller designs. The effect on the SR-1 is 
negligible, the efficiency is increased for the SR-2 using the interpolation 
model, and the efficiency is reduced or does not change for the SR-3 when the 
interpolation model is used. To demonstrate the sensitivity of the interpola- 
tion function model to the exponent used in the model the exponent was varied 
from 1.0 to 2.0. This variation resulted in an inconsistent behavior with 
increasing exponent for the different propeller designs as shown in figure 44. 
Generally, it appears that the use of the interpolation model in place of the 
correction model will increase the predicted efficiencies a small amount. 


Effect of Viscous Nacelle Induced Inflow 

An initial application of the analysis was performed using the SR-3 
propeller/nacelle configuration to provide a performance map of propeller 
efficiency versus propeller rotational speed for selected blade angles at a 
0.8 Mach number condition. The predictions were obtained using both the 
inviscid and viscous nacelle induced influence in the solution procedure for 
three different blade angles at the 8440 rpm condition. The differences 
between the results obtained using the inviscid or viscous nacelle induced 
inflow are so small for the conditions studied that they are not worth 
presenting. As a result, the need to include the viscous nacelle induced 
inflow solution in the prediction of high speed propeller designs such as 
the SR-3 at design flight conditions does not appear to be warranted and in 
general it is not used in the studies to follow. 


Effect of Centrifugal Blade Twist 

After these SR-3 cases had been run, the SR-3 propeller design was found 
to have a measurable steady-state torsional deflection due to centrifugal 
effects. The measured deflections were for the 8440 rpm condition and are 
presented in figure 45 in the form of a comparison between the aerodynamic 
blade angle distribution with and without the centrifugal effect. As a result 
of this finding, the cases were rerun with the measured torsional deflections 
included in the analysis. In addition, the other rpm conditions were run 
including the measured torsional deflections scaled by the square of the ratio 
of the actual rpm to the reference 8440 rpm. The results of this application 
are presented in figure 46. At the time of the performance of this activity, 
most of the effort had been directed towards the SR-3 configuration, but not 
all combinations of modeling features had been studied. And as a result, the 
combination of modeling features used to obtain the results presented in 
figure 46 are not necessarily the most appropriate. The actual combination 
used for this study was the AIR24 (NACA data) without cascade data and a tip 
loss model, with the compressible wake model. From the results presented in 
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figure 46 it is evident that the effect of including torsional deflections due 
to dynamic or aerodynamic effects can significantly alter the performance 
predictions. The inclusion of the measured torsional deflections of over one 
degree in the working sections of the blade for the SR-3 propeller resulted in 
an increase in the predicted efficiencies of 2.5 to 3.0 points. For all 
studies presented hereafter, this torsional deflection is included for the SR- 
3 configuration. For the other configurations the effect is not known, but it 
is assumed to be small due to the severe swept angles. 


Comparison With Data 


In the previous sections the results of applying the analysis with 
different aerodynamic modeling features was presented. From these studies a 
selected combination of modeling features were chosen for further investiga- 
tion and comparison with test data. The unpublished test data was provided by 
NASA Lewis Research Center and was obtained in the NASA Lewis Research 
Center’s test facilities. Initially the analysis was run including the 
nacelle’s induced influence from both the inviscid and viscous solutions in 
the propeller performance portion of the ‘ analysis . These early cases were 
run on the NASA Lewis Research Center’s Univac computer system. 


As noted earlier, the results from these cases supported two major 
conclusions concerning the use of the analysis made during the early applica- 
tion and code development work. First, the use of either the inviscid or 
viscous flow solution resulted in essentially identical propeller performance 
solutions for the high speed conditions investigated. This result is not 
surprising in view of the assumption made in the nacelle flow field solution; 
that the flow field away from the body surface can be characterized by the 
inviscid compressible flow solution. Since a very large fraction of the 
propeller blade is operating in this region, the small differences in the 
local flow field near the nacelle surface due to viscous effects have very 
little effect in the propeller performance solution. Second, the need to 
cycle between the nacelle and propeller solution portions of the analysis is 
unnecessary for the high speed conditions studied. Again, this result is not 
surprising since the effect of the propeller for the conditions being studied 
is to induce an incremental increase in the flow around the nacelle body (on a 
momentum basis) of less than four percent relative to the freestream velocity. 
It should be noted that this method of using the PANPER analysis may not be 
applicable to low speed conditions, and flight conditions where the local flow 
is nearly choked or nearly sonic since the inclusion of the blade forces in 
the nacelle solution may be important. 
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Based upon the results of these cases the analysis was rerun at UTRC 
using the nacelle induced inflow as prescribed input to the propeller portion 
of the analysis. In this manner the effects of the different aerodynamic 
modeling features were studied without the computational expense of running 
the nacelle portion of the analysis. Thus a large matrix of results was 
obtained for the different modeling features, flight conditions, and propeller 
types for comparison with test data. The propeller configurations studied 
were those noted earlier along with a variation in the spinner shape, conic 
verses area ruled. The flight conditions studied consisted of three Mach 
number variations at a fixed advance ratio and power coefficient level and 
three advance ratio variations at a fixed Mach number and power level. A 
variation in the source of the isolated airfoil data was made, the Manoni 
(AIR23) and the NACA data (AIR24) . And finally eight different combinations 
of aerodynamic modeling features were included in this study. The effect of 
using cascade data (NASA data with flat plate corrections in the transition 
region) versus isolated data alone was investigated. The use of the 
compressible wake model (Mach cone region of influence model) was compared 
with the conventional incompressible wake model. And the effect of the tip 
loss model (Evvard) was studied. In total this comparison with test results 
included 384 separate cases to establish the sensitivity of the performance 
predictions to selected modeling features. The results of this comparison are 
presented in the following sections and are also provided in Tables I, II, and 
III. 


Variation with Mach Number 

The effect of varying flight speed (Mach number) at a fixed advance 
ratio and fixed power coefficient on propeller efficiency are presented in 
figures 47 to 53 for the different isolated airfoil data sources, modeling 
features, and propeller types. These results are also presented in tabular 
form in Table I. For all of these conditions studied the area ruled spinner 
geometries were used to calculate the nacelle influence. The three propeller 
geometries used were the SR-1M, SR-2, and the SR-3 designs. Thus the results 
presented in these figures also represent the effect of different blade sweep 
along with twist, camber, and chord distribution variations. Figure 47 
presents the unpublished NASA test results for the three different 
propeller/nacelle configurations studied. These test results indicate that 
the swept blade designs (SR-1M and SR-3) are more efficient than the unswept 
design (SR-2). However, the moderately swept design (SR-1M) is more efficient 
at the lower Mach number conditions than the highly swept design (SR-3). As 
the Mach number increases, the highly swept design's performance increases 
until it becomes the most efficient design at the highest Mach number 
studied. 


c 
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One of the effects of increasing the Mach number is to increase the size 
of the supersonic region on the propeller blades. To determine the sensi- 
tivity of the performance predictions to the aerodynamic modeling features 
used to model the compressible flow effects, several of these features were 
varied in this study. The compressible wake model incorporated in the 
analysis to account for compressibility effects on the wake induced flow field 
is studied along with the tip loss models in this comparison. The effect of 
different sources of airfoil data (some airfoil) is also studied in conjunc- 
tion with the influence of cascade data. The comparisons of the different 
prediction models with the test results are presented below proceeding from 
the unswept to the highly swept design. 

SR-2 Propeller Design 

In figure 48, the effect of the tip loss and compressible wake model are 
shown without the use of the cascade data on the unswept SR-2 propeller for 
both sets of isolated airfoil data. The most obvious difference is that the 
use of the NACA airfoil data (AIR24) for the 16 series airfoil characteristics 
consistently results in a higher predicted efficiency. At the higher Mach 
number (0.8) where the effect of the compressible wake and tip loss models are 
significant, the effect of the use of tip loss model is seen to reduce the 
propeller efficiency. The effect of the use of the compressible wake model is 
to increase the propeller efficiency. In figure 49 the same results are 
presented using the cascade data. The effects of the different aerodynamic 
modeling features are identical to those in figure 48. However, the effect of 
the different airfoil data sources on the efficiency is relatively small and 
does not show a consistent trend with Mach number. This implies that the 
major differences in the airfoil data sources are at the higher chord-to-gap 
ratios, thickness ratios, and or the lower Mach numbers associated with the 
inboard sections of the blades of the SR-2 propeller. A careful comparison of 
the results in figures 48 and 49 obtained for the use of the cascade data 
versus not using the data shows that the cases using the NACA airfoil data 
(AIR24) is most affected by the use of the cascade data. The use of the 
Manoni Airfoil data (AIR23) is only slightly affected by replacing the inboard 
airfoil data with the cascade data. From the results presented in both figures 
for the SR-2 design, it is apparent that the best overall prediction of the 
character of the test results with Mach number are generally obtained in the 
cases which use the compressible wake model. The efficiency predictions, 
based on the compressible wake model, with or without tip loss, are within 2 
percent of the test data. 

SR-1M Propeller Design 

The same format is used for presentation of the comparisons of predicted 
results with data for the moderately swept propeller design (SR-1M) as was 
used for the unswept design. The comparisons are presented in figures 50 and 
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51 for the non-cascade and cascade data cases respectively. The results are 
essentially identical to the SR-2 results with the following exceptions. The 
use of the different isolated airfoil data set results in a larger difference 
in predicted efficiency. The AIR24 data again yields the higher efficiencies 
for the non-cascade data cases (figure 50). For the cases which included the 
cascade data (figure 51), sensitivity to airfoil data is again reduced. 
Generally, similar results are achieved with the various modeling options at 
the 0.6 and 0.7 Mach numbers. Significant differences are evident at the 0.8 
Mach number. 

SR-3 Propeller Design 

Similar results are presented in figures 52 and 53 for the highly swept 
SR-3 blade design. Comparing predicted versus test values of propeller 
efficiency for the non-cascade data cases in figure 52 shows that the same 
behavior with airfoil data source as was noted in the other propeller designs 
has occurred. The difference in efficiencies between the predicted values is 
even greater than for the SR-1M cases. However, the difference in the 
predicted results using the different airfoil data sources for the cascade 
data cases is again small. Inclusion of the compressible wake model again 
increases the efficiency predictions and the inclusion of the tip loss model 
reduces the predicted efficiencies for both the cascade and non-cascade 
results. For this propeller configuration the character of the measured test 
results is well predicted using most of the models investigated. However, the 
simultaneous use of both the compressible wake and the tip loss models with 
AIR24 airfoil data results in a large exaggeration of the predicted loss of 
efficiency with increasing Mach number compared with the behavior exhibited in 
the test results for this SR-3 propeller design. The compressible wake model 
without tip loss results in good correlation for the SR-3 propeller. 
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Variation With Advance Ratio 


The effect of varying the propeller advance ratio for a fixed power 
output at a fixed Mach number of 0.8 on the measured efficiency is shown in 
figures 54 to 60 for three propeller/nacelle configurations. For this study 
the area ruled spinners were used for the three different propeller designs; 
the SR-1M, SR-2, and SR-3. The power output level studied corresponded to the 
design value at the 3.07 advance ratio condition. As illustrated by the test 
results presented in figure 54, the effect of increasing the advance ratio at 
a fixed power level is to reduce the propeller efficiency. For the fixed 
free-stream Mach number, increasing the advance ratio corresponds to 
decreasing the propeller rotational speed, and when done at a fixed power 
level the increase in advance ratio also corresponds to an increase in the 
propeller blade pitch angle. This increase in blade angle reduces the lift- 
to-drag ratio that the blade sections operate at and the result is reduced 
propeller efficiency. At this Mach number (0.8), the effect of increasing 
blade sweep is to increase the efficiency. This behavior is consistent over 
the range of advance ratio studied, but at the higher advance ratios the 
difference between the highly swept (SR-3) and the moderately swept (SR-1M) 
designs vanishes because the SR-3 results exhibit a larger negative gradient 
with advance ratio. The prediction capability of the analysis is demonstrated 
in figures 55 to 60 for this study, where the predicted propeller efficiencies 
are compared with the measured results using the same format as was used in 
the previous sections. The predicted efficiencies for this study are also 
presented in tabular form in Table II. 

SR-2 Propeller Design 

In figure 55 the results of applying the analysis to the unswept 
propeller design are presented for the cases run without cascade data. The 
effect of the use of the different isolated airfoil data sets is seen to be 
similar to the behavior noted in the studies performed for the variation with 
Mach number. That is, the use of the AIR24 data results in consistently 
higher predicted efficiencies over the advance ratio range investigated. 

Again, the inclusion of the compressible wake model results in higher 
efficiencies and the inclusion of the tip loss model results in lower 
efficiencies. In figure 56, the same comparisons are presented including the 
use of the cascade data. Use of the cascade data again reduces the differ- 
ences in the predicted efficiencies between results for the two isolated 
airfoil data sets to relatively small values. In general, all of the results 
from the prediction models studied follow the trends of the test results, but 
the overall levels are not always well predicted. The general character of 
the test results seems to be well predicted in most cases. The best correla- 
tion with test results is obtained with the compressible wake model. The 
predictions using the combination of the compressible wake and tip loss models 
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with or without cascade data show an unrealistically strong sensitivity to the 
compressibility effects associated with increasing tip Mach number (lowering 
advance ratio) compared with the behavior of the test results. 

SR-lM Propeller Design 

In figures 57 and 58, the results of applying the modeling features and 
airfoil characteristics to the moderately swept blade design (SR-lM) are 
presented with and without the use of cascade data respectively. The 
predicted results exhibit essentially the same characteristics for the 
different modeling features as exhibited by the unswept propeller results 
described above. The general character of the test results is reproduced 
using the prediction analysis with all of the different models, however, the 
gradient of the curves with advance ratio is usually larger than the test 
results indicate. Generally better correlation was again obtained with the 
compressible wake model. 

SR-3 Propeller Design 

Application of the same set of modeling features to the highly swept SR-3 
propeller design is presented in figures 59 and 60. Again, the overall 
behavior of the predicted results with the different modeling features is 
consistent with the results from the other propeller designs. The most 
obvious difference between these predicted results and the results obtained 
for the SR-2 and SR-lM propellers is the larger negative gradients with 
advance ratio. The increased gradient is consistent with the test results 
which exhibit the same behavior (figure 54), but the predicted slopes are 
significantly greater than the test results indicate. This behavior is 
exhibited for all of the combinations of the modeling features. 


Variation With Blade Twist and Spinner Shape 

The effect of spinner design and blade twist on the measured propeller 
efficiency as measured is presented in figure 61. The results presented in 
this figure are for a propeller blade twist variation (SR-lM versus SR-1) 
using the conic spinner geometry, and a spinner design variation (conic versus 
area ruled) for the SR-lM propeller. The configurations studied were run at a 
fixed Mach number (0.8) and a fixed power output level equal to the design 
power as was the case for the previous study. From the test results (figure 
61) it can be seen that the effect of decreasing the blade twist (SR-lM in 
place of the SR-1) along with increasing the blade camber on the tip region of 
the propeller blade is to improve the efficiency over most of the advance 
ratio range investigated. The effect of using the area ruled spinner in place 
of the conical geometry also results in an increase in the measured propeller 
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efficiency over the complete advance ratio range. The results of the applica- 
tion of the various modeling features of the analysis are shown in figures 57, 
58, and 62 to 65 in the same format as was used for the previous investiga- 
tion. The tabulation of the predicted efficiencies is provided in Table III. 

SR-1M Propeller With Area Ruled Spinner 

The results for this propeller design have been presented in the previous 
study (figures 57 and 58). 

SR-1M Propeller With Conic Spinner 

The results of the use of the analysis to predict the performance of 
the SR-1M propeller with the conic spinner for the different modeling feature 
is presented in figures 62 and 63. In figure 62 the efficiencies using the 
different modeling features are presented without using the cascade data, 
while in figure 63 the results for the same combination of modeling features 
are presented with the cascade data. Comparing the results presented in 
figures 62 and 63 with those presented in figures 57 and 58, it is noted that 
the increase in efficiency of the area ruled spinner over the conic spinner, 
indicated by the test data, is not predicted by the PANPER analysis. Relative 
to all preceding theory-test comparisons, the conic spinner predictions 
exhibit the poorest correlation with test data. 

SR-1 Propeller With Conic Spinner 

Results are shown in figures 64 and 65 for the SR-1 propeller design and 
the conic spinner. With this configuration there is an improved prediction of 
the gradient of efficiency with advance ratio. Comparing the results 
presented in figures 64 and 65 with those presented in figures 62 and 63, it 
is noted that increased efficiency with increased blade twist (SR-1M vs. SR-1) 
is predicted. This is the general trend of the test data (figure 61), 
although the magnitude of the difference is overpredicted, especially at the 
lower advance ratio where the measured efficiencies merge. 


Correlation of Efficiency at Prescribed Power Levels 

The preceding comparisons of predicted propeller efficiency with test 
data at prescribed power coefficients or power loadings, presented in figures 
48 through 65, indicate that the predictions are sensitive to the choice of 
aerodynamic modeling features and airfoil data. The test data trends of 
decreasing efficiency with increasing free-stream Mach number and increasing 
advance ratio are generally predicted by most of the modeling options. The 
general influence of twist is overpredicted, and the predicted influence of 
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sweep and nacelle design (conic vs. area ruled) is not generally consistent 
with test. Sensitivity to airfoil data selection is apparent without the use 
of cascade data, but the differences with airfoil data are small when cascade 
data is used. Overall, it is difficult to select a "best" modeling option, 
because there is not one that demonstrates consistent correlation throughout. 
The use of the compressible wake option generally results in improved 
correlation. The results of the use of tip loss are inconsistent between 
propeller designs. The modeling option combination of compressible wake with 
no tip loss was selected as the primary modeling option for the forthcoming 
performance map predictions. With the notable exception of the SR-1M, conic 
nacelle predictions, the predicted apparent efficiencies with this modeling 
option are generally within approximately ± 2 percent of the test data when an 
average of the A1R23 - AIR24 predictions is used. Considering the limitations 
of lifting line theory, as applied to the severe propfan design features and 
operating conditions, the sensitivity to the airfoil data and modeling 
options, and the accuracy of the test data, the degree of correlation is not 
unreasonable. 


Performance Map Prediction and Comparison with Test Data 

A final comparison with data is presented in this section for the single 
propeller mode of operation. In this section predicted performance maps are 
compared with measured performance maps for the SR-1M, SR-2 , and SR-3 . 
propeller configurations at the 0.6 and 0.8 Mach number conditions of the 
previous sections. The unpublished test data was supplied by NASA Lewis in 
graphical form for the three configurations at these Mach numbers. The maps 
are plots of propeller efficiency (net) and power coefficient as functions of 
advance ratio for several fixed blade angles. For this study, the analysis 
was run for three different comparisons using different combinations of 
aerodynamic modeling features and several blade angles to illustrate the range 
of correlation that can be obtained. For modeling option 1, the aerodynamic 
modeling features used consisted of the AIR24 data without cascade data, no 
compressible wake, and no tip loss models. Modeling option 1 was applied to 
the SR-2 and the SR-1M configurations. For modeling option 2 the compressible 
wake was used without the tip loss model, and the AIR24 data was used in 
conjunction with the cascade data and transition region interpolation feature. 
For modeling options 1 and 2, the PANPER analysis was run at measured blade 
angles. To illustrate the trending capabilities of the analysis, modeling 
option 3 was applied using the same combination of modeling features as 
modeling option 2, but the blade angles were adjusted by a constant value so 
that the measured and predicted power coefficients matched at one Cp, J point 
for each Mach number condition for each configuration. The point chosen was 
generally near the design point condition. Modeling options 2 and 3 were 
applied to all three propeller conf iguations . Both the test and predicted 
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results are presented in figures 66 to 71 for the SR-2 , SR-1M, SR-3 respec- 
tively. The predicted results are also provided in Tables IV to VI for the 
same configurations respectively. These results are described for each 
configuration in the following sections. 

^R^2_Res_ulJts 

The performance predictions for the unswept SR-2 propeller, using the 
three modeling options, are presented in figures 66 and 67. At the 0.6 Mach 
number (figure 66a), the predicted efficiencies for modeling option 1 are 
higher than the test values and the peaks occur at greater advance ratio for 
each blade angle studied. The predicted power coefficients are also higher 
than the test values for the same blade angle. The difference in power 
coefficient is consistent with the fact that the analysis has been unable to 
correlate with measured efficiency and power at the measured blade angle, 
while efficiency correlations have been much better when trimmed in blade 
angle to match Cp. The comparison obtained with modeling option 2 for this 
Mach number (figure 66b) is very similar. Peak efficiencies are not well 
predicted in either amplitude or advance ratio, while the power coefficient 
levels are too high. However, if the blade angle is adjusted to match advance 
ratio and power level for one point on the test data (modeling option 3), the 
results are greatly improved (figure 66c). Peak efficiencies are better 
predicted in both absolute value and in terms of the advance ratio (at which 
they occur) as compared with the other modeling options. At the 0.8 Mach 
number (figure 67a), the predicted efficiencies for modeling option 1 were 
closer to the test values than the 0.6 Mach number results, but the peak 
values were still predicted to occur at greater advance ratio than the test 
results indicate. For modeling option 2 (figure 67b), the comparison is 
slightly degraded. The results obtained using modeling option 3 (figure 67c) 
are improved in terms of the predicted power coefficient behavior as a 
function of blade angle and advance ratio. But, the predicted propeller 
efficiencies are now higher than the test results and the predictions obtained 
using modeling option 1. 

SR-lM jtes_u_lts_ 

The performance predictions for the SR-1M propeller, using the three 
modeling options, are presented in figures 68 and 69. For the 0.6 Mach number 
condition (figure 68a), the degree of correlation is similar to the SR-2 
correlation using modeling option 1. No significant improvements were 
obtained at this Mach number for modeling option 2 (figure 68b). When 
modeling option 3 was used, the correlation (figure 68c) was again improved, 
efficiencies and power levels were better predicted. For the 0.8 Mach number 
comparison with modeling option 1 (figure 69a), the SR-1M predicted results 
were higher than the test results for both the efficiency and power 
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coefficient comparisons. The peak efficiencies were also predicted to occur 
at higher advance ratio than the test results. Again, no significant improve- 
ments were noted when modeling option 2 (figure 69b) was used. The comparison 
based on modeling option 3 (figure 69c) indicated some improvement in power 
coefficient correlation, but little improvement in predicted efficiency. 

SR-3_RjJ sji l_t^s 

The performance predictions for the SR-3 propeller are presented in 
figures 70 and 71. At the 0.6 Mach number condition (figure 70a), the 
correlation using modeling option 2, for both the efficiency and power 
coefficient, is very good relative to the other configurations. The results 
obtained for modeling option 3 are slightly different than for the other 
options because the trim point is slightly different (figure 70b). At the 0.8 
Mach number condition (figure 71a), both the power coefficient and the 
efficiency are over-predicted for modeling option 2, although the peak 
efficiency predictions do not occur at higher advance ratio as was the case 
for the other two configurations. The correlations obtained for modeling 
option 3 (figure 71b) for this Mach number are not as good as those for the 
0.6 Mach number condition. Overall the SR-3 correlation is the best of the 
three configurations studied. 

Performance Map Discussion 

These performance map comparisons illustrate that the analysis can be 
tuned to provide useful trending results if a reference point is known in 
advance (modeling option 3). The correlations for the SR-2 and SR-lM could be 
improved by selecting different sets of aerodynamic modeling options which are 
best suited for their particular geometries and adjusting the blade angle to 
match power coefficient at one point. However, without this type of informa- 
tion in advance, the predictions are sensitive to the selection of modeling 
options and the results can vary significantly. 
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WIND TUNNEL APPLICATION 


A major task supported by the contract activity was to apply the analysis 
to several wind tunnels and propeller/nacelle configurations to study the 
effect of the tunnel-propeller/nacelle geometry on the predicted performance 
of the propeller. The tunnels selected for the study were the U.S. Tullahoma 
tunnel and the French Modane tunnel. Details concerning these tunnel 
facilities can be found in references 70 and 71 for the Tullahoma and Modane 
tunnels respectively. The key geometry feature used in the modeling of the 
tunnel with an axisymmetric code is the wall diameter. The tunnel wall 
diameter for the Tullahoma tunnel is 5.6 meters (18.054 ft) and the Modane 
tunnel diameter is 8.0 meters (26 ft). Several SR-3 propeller/nacelle config- 
urations were studied with propeller diameters of 3.08, 2.76 and 2.29 meters 
(11, 9, and 7.465 ft) and nacelle diameters of 1.27, 1.04, and 0.85 meters 
respectively (4.122, 3.373, and 2.754 ft). The resulting propeller to nacelle 
diameter ratios are essentially constant so that the geometric variations 
studied represent the effect of propeller to tunnel wall diameter ratios. The 
operating conditions studied consisted of variations in Mach number and power 
coefficient. Since the Tullahoma tunnel can be operated with tunnel wall 
bleed, most of the conditions were run for this tunnel with the wall bleed 
model. Because the actual bleed rate is unknown, selected values were used in 
the predictions to study the influence on the propeller performance. The 
cases studied are summarized in Table VII, which includes the predicted 
propeller efficiencies. The application of the analysis to these conditions 
was made during a period of the contract activity when the most appropriate 
combination of modeling features was felt to be the AIR24 isolated airfoil 
data used without any cascade data, with the compressible wake model, and 
without the tip loss model. 


Effect of Tunnel/Propeller Geometry 

The matrix of conditions and configurations studied (Table VII) is not 
complete in terms of all of the possible variations that could be studied. 
However, there is sufficient information to identify predicted trends due to 
the effect of the nacelle to tunnel diameter ratio on the performance of the 
propeller. In figure 72 the effect of this geometric ratio is illustrated, 
where the apparent efficiency is plotted against the ratio of propeller to 
tunnel diameter for the 0.8 Mach number conditions. In addition, the nacelle 
to tunnel wall diameter ratio is presented on the plot. When the ratio is 
zero, the condition corresponds to free air (no tunnel) and when the ratio is 
equal to unity, the propeller would completely fill the tunnel. The results 
presented in this figure are for the non-bleed configurations at the various 
power coefficients levels studied. Within the limitations in terms of the 
range of the data, the effect of increasing the propeller diameter relative to 
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the tunnel diameter is shown to be relatively insensitive up to about the 
ratio of about 0,3 depending on the loading condition. Beyond this region the 
effect of increasing the ratio is to reduce the predicted efficiencies. A 
case was run for a value of this ratio of 0.61 and the nacelle portion of the 
code failed to run completely through the solution procedure. The failure was 
identified as being caused by the predicted choking of the flow in the sense 
that a mass flow balance could not be maintained through the tunnel with this 
particular nacelle diameter size. This numerically predicted result indicated 
that the use of this configuration would be inappropriate to consider for 
testing in the actual tunnel at the 0.8 Mach number. The effect of power 
loading in relation to the nacelle to tunnel diameter ratio variation is also 
shown in this figure. The predicted behavior of the propeller efficiency with 
decreasing power loading is shown to increase regardless of the nacelle to 
tunnel ratio. This result is consistent with the results demonstrated in 
figure 71 at the corresponding advance ratio (2.9). There is an indication 
that the effect of the nacelle to tunnel ratio is greater for the higher 
loading conditions, although the data is too limited to make a strong 
statement. The influence of tunnel size on the predicted propeller perfor- 
mance is shown in figure 73 for a fixed propeller/nacelle geometry, power 
coefficient, and propeller tip speed as a function of Mach number, again for 
the non-bleed conditions. The efficiency is seen to increase with increasing 
Mach number and increase with increasing tunnel diameter. The free air result 
shown is slightly lower in efficiency than the largest tunnel result, but the 
difference may be attributed to the exactness of the modeling of the free air 
condition with the analysis or to the fact that the bleed model as incor- 
porated in the analysis should include an iteration between the calculation of 
the local bleed rate and the local pressure field solution. 


Effect of Wall Bleed 

The analysis currently includes the treatment of wall bleed in the 
viscous nacelle flow solution portion of the computer program. For thi^^study, 
the bleed rate is defined as the total mass loss or gained in the tunnel 
(positive bleed or negative bleed, respectively) through the tunnel wall 
divided by the initial mass flow. Since the bleed model is only applied over 
the region of the tunnel wall modeled in the viscous flow solution procedure, 
this definition for the bleed rate is non-standard. This region starts from 
the axial location on the tunnel wall which corresponds to the start of the 
nacelle and extends to the end of the computational region. Thus a large 
fraction of the complete tunnel wall is neglected. The analysis was run using 
various bleed rate values by prescribing the plenum stagnation pressure and 
temperature and the hole size and discharge coefficient. The actual bleed 
rate varies along the wall surface and is a function of the solution, thus it 
is not possible to define the bleed rate a priori. The bleed rate is 
calculated once the inviscid solution is known to define the value for use in 
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the viscous solution. Thus the values used for this study are not necessarily 
realistic values that would be used in an actual test. The intention of this 
study is to show the sensitivity and capability of the analysis. The results 
presented in Table VII demonstrate the effect of the use of the wall bleed 
model on the predicted propeller efficiencies for a fixed propeller/nacelle 
and tunnel wall configuration. The variation in bleed rate was from - 2.0 to 
5.0 percent. The variation in the predicted efficiencies with bleed rate is 
seen to increase with the inclusion of either positive or negative bleed, with 
the exception of one condition which had a larger tunnel to nacelle diameter 
ratio. The results are plotted in figure 74 for the 0.8 Mach number with the 
nacelle tunnel diameter ratio of 0.499. An interesting result shown in this 
figure is that the small (essentially zero) bleed rate conditions result in an 
increase in the predicted efficiency as compared with the no bleed model 
values shown on this figure (solid wall). It may be that the use of the wall 
bleed model with such low rates, either positive or negative, allows the 
solution to be locally unconstrained compared with the conventional zero 
normal flow tunnel wall boundary conditions. In fact, for the lower negative 
bleed rates nearly zero (-0.1 percent net loss) the local bleed is positive 
(mass gained) for some distance along the wall upstream of the propeller blade 
row. Thus the perforated tunnel wall acts like a partial free air boundary 
for low bleed rates. This result is not really surprising in view of the fact 
that one of the functions of the perforated tunnel wall is to alleviate the 
tunnel wall influence on the tunnel flow. The free air results for two of the 
loading conditions are also shown on this figure. These results illustrate 
that the use of the bleed model in the analysis does tend to result in 
predicted propeller efficiencies nearly equal to the free air results. 

In general, the use of the feature to model the influence of tunnel 
walls, with or without wall bleed, shows definable trends with the parameters 
of interest. This feature could be used to study the potential influence of 
particular model to tunnel size configurations to determine the merits of the 
selected configurations and the interference effects. With this capability it 
is possible to numerically define the relative differences between different 
tunnel facilities for comparing the test results obtained in the different 
tunnels using the same models. 
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COAXIAL PROPELLER PERFORMANCE PREDICTIONS 


The analysis is capable of treating counter-rotating propeller/nacelle 
configurations. In order to demonstrate this capability and to provide a 
preliminary assessment of the predictive accuracy of the analysis in this mode 
of operation, the analysis was applied to three propeller/nacelle configura- 
tions. The flight speeds consisted of three moderate to high Mach number 
conditions (0.6, 0.7, and 0.8) for two of the propeller/nacelle geometries and 
one low Mach number condition (0.04) for the remaining propeller/nacelle 
geometry. The conditions and geometries studied are defined in the following 
sections . 


Counter-Rotation Application - Low Speed 

The analysis was applied to the propeller/nacelle geometry described in 
reference 72 for three combinations of blade angle. For each blade angle 
setting, the propeller advance ratio was varied to provide a performance map 
of efficiencies and power coefficients. Initially, the AIR23 airfoil data 
base was used because of the limits on the range of design lift coefficient in 
the AIR24 airfoil data base. The propeller design lift coefficients varied 
from 0.1 to 0.3 while the AIR24 data was limited from 0.0 to 0.2. The 
predicted performance results of this application are presented in Table VIII 
for each condition. The performance predictions were obtained using the 
combination of modeling features deemed best for the single propeller applica- 
tion. No tip loss or compressible wake models were used, and because the 
propeller design was a 3x2 (three blades per blade row) the use of cascade 
data was clearly unnecessary. The nacelle and tunnel wall influence was 
calculated from a single propeller/nacelle case at the lowest Mach number for 
which the nacelle portion of the analysis would run (0.2). The resulting 
nacelle induced influence was used as input to the propeller portion of the 
analysis for all of the conditions which were run. The predictions are 
compared in figures 75 to 77 with the measured results obtained from figure 8 
and figure 12 of the above reference. 

In figure 75, figure 8 of reference 72, a comparison of the predicted and 
measured power coefficients versus advance ratio are presented for three blade 
angles (55, 60, and 65 degrees). Included on the predicted curves are the 
predicted efficiencies at the discrete advance ratio points for which the 
analysis was run. The power coefficients tend to be over-predicted at the 
lower advance ratios and under-predicted at the higher advance ratios when 
compared with the measured values. The predicted efficiencies tend to be 
lower than the measured values. Figure 76 presents a comparison of the 
propeller thrust coefficient as a function of advance ratio for the three 
blade angles. The correlation between predicted and measured thrust results 
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is better than the power coefficient results, although the difference between 
measured and predicted values increases with increasing blade angle and the 
peak values are progressively under-predicted with increasinf blade angle. 

Thus in general, the differences between the measured and predicted results 
tend to become larger with increasing blade angle. These differences could be 
attributable to two problems. First, the propeller blades were made of wood 
and thus the blade angles and twist distributions used in the test are of 
questionable accuracy. This problem was noted in the reference (pages 7 and 
9). The question of retwisting the blades due to centrifugal and aerodynamic 
loading of the thin wood blades may also be important. The comparisons 
exhibit the trend of increasing difference as the blade angle and propeller 
rpm increase (increased blade airloading) which is symptomatic of the blade 
twisting with increased loading. As shown in an earlier section of this 
report, the effect of centrifugal loading on the blade twist is important and 
the sensitivity of the performance predictions to blade angle is quite strong. 
Second, the airfoil data base may not adequately represent the actual airfoil 
characteristics. For the higher power coefficients the range in angle of 
attack of the airfoil data tables was exceeded introducing an artificial stall 
like behavior into the solution as noted in Table VIII. This probably 
explains why the peak thrust coefficients are under-predicted. 

In figure 77, a comparison of the section blade thrust distributions are 
shown at two blade angles (55 and 65 degrees). The overall trending is well 
predicted. The largest differences occur outboard on the propeller for the 55 
degree blade angle and inboard for the 65 degree blade angle. Again, this may 
be due to twisting of the blade with aerodynamic load increasing as the blade 
angle is increased. 

To further illustrate the sensivity of the solution process to small 
variations in blade twist and airfoil section properties associated with 
testing uncertainies, the analysis was rerun using the AIR24 data base. With 
this data base, if the section design lift coefficient range (0.0 to 0.2) was 
exceeded, the lift and drag coefficient values at the bounds (0.0 or 0.2) were 
used. One blade angle (60.0 degrees) was studied using only the change in 
airfoil data base and these results are also included in figures 75 and 76. 

The predicted efficiency was increased significantly with a slight drop in 
predicted power coefficient in the range of advance ratio for which the range 
in angle of attack in the airfoil tables was not exceeded. The results of 
increasing the blade angle by one-half degree using the AIR24 airfoil data 
base are also included in figures 75 and 76. These results show the need to 
have very accurate blade angle measurements for the test conditions; and by 
inference, an accurate measure of the twist distribution. 
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Counter-Rotation Application - Moderate to High Speeds 


The analysis was applied to two propeller/nacelle configurations 
consisting of different spinner designs, A and B, described in reference 73 
and shown graphically in figure 78. For the A spinner, two blade angles (60.0 
and 65.0 degrees) were run for three Mach number variations (0.6,* 0.7, and 
0.8). For each of these conditions an advance ratio variation was made to 
provide performance maps of predicted efficiencies and power coefficients. 

For the B spinner only one Mach number-blade angle condition (0.7,65.0 deg.) 
was investigated. The predicted performance results for these conditions are 
presented in Table IX. The aerodynamic modeling options used for this study 
were the compressible wake model with the tip loss model. The influence of 
the nacelle and tunnel wall was included in the propeller portion of the 
analysis in the same manner as noted above, by applying the nacelle portion of 
the analysis once to each Mach number condition. The spinner designs resulted 
in inflow profiles which were essentially uniform for spinner A, while the 
spinner B inflow profiles exhibited deceleration and acceleration of the flow 
at the roots of the front and rear propeller blades respectively. This 
comparison is shown in figure 79. Again, the AIR23 airfoil tables were used 
because the blade design lift coefficients were off the tabulated range in the 
data base. The performance predictions are compared in figures 80 to 83 with 
the measured test results obtained from figures 8d, 8e, 8g, 9b, 10b, 14d, 14e, 
14g, 15b, and 16b of the reference noted above. 

The test results were obtained with different front and rear blade 
angles. The values noted above are the front propeller values. Two different 
cases are presented; one with a differential blade angle between the front and 
rear blade angles which varied with pitch setting and which corresponded to 
equal power loading on both propellers (assummed optimum) near the advance 
ratio corresponding to peak efficiency, and one with a fixed value of 0.8 
degrees. The details of the values used and the propeller designs are 
available in reference 73. The optimum differential blade angles were blade 
angle dependent. For the 60.0 degree blade angle condition (front) the value 
was 2.0 degrees, and for the 65.0 degree blade angle condition the value was 
approximately 2.5 degrees. 

In figure 80, the effect of Mach number and blade angle is illustrated 
for the optimum differential blade angle. For this comparison four blades 
were used on each propeller. The predicted power and thrust levels are 
greater than the measured values for all blade angles and Mach numbers 
investigated. The efficiencies are under-predicted at low advance ratio and 
over-predicted at higher advance ratio for all of the conditions studied. 
Increasing Mach number is seen to decrease the predicted efficiency which is 
consistent with the trending of the measured results. A comparison of the 
variation in the power loading between front and rear propellers for these 
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conditions is shown in figure 81. Although the level is over-predicted, the 
trends are well predicted for the 0.6 and 0.7 Mach number conditions. At the 
0.8 Mach number the trending is not predicted as well, in so much as the 
advance ratio corresponding to equal power loading is not predicted 
correctly. 

The results shown in figure 82 illustrate the effect of the spinner 
design on the propeller performance for three bladed propellers operating with 
a differential blade angle of 0.8 degrees at the 0.7 Mach number. The spinner 
induced inflows, as predicted by the analysis, were shown in figure 79 to be 
significantly different. Both the measured and predicted results show 
slightly higher performance for Spinner A, which had essentially no effect on 
the flow field. However, the difference in predicted efficiency for the two 
spinners is very small. Overall, the predicted thrust and power levels for 
both spinner designs were greater than those measured. It should be noted 
that for these three bladed propeller configurations the range in angle of 
attack on the airfoil data base was exceeded for most conditions. The corres- 
ponding variation of power loading between front and rear propellers for this 
study is shown in figure 83. Again, the levels are over-predicted but the 
trends are well predicted. 

The results of this coaxial propeller comparison, based on test data from 
reference 73, show generally adequate correlation in terms of predicted versus 
measured trends, but the over all performance (power and thrust) is over- 
predicted. This discrepancy has not been adequately resolved. Some of the 
differences could be due to airfoil characteristics not being modeled 
correctly with the existing data base, and differences in the measured blade 
angles and blade twist could account for some of the differences. The fact 
that the trends are well predicted for the lower Mach numbers studied 
indicates some problems with either the airfoil data base or that the 
selection of aerodynamic modeling features used in this study are 
inappropriate. As was the case for the low speed study, some conditions 
resulted in stalled airfoil sections as noted in Table IX. As shown in the 
previous low speed study, the use of different airfoil data helped the 
correlation. The use of the AIR24 data was not investigated for this moderate 
speed study because the range in design lift coefficient for this airfoil data 
base (0.0 to 0.2) was so far removed from that associated with the test 
geometry range (0.1 to 0.5) that it would be a meaningless study. 
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CONCLUDING REMARKS AND RECOMMENDATIONS 


The Propeller/Nacelle Aerodynamic Performance Prediction Analysis 
(PANPER) was developed to provide an efficient propeller performance and pro- 
peller/nacelle airflow prediction code. This was performed by combining two 
existing methods to treat the influence of the individual components (pro- 
peller and nacelle) on each other by super- imposing their respective 
influences in the individual solution processes where appropriate within the 
constraints of the analyses. The propeller portion of the analysis predicts 
the propeller performance and the nacelle portion of the analysis predicts the 
steady axisymmetric airflow. Since the propeller analysis uses inviscid, 
incompressible, lifting-line theory with some ad-hoc features to approximate 
many of the compressible and three-dimensional aspects of the flow problem, 
there are many modeling features which can be used to tune the predictions. 

As a result, a large number of cases have been run to study the sensitivity of 
the predictions to flight conditions and combinations of these modeling 
features. These cases do not represent all possible combinations of these 
features, but they do serve to demonstrate the sensitivities of the predic- 
tions to the assumptions. The nacelle analysis was originally developed and 
successfully validated for compressible, viscous, internal flow problems. The 
modification of the original analysis to treat external flows was a signif- 
icant extension of the methodology. The comparisons shown between the pre- 
dicted and measured results represent an attempt to test the assumption of 
steady axisymmetric flow for the prediction of a time averaged unsteady flow 
and to ascertain the success of the method to this application. 

The use of the PANPER analysis for the prediction of the aerodynamic 
performance of high speed propellers (Propfans) has been demonstrated. The 
analysis was applied to a large variety of operating conditions applicable to 
the high speed propeller designs. The propeller operating conditions studied 
included windmilling to high speed conditions with supersonic tip speeds. The 
designs used for the studies included both single and coaxial configurations. 
The analysis was applied to configurations including wind tunnel wall geometry 
with and without wall bleed effects. Comparisons with data were made wherever 
possible to illustrate the degree of correlation possible which can be 
obtained using this analysis and to define the limitations of the use of the 
analysis. Based on these studies the following conclusions have been made. 


Nacelle Flow Induced Velocities 

1. The use of the analysis for propeller performance predictions 
requires that the nacelle induced flow field be accurately predicted. The 
limited comparisons with data indicate that the nacelle portion of the analy- 
sis does a reasonable job of predicting the velocity field upstream of the 
non-lifting blades. Since the general shape of the predicted inflow profile 
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is very similar to the measured profile, the differences between the test 
results and the predictions would appear to be due to inaccuracies in the 
assumed freestream conditions. It should be noted that inaccuracies in the 
assumed freestream velocity will strongly affect the correlations made at 
fixed blade angle. 

2. Inflow predictions at the blade row from the nacelle analysis have 
exhibited strong sensitivity to blade thickness modeling and the chordwise 
location at which the prediction is made (leading to trailing edge). This can 
have a significant effect on the effective section pitch distribution, and 
therefore the predicted performance. The inclusion of blade thickness effects 
in the predicted velocity profiles was shown to be necessary for good correla- 
tion with measured values upstream of the blade rows. However, the inclusion 
of blade thickness effects in the analysis to predict the nacelle induced 
velocity at the blade row introduces the question of the use of cascade air- 
foil data in the propeller analysis portion, since this data also contains the 
effect of blade thickness. 


Single Propeller Performance 

1. Based on the single propeller designs studied, the analysis is 
generally capable of predicting the high speed propeller performance to within 
two efficiency points about the design conditions when trimmed in blade angle 
to the measured power. However, the selection of aerodynamic modeling 
features and the accuracy of measured airfoil characteristics and actual blade 
pitch distributions strongly affect the correlation results as demonstrated; 
this correlation was obtained knowing the test results beforehand. This sen- 
sitivity is due in part to the transonic nature of the problem and also the 
low aspect ratios and swept geometries of the propeller blades for which the 
lifting-line concept is not rigorously applicable (although desirable from a 
computational efficiency viewpoint). 

2. Comparisons of predicted and measured results at fixed blade angle 
show significant differences in the degree of correlation. Best correlation 
is obtained near the design operating conditions. Increasing either the Mach 
number or the blade sweep (via propeller designs, SR-2 , SR-1M, and SR-3) 
degrades the degree of correlation for the fixed blade angle comparisons. 

This degradation of correlation can be attributed to some degree to the above 
noted problems. 

3. Selection of the most appropriate aerodynamic modeling features and 
airfoil data bases which yielded the best correlation with the test data 
varied with each type of propeller design, and to some extent with the type of 
comparison being made (i.e., Mach number variation, fixed blade angle, etc.). 
As a result, the application of the analysis to propeller designs for which no 
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test data exists should be made with care, and the absolute level of the 
predictions viewed with discretion. 

4. The results obtained illustrated a tendency to degrade in the level 
of correlation with increasing blade sweep at the design Mach number (0.8) as 
the advance ratio was varied. For this condition, the variation in advance 
ratio also represented variation in the tip Mach number. Since the effect of 
increasing blade sweep is to alleviate compressibility effects, this result 
may indicate that the analysis (lifting-line) has been extended in application 
beyond an acceptable limit. 

5. The analysis can be used for performance trending. However, the 
prediction of the occurence of peak efficiency as a function of blade angle 
and advance ratio is not sufficient for final design work. In addition, any 
significant change in propeller design was shown to have significant impact on 
the ability to correlate with test results for a given set of modeling 
options. This implies that trending predictions across propeller designs 
should be viewed carefully in terms of absolute levels. 


Coaxial Propeller Performance 

1. The correlation obtained at low speed using the analysis in the co- 
axial mode of operation at fixed blade angle was generally acceptable when the 
blade angle was low enough to avoid stalled operation of the blade sections. 
The actual correlations varied from over- to under-prediction of the power 
coefficient with increasing advance ratio and good correlation for the thrust 
predictions. The correlation on the propeller efficiency was within four 
efficiency points for the lower blade angle conditions. Differences between 
test and theory are believed to be partially due to inaccuracies in the blade 
angle and pitch distribution measurements and airfoil characteristics. 

2. The correlations obtained for the high speed coaxial propeller, run 
at fixed blade angle, indicated an over-prediction of both thrust and power 
levels for all advance ratios. Comparisons of test and predicted efficiencies 
show an under-prediction at lower advance ratio and an over- predict ion at 
higher advance ratio. On a front to rear propeller basis, the proportioning 
of front and rear loading is shown to trend correctly with increasing advance 
ratio. Again, the differences between test and theory may be attributable to 
all of the problems noted above, with the exception of low aspect ratio and 
sweep . 
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Wind Tunnel Effects 


The use of the analysis to analyze the effect of wind tunnel configura- 
tion was demonstrated. The results indicate that the effects of the wind 
tunnel influence can be predicted using the analysis. Decreasing the tunnel 
size vas shown to affect the propeller performance in a consistent manner and 
predict the onset of choked flow at the nacelle for too small a tunnel 
diameter. The predicted results were also shown to be propeller loading 
dependent. Although no test results were available for direct correlation, 
this aspect of the analysis may prove to be a useful tool for evaluation of 
propeller/tunnel configurations in advance of the testing procedure and to 
provide estimates of the tunnel effects on the propeller results to correct 
the performance to free air conditions. 
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Recommendations 


Based on the results obtained under this investigation and the remarks 
noted above, the following recommendations are presented. 

1. The use of the lifting-line propeller portion of the analysis should 
be limited to initial trending analysis. The sensitivity of the results to 
the aerodynamic modeling features and the questionable application of lifting- 
line theory to the highly three dimensional flow problems associated with the 
high speed propeller configurations result in predictions where the absolute 
level must be suspect. The degree of accuracy needed for current design work 
is too high for this type of analysis since efficiencies of one half of a 
point can be significant in terms of the life-cycle cost of a new design. 

2. When using the analysis for trending applications, the results 
obtained should be viewed with care at the higher Mach number and tip speed 
conditions which result in transonic and supersonic tip speeds. This is due 
to the sensitivity of the results to both the aerodynamic model and the need 
for accurate airfoil characteristics. 
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APPENDIX A 


COORDINATE TRANSFORMATION RELATIONSHIPS FOR 
VARIOUS QUANTITIES IN THE ANALYSIS 


Unit Vector Review 

This appendix assumes that the reader is familiar with vector properties, 
however a quick overview of the properties of unit vectors in 3-dimensions is 
given as an introduction because they are used throughout this appendix. 

Given a three-dimensional vector V, whose orthogonal unit directions 
are ij, and i^, respectively, 


V = Vjij + V 0 io + V^i 


2*2 


3 3 


then the unit vector for this vector is 


- + a i 0 + a i 0 

v Vi 1 V 2 2 V 3 3 


■ V’ 


whose direction cosines are 


V 1 


Vj/V 


v 2 /v 

V 3 /V 


where 


V = 



+ 



1/2 


such that 


V - vi v 
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The dot product of two unit vectors i x and i y is then, 


l *S ‘ W “*3^3 

The cross product of the same two vectors is then 

T X X V *1*1 + ^2*2 + hh 

where 

% * «x 2 “y 3 ‘ 

h = °Sc3 a yi ' % a y 3 
% = **1^2 " °^ K 2 < Vl 

Note that the cross product of orthogonal unit vectors is a unit vector in 
itself. 
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Input Chordwise Unit Vectors 


In the input coordinate system, (x,y,z) the input unit chordwise vector 
from figure 42 is: 


= cos0i + sin8i„ 


where 0 is the local input blade angle including built-in pitch. 

For a swept blade, a point (A) along the blade reference line (1/4 chord- 
line) is displaced in the x, y plane by an angle 6 as seen in figure 42, thus 
the input unit chordwise vector becomes 


I 



a c 


I 



z 


in the cylindrical coordinate system (r, <f>, z) where the input chordwise 
direction cosines in the cylindrical coordinate system are: 


a c 


cos 8 sin6 
cos 8 cos 6 
sin8 
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Input Thickness Unit Vector 

In the input coordinate system, the input unit thickness vector is 



-sinGiy + 


cos 0i_ 
2 


From figure 42, in a manner similar to that of the input chordwise unit 
vector, the input thickness unit vector in the cylindrical coordinate system 
is found to be: 



where the direction cosines are: 


sinG sin6 


sin6 cos 6 


cos 0 



129 



Spanwise Unit Vector 


Consider a segment of the blade reference line centered about point A 
in figure 43, whose end points are denoted by the subscripts 1 and 2, hub 
to tip respectively. The unit vector for the spanwise direction at point A, 
positive from hub to tip, can be approximated by the spanwise unit vector of 
this segment. 


T - x 2lHl T + llZ Li 
8 d x d y d z 


where 


d = [(x 2 - Xl ) 2 + (y 2 -yi ) 2 + (Z2 -Z i) 2 ] 12 


In cylindrical coordinates the spanwise unit vector and direction cosines 
become respectively 


where 


is = Vr + V* + <V Z 

* (r 2 -r 1 )cosC/d 
a s = ( r 2 +r i) sin 5/ d 

“ ( z 2~ z l ) /d 
z 


d * [ r l + r 2 “ 2r 1 r 2 cos(2^) + (z 2 -Zj) 2 ] 


1/2 


5 



(« 2“«1 ) 


Where the following approximations are made for the lag angle and radius at 
point A. 


<A - 2 

r A ’ j ( r l +r 2 1 
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Normalwise Unit Vector 


Since the input chordwise unit vector and the spanwise unit vector both 
lie in the blade element plane, the normalwise unit vector is just the unit 
vector of the cross product of the spanwise and input chordwise vector. 


i 

n 


x 




x 





The direction cosines for the normalwise unit vector in the cylindrical 
coordinate system are 


% ' 7 ( “>zV 

( VV V c i r ’ 


where 
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Chordwise Unit Vector 


In the blade element coordinate system the chordwise direction must be 
normal to the spanwise direction and the norraalwise direction to define an 
orthogonal system, thus the unit chordwise vector is 


i = i x i = a l + a l . + a 1 
c n s C r r <j> c z z 


The direction cosines in the cylindrical coordinate system are then: 


a = a a -a a 
n x s„ n„ s 4 


> _ II— O A 

z z 9 


a c = a n a « ~ a n a s 
c <|> n z s r "r s z 


W V-r 
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Blade Geometry Quantities in the Blade 
Element Coordinate System 

Because it is necessary to calculate the aerodynamic solution in the 
blade element coordinate system, the local input blade section properties must 
be transformed to the blade element system. 

The blade element chord length is obtained from the input blade chord by 
use of the dot product relationship for unit vectors. 

c ■ c iV T c ' ‘iK'V ° C + < V ac *° ei * ) 


And in a similar manner the blade element thickness is found to be 

c = + “nA + ) 

r x r ♦ 1 4> z L z 
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Velocity Vectors in the Blade 
Element Coordinate System 

The velocities in the cylindrical coordinate system must be transformed 
to the blade element system for computational purposes. In the cylindrical 
coordinate system the velocity at any point is: 


V - Vi v = v r i r + + V z i z 


the unit vector and direction cosines are then; respectively, 


i = ot i +a lx + ex i 

V v r 1 r <f v z z 


a - V /V 

v r r 


a = Vx/V 

v$ ♦ 


% - v v 


where 


V - (v* + vj + v\) 


1/2 


The blade element velocity components are then obtained by use of the dot 
product relationships as: 


Vi *i = vfa a+aa+aa] 
V1 v l s Vl V S, %% V. S, > 


r r r 


z z 


V c ■ Vl v-ic ■ V K r °c/ 

\ ' V V T n - V (\ r V \\* \\ ] 


Thus the velocity vector in the blade element coordinate system becomes 


Vs * Ve * 


Vi 
n n 


whose direction cosines are the respective bracketed quantities in the 
velocity component relationships above. 
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Force Vectors 


Once the blade forces per unit area have been found in the blade element 
coordinate systems, they must be transformed to the cylindrical coordinate 
system to obtain thrust and torque. 

The blade force per unit area and direction cosines are respectively 


h = “F^s + “fA + “F^n 

sen 

°P 8 - V F 

% ' F c /F 

■¥„ - V p 

where 

( o 2 2 ^ 

F = + F“ J 

^ s c n J 


using the dot product relationship the component forces per unit area in the 
cylindrical coordinate system are found to be 


p i F ’i r ’ F (“F c a c r * “F.’V 

P V T * ■ p (%v % v x's’ 
F V T * ‘ <*, V “F„%) 


The thrust per unit area is then just 

T = F 


z 


and the torque per unit area is just 

Q “ rF * 
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APPENDIX B 


BIOT-SAVART RELATIONSHIPS FOR THE INDUCED 
VELOCITY OF A CONSTANT STRENGTH VORTEX SEGMENT 

The relationships shown here are for a straight line constant strength 
vortex segment in cylindrical and cartesian coordinate systems which can then 
be transformed to any desired coordinate system. Consider a straight line 
vortex segment in a cylindrical coordinate system with end points (r,<f>,z) A 
and (r,<J>,z) B as shown in figure 49 with a circulation strength I\ From 
reference 15, the induced velocity components at a field point P, from the 
Biot-Savart Law, in a cartesian coordinate system coincident with the 
cylindrical coordinate system are, 

Vi x = K ^ y p" y A^ z p" z B^^p- z AKyp-y B )] 
v iy'^ K Uz p -z A )(x p -x B Hx p -x A )(z p -z B )] 

Vi z " 4 ^r K ^v x A^yp‘yB^^y P -yAKx p -x B )] 

K = AP+BP 1 

(ap) (bp) (ap) (bp)+ (x p -x A ) (x p -x B )+ (y p -y A )(y p -y B ) + (y z a^ z p" z b) 

AP - [(x p -, A ) 2 ♦ (y p -y A ) 2 ♦ ( V x a ) 2 ] I/2 

bp - [(x p -x B ) 2 + (y p -y B ) 2 ♦ ( V z b ) 2 ] 1/2 

where all o£ the distances have been nondimensionalized by the reference 
length R (propeller radius). 
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Applying the standard relationships between a cartesian and cylindrical 
coordinate system, 


x * r cos<t> 
y = r sin<J> 
z = z 

The induced velocity components, after some algebriac manipulation, in the 
cylindrical coordinate systems are, 

V ir - r*G*[(z p -z A )r B sin(4» p -^ B )-(z p -z B )r A sin(4> p -<|) A )]«*r*GC r 
Vi = r *G*[(z p -z A )r B cos(4. p -(J. B )-(z p -z B )r A cos(<fr p -«}> A )]=r«GC 4) 

V iz - r*G‘[(r p -r A )sin(<|> p -«J) A )-(r p -r B )sin(<t. p -<|) B )+r A r B 8in((). A -(J> B )]*r*GC z 

G . -1- AP + BP 
4itR (AP) (BP) 

1 

(AP ) (BP)+r p -r p r A cos (<t> p -$ A )-r p r B cos ( , l , p “ < l , B ) + r A r B cos (<|> A -<)> B )+ (z p -z A )(z p -z B ) 

AP - [ r p _ 2r p r A cos (♦ p -i)> A )+r A +(z p -z A ) 2 ] 1/2 
BP = [r 2 -2r p r B co8(4> p -^ B )+r 2 +(z p -z B ) 2 ) 1/2 

The geometric influence coefficients (GC r , GC^, GC Z ) can then be transformed 
to the blade element coordinate system for the field point under considera- 
tion. 
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Thus by the appropriate transformation the induced velocity components in the 
blade element coordinate system become 


v. * T GC 
m n 


v. * r GC 
1C c 


V. * r GC 
IS s 


The transformation required for this last step are shown in Appendix A. 


I 


i 


I 
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APPENDIX C 


a 


«I» b I 

A 


A* 


b 


c 

C 


c d» C D 

So 

Cd cu 


<T> 


CL2 


°D. 


c,,c 


1»°L 


LIST OF SYMBOLS 


linearized lift curve slope (per radian), or transformation 
variable 


Schwartz-Christof fel parameters 

cascade correlation parameter or generalized wake contraction 
limit 

Van Driest constant (26.0) 

number of blades or cascade correlation parameter for exit 
deviation angle (deg) 

airfoil chord (ft) 

Vector of aerodynamic quantities representing the constant 
vector in a system of simultaneous linear equations 

section drag coefficient 

minimum drag coefficient in the "drag bucket" 
drag due to lift 


dc D 


(t/cr 1 / 3 , 


generalized drag due to lift 


section wave drag coefficient (equal to minimum drag 
coefficient minus subsonic minimum drag coefficient) 

■ (C D ) (t/c)-5/ 3 , generalized wave drag 
P 

drag coefficient at zero angle of attack 

extrapolated minimum drag coefficient from the non- "drag 
bucket" region 

wall friction drag coefficient 

combined nonlinear correction term added to matrix constant 
section lift coefficient 
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c 


c 


m 

P 


C 


P r 


C V 


D 


e 


design lift coefficient 

design lift coefficient for a reference airfoil 

lift coefficient at zero angle of attack 

lift coefficient at zero angle of attack for a reference 
airfoil 

maximum lift coefficient 

section lift curve slope (per radian) 

dCj j /o 

* (t/c) ' , generalized lift curve slope 

lift curve slope of 2n 

function of Mach number and velocity 

o O 

specific heat pressure (ft /sec /deg R) 

pressure drag coefficient 

specific heat volume ( ftVsec /deg R) 

nonlinear correction term. added to the linearized velocity 
(fps) 

axial compressor force coefficient 

nonlinear correction term added to the linearized lift 
coefficient expression 

nonlinear correction term added to the linearized angle of 
attack expression (radians), or rotational compressor force 
coefficient 

cascade diffusion parameter, or drag force (lbf or 
lb f / ft 2 ) 

error 
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f 

F 

F v (x) 

g 


force /span (lb^ /ft) 

compressor force coefficient corrected for the absence of the 
induced velocity 

distribution of noninduced inflow velocity ratio 
blade gap (ft) 


GC 


geometric influence coefficient (per ft) 


GC* geometric influence coefficient summed over all blades for all 

segments of a particular filament (per ft)(equation 4) 

GC geometric influence coefficients differenced to relate to the 

bound circulation (per ft) 


h 


camber, maximum displacement of profile mean line from chord 
line, or duct height (ft) 


h/t 


generalized camber parameter 


unit direction vector or cascade incidence angle (degrees), 
or /-I 



j 

k 

K 


minimum loss incidence angle (degrees) 

cascade correlation parameter for minimum loss incidence angle 
(degree) 

entropy 

summation indices on propeller position or blade station 
summation indices on trailing filament number 
cascade correction factor for isolated airfoil 


^it cascade correlation parameter for minimum loss incidence angle 

K K 

l rM Q* 2 r- 2 generalized wake transport parameters for extension of the 
inboard sheet to r = 1 


K K 

l r- 2 * 2 r- j generalized wake transport parameters for extension of the 
inboard sheet to r «*1 
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cascade correlation parameters for deviation angle 

summation indices on segment number of a trailing filament 

lift per unit span (lb^/ft), or duct wall length, .or number 
of propeller positions 

summation indices on blade number 

Mach number or matrix of aerodynamic and geometric influence 
coefficients for a system of simultaneous linear equations 

critical Mach number 

cascade correlation parameter for deviation angle 

cascade correlation parameter for minimum loss incident angle, 
or normal coordinate (dimensionless) 

empirical cascade correlation parameter 

empirical cascade correlation parameter 
Reynolds number, r r o r U r /y r (dimensionless ) 
static pressure (lb^/ft 2 ), or power (ft lb^/sec) 

o 

Base Pressure (lb^/ft ) 
total pressure (lbj/ft ) 
turbulent Prandtl number (dimensionless) 
heat flux 

radial coordinate in the cylindrical coordinate system 

blade tip radius (ft) 

gas constant ( ft 2 /sec 2 /deg R) 

radial coordinates of selected generalized wake features 
nondimensionalized by the blade radius 
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st reamwise coordinate (dimensionless) 
maximum airfoil thickness (ft) 

compressor force ( lb £ ) » or temperature (deg R) , or. thrust 
(lb f ) 

total temperature (deg R) 
generalized thickness parameter 

magnitude of velocity, or velocity component (fps) 

total cascade velocity including induced velocities (fps) 

axial cascade velocity component (fps) 

rotational cascade velocity component (fps) 

universal velocity (dimensionless) 
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velocity as defined where used in text (fps) 

momentum induced velocity (fps) 

wake transport velocity (fps) 

complex coordinates in duct plane (r + iz) 
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system 
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system 
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V 

compressor momentum induced velocity (fps) 

a 

angle of attack or cascade flow angle or swirl angle to axis 
(deg) 

a e 

angle between U e and axial direction (deg) 

a i 

cascade incidence angle (deg) 

a min 

angle of attack for Q (deg) 

a 

max 

angle of attack for CL max (deg) 

a s 

cascade stagger angle (deg) 

a o 

angle of attack for zero lift (deg) 

Subscripts 


B 

due to bound vortex 

c 

blade element chordwise direction 

fr 

due to skin friction force 

H 

hub conditions 

i 

of a praticular inflow station 

i» j.P 

summation indices 

n 

blade element normalwise direction, or in the direction of the 
normal 

pr 

due to pressure force 

r 

reference conditions, based on structural sea level atmos- 
pheric conditions for all thermodynamic quantities. The 
reference r f is the inlet outer radius, and the velocity is 
the mean inlet velocity. 

s 

blade element spanwise direction, or in the streamwise direc- 
tion 

T 

tip condition 
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Tip 


blade tip 


Superscripts 


due to trailing wake 


wall conditions 


freestream value, or stagnation conditions 


inlet of the blade row 


exit of the blade row 


freestream or edge of boundary layer 


iteration number 


mean or average quantity or indicates vector 
corrected exit flow angle due to back pressure (deg) 


n 

angle of attack for U D 0 (deg) 

metal angle of blade camberline at trailing edge relative to 
the stagger angle (deg) 

flow swirl angle (deg) 

ratio of specific heats or boundary layer thickness (ft) 

circulation of vortex segment (ft 2 /sec) 

boundary layer displacement thickness (ft) 

flow deviation angle from the metal angle 

cascade correlation parameters for deviation angle (deg) 

boundary layer thickness, 6/r f (dimensionless) 

boundary layer displacement thickness, 6*/r r (dimensionless) 

empirical correction factor for lift curve slope 

= (M 2 -l) (t/c)- 2/3 , generalized Mach number 
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von Karman constant (0.41) 
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pressure ratio, p/p r 

density (slugs/ft 3 ) 
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stress^ t A w (dimensionless) 

reference coordinate in the cylindrical coordinate system 
(rad), or inflow angle (rad) 

blade dissipation function, rate of dissipation of energy of 
the propeller boundary layer 

effective camber angle (deg) 

empirical curve fit relationship 

Clauser constant (0.016) (dimensionless) or transformation 
function 

stream function (slugs/ft) 

wake age or wake azimuth position (deg or rad) - see section 
entitled "Generalized Wake Model" 

propeller rotational speed (rad/sec) 
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Table I. Predicted Propeller Performance Results - Mach Number Variation Study 
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Table II. Predicted Propeller Performance - Advance Ratio Variation Study 
M = 0.8, 100% Design Power Loading 
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Table III. Predicted Propeller Performance Results - Twist and Spinner 
Variation Study, M = 0.8, 100% Design Power Loading 
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Table Vc. Performance Map Predictions, SR-1M, Modeling Option 
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Table Via. Performance Map Predictions, SR-3, Modeling Option 
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Table VIb. Performance Map Predictions, SR-3, Modeling Option 
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Table VII. Predicted Propeller Performance Results - Wind Tunnel Study 
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Table VII. Continued 
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Table VIII. Predicted Low Speed Coaxial Performance Results 
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Table IX. Predicted Moderate Speed Coaxial Propeller Performance 
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Figure 1 — Conceptual Model of a High Speed Propeller-Nacelle Configuration 


COORDINATE IN PLANE OF ROTATION (PLAN VIEW) , y 




Figure 2. Variation of Blade Geometry of a Model Propeller with Sweep, Axial and Inplane 
Views at Zero and 60 Deg. Pitch 
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Figure 3. Simplified Structure of Major Program Tasks 
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Figure 5. Cylindrical and Cartesian Coordinate Systems for Propeller Geometry 
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Figure 6. Blade Element Coordinate System 
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Figure 7. Representation of Propeller Blade and Wake by Bound and Trailing Vortex Segments 
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LOCAL BLADE ELEMENT VELOCITY VECTOR AND FORCE DIAGRAM 

Figure 8. Aerodynamic Related Quantities in the Local Blade Element Coordinate System 
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Figure 9. Vector Oiagram of Drag Forcee in the Blade Element Coordinate Syetem 
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Figure 10. Evaard Tip Relief Region on Propeller Blade 










Figure 11. Definition of Surface Regions and Coordinate System for a Swept Wing with 
Subsonic Leading Edge and Supersonic Trailing Edge 
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Figure 12. Geometric and Velocity Quantities for the Cascade Airfoil Correlation 
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Figure 15. Classical Wake Vortex Model 
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Figure 16. Modified Classical Wake Vortex Model 
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Figure 17. Characteristics of the Generalized Wake Model 
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Figure 20. Induced Velocity Variations for Vortex Segment with Finite Core 
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MACH CONE SURFACE 
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Figure 21. Mach Cone Caused by Propeller Tip Operating at Mach Number Greater than 1.0 










PROPELLER POSITION 3 



Figure 23. Relative Blade and Wake Geometry for a Two Bladed Propeller for a Azimuth 
Increment Variation of 30 Degrees 
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Figure 24. Cylindrical and Cartesian Coordinate Systems for Coaxial Counter Rotating 
Propeller Configurations 
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Figure 25. Conformal Mapping of Duct 
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Figure 28. Flow Diagram of the Combined Propeller-Nacelle Analysis 
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Figure 32. Computed Inviscid Axial Velocity at the Blade Leading Edge and Quarter Chord 
of theSR*3 Propeller at Moo= 0.8 
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Figure 34. Comparison of Spinner Profiles 
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Figure 38 Effect of Compressible Flow Tip Models on Propeller Efficiencies for the SR*2 
Propeller, Area Ruled Spinner. 
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Figure 39 Effect of Compressible Flow Tip Loss Models on Propeller 
Efficiencies for the SR-1 Propeller, Area Ruled Spinner. 
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Figure 40 Effect of Compressible Flow Tip Loss Models on Propeller 
Efficiencies for the SR-3 Propeller, Area Ruled Spinner. 
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Figure 41. Effect of Cascade Data Correlation Models on SR-1 Propeller Efficiency, 
Fixed Blade Angle (60.2°), 0.8 Mach Number 
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Figure 42 Effect of Cascade Data Correlation Models on SR-1 Propeller 
Power Coefficient, Fixed Blade Angle (60.2°), 0.8 Mach Number 
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Figure 43 Effect of Cascade Data Correlation Models on SR*1 Propeller 
Thrust Coefficient, Fixed Blade Angle (60.2°),0.8 Mach Number 
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Figure 45 Centrifugal Effect on Blade Twist 
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Figure 46. Effect of Centrifugal Twist on Propeller Efficiency — SR-3 
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Figure 48. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-2 Propeller with Area Ruled Spinner — Mach Number Variation, 





Figure 49. Effect of Airfoil Data Source and Aerodynamic Modeling Features with Cascade 
Data on SR-2 Propeller with Area Ruled Spinner — Mach Number Variation, 
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Figure 50. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-1M Propeller with Area Ruled Spinner — Mach Number Variation, 
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Rgure 51. Effect of Airfoil Data Source and Aerodynamic Modeling Features with Cascade 
Data on SR-1M Propeller with Area Ruled Spinner — Mach Number 
Variation, J = 3.07, C D = 1.7 
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Figure 52. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-3 Propeller with Area Ruled Spinner — Mach Number Variation, 
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Figure 53. Effect of Airfoil Data Source and Aerodynamic Modeling Features with Cascade 
Data on SR-3 Propeller with Area Ruled Spinner — Mach Number 
Variation, J =3.07, C D = 1.7 
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Figure 55. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-2 Propeller with Area Ruled Spinner — Advance Ratio Variation, 

M =0.8, 100% Design Power Loading 
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Figure 56. Effect of Airfoil Data Source and Aerodynamic Modeling Features with Cascade 
Data on SR-2 Propeller with Area Ruled Spinner — Advance Ratio Variation, 

M =0.8, 100% Design Power Loading 
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Figure 57. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-1M Propeller with Area Ruled Spinner — Advance Ratio Variation, 

M * 0.8, 100% Design Power Loading 
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Figure 58. Effect of Airfoil Data Source and Aerodynamic Modeling Features with Cascade 
Data on SR-1M Propeller with Area Ruled Spinner — Advance Ratio Variation, 

M * 0.80, 100% Design Power Loading 
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Figure 59. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-3 Propeller with Area Ruled Spinner — Advance Ratio Variation, 

M * 0.8, 100% Design Power Loading 
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Figure 60. Effect of Airfoil Data Source and Aerodynamic Modeling Features with Cascade 
Data on SR*3 Propeller with Area Ruled Spinner — Advance Ratio Variation, 

M = 0.8, 100% Design Power Loading 
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Figure 62. Effect of Airfoil Data Source and Aerodynamic Modeling Features without Cascade 
Data on SR-1M Propeller with Conic Spinner — Advance Ratio Variation, 

M =0.80, 100% Design Power Loading 
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Figure 63. Effect of Airfoil Deta Source and Aerodynamic Modeling Features with Cascade 
Data on SIMM Propeller with Conic Spinner — Advance Ratio Variation, 

M *0.8, 100% Design Power Loading 
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Figure 64. Effect of Airfoil Dete Source and Aerodynamic Modeling Features without Cascade 
Data on SR-1 Propeller with Conic Spinner — Advance Ratio Variation, 

M =0.8,100% Design Power Loading 
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Figure 65. Effect of Airfoil Date Source end Aerodynamic Modeling Featuree with Cascade 
Data on SR-1 Propeller with Conic Spinner — Advance Ratio Variation, 

M * 0.8. 100% Design Power Loading 
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Figure 66a. Comparison of Measured and Predicted SR-2 Propeller Performance, 
-Modeling Option 1 
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Figure 66b. Comparison of Measured and Predicted SR-2 Propeller Performance, 
•Modeling Option 2 
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Figure 66c. Comparison of Measured and Predicted SR-2 Propeller Performance, 
•Modeling Option 3 
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Figure 67a. Comparison of Measured and Predicted SR-2 Propeller Performance, M = 0.8 
-Modeling Option 1 
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Figure 67b. Comparison of Measured and Predicted SR-2 Propeller Performance, M = 0 
•Modeling Option 2 
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Figure 67c. Comparison of Measured and Predicted SR-2 Propeller Performance, 
•Modeling Option 3 
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gure 68b. Comparison of Measured and Predicted SR-1M Propeller Performance M =0.6 
•Modeling Option 2 
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Figure 69a. Comparison of Measured and Predicted SR-1M Propeller Performance, 
•Modeling Option 1 
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Figure 69b. Comparison of Measured and Predicted SR-1M Propeller Performance, M =0 
•Modeling Option 2 
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Figure 70b. Comparison of Measured and Predicted SR*3 Propeller Performance, 
•Modeling Option 3 
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Figure 71a. Comparison of Measured and Predicted SR-3 Propeller Performance, M =0.( 
•Modeling Option 2 
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Figure 71b. Comparison of Measured and Predicted SR-3 Propeller Performance, 
•Modeling Option 3 
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Figure 72. Predicted Effect of Propeller/Tunnel SizeonSR-3 Propeller Performance 
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Figure 73. Influence of Tunnel Size on Predicted Propeller Performance for a Fixed 
Propeller/Nacelle Configuration (SR-3) 
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Figure 75. 
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Figure 76. Comparison of Predicted and Measured Thrust Coefficient for Coaxial 
Propeller — Low Speed (Ref. 72) 
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Figure 77. Continued 
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Figure 78. Spinner Planforms for Moderate Speed Coaxial Propeller Study (Ref. 73) 
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Figure 70. Comparison of Predicted Nacelle Induced Inflow Profiles for 
Type A ft B Nacelles of Ref. 73 
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Figure 80 Concluded 
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figure 81 Comparison of Predicted and Measured Power Coefficient Variation Between 
Front and Rear Propellers, Optimum a/ 3 (Test Data From Fig. 14 of Ref. 73) 
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Figure 82. Effect of Spinner Design on the Comparison of Predicted and Measured 
Propeller Performance, A/3 = 0.8°, M =0.70 (Ref. 73) 
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Figure 83. Comparison of Predicted and Measured Power Coefficient Variation Between Front 
and Rear Propellers, A/3 =0.8°, M =0.70 (Test Data From Figs. 15, 16 of Ref. 73) 
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Figure 86. Field Point and Vortex Segment in Cartesian and Cylindrical Coordinate Space 
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